Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,424,468 A→C 7.3% intergenic (+82/+48) ulaF → / ← yjfY L‑ribulose 5‑phosphate 4‑epimerase/YhcN family protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,424,4680AC 7.3% 391.5 / 15.4 117intergenic (+82/+48)ulaF/yjfYL‑ribulose 5‑phosphate 4‑epimerase/YhcN family protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (62/46);  new base C (3/6);  total (65/52)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 119 total aligned reads displayed.

CTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTC  >  NC_000913/4424323‑4424617
                                                                                                                                                  |                                                                                                                                                      
cTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑gg                                                                                                                                                    <  2:282843/149‑1 (MQ=255)
 ttCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑ggg                                                                                                                                                   >  3:82698/1‑149 (MQ=255)
 ttCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑ggg                                                                                                                                                   >  2:287490/1‑149 (MQ=255)
   ccTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑ggggg                                                                                                                                                 <  2:474956/149‑1 (MQ=255)
   ccTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATACCGCGGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCA‑ggggg                                                                                                                                                 <  4:431959/149‑1 (MQ=255)
    cTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCA‑GGGGgc                                                                                                                                                <  4:257537/149‑1 (MQ=255)
     tgatgaATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGgca                                                                                                                                               >  2:488711/1‑149 (MQ=255)
      gatgaATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAAAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTAGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCA‑GGGGgctt                                                                                                                                              <  2:172241/149‑3 (MQ=255)
       atgaATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGgcaaa                                                                                                                                             <  4:497812/149‑1 (MQ=255)
       atgaAAAAAAACTTAATGCGTAAACACGGTCCTACCGCTTATTACGGGCAAAAGTAGAACACGCGCTCCGAAAATTTCCTTCCCCGGGAGATAACTGGTCTAATTCCGCACCCGTTTTTA‑AAAAAAAAGCCCCCTTCGCA‑GGGGgcttt                                                                                                                                             <  4:24114/149‑4 (MQ=255)
           aTAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCA‑GGGGgcttttttt                                                                                                                                         <  2:410888/149‑8 (MQ=255)
           aTAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCA‑GGGGgcttttttt                                                                                                                                         <  3:217460/149‑8 (MQ=255)
            tAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGgcaaagcaaa                                                                                                                                        <  4:478548/149‑1 (MQ=255)
            tAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGgcaaagcaaa                                                                                                                                        >  4:194790/1‑149 (MQ=255)
            tAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGgcaaagcaaa                                                                                                                                        >  4:313913/1‑149 (MQ=255)
            tAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTACTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGgcaaagcaac                                                                                                                                        >  1:640625/1‑148 (MQ=255)
               acacTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTa                                                                                                                                     >  3:293362/1‑149 (MQ=255)
                 acTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATg                                                                                                                                   <  1:57126/149‑1 (MQ=255)
                   ttCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGc                                                                                                                                 >  4:152434/1‑149 (MQ=255)
                    tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCa                                                                                                                                >  1:344101/1‑149 (MQ=255)
                     cATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCaa                                                                                                                               >  3:832959/1‑149 (MQ=255)
                      aTGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAAt                                                                                                                              >  2:645057/1‑149 (MQ=255)
                      aTGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAAt                                                                                                                              >  2:429352/1‑149 (MQ=255)
                      aTGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAAt                                                                                                                              <  2:132815/149‑1 (MQ=255)
                      aTGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAAt                                                                                                                              <  1:363356/149‑1 (MQ=255)
                       tGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATg                                                                                                                             >  1:165931/1‑149 (MQ=255)
                       tGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATg                                                                                                                             >  4:423523/1‑149 (MQ=255)
                           tAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTc                                                                                                                         >  2:700590/1‑149 (MQ=255)
                           tAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTc                                                                                                                         >  3:98434/1‑149 (MQ=255)
                             aaCACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGt                                                                                                                       >  3:278505/1‑149 (MQ=255)
                              acacGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGtt                                                                                                                      <  2:696707/149‑1 (MQ=255)
                              acacGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGtt                                                                                                                      <  3:769171/149‑1 (MQ=255)
                                    tCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat                                                                                                                >  3:832574/1‑149 (MQ=255)
                                       tAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATAcc                                                                                                             >  2:831312/1‑149 (MQ=255)
                                           gCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑ggg                                                                                                                                                   >  1:582214/1‑107 (MQ=255)
                                           gCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑ggg                                                                                                                                                   <  2:582214/107‑1 (MQ=255)
                                           gCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt                                                                                                         >  2:879343/1‑149 (MQ=255)
                                           gCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt                                                                                                         >  3:581397/1‑149 (MQ=255)
                                           gCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt                                                                                                         <  1:648523/149‑1 (MQ=255)
                                             ttattaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGc                                                                                                       <  2:424635/149‑1 (MQ=255)
                                                 taCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGc                                                                                                   <  3:423523/149‑1 (MQ=255)
                                                    gggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGAt                                                                                                <  3:267408/149‑1 (MQ=255)
                                                      gCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATaa                                                                                              >  3:55082/1‑149 (MQ=255)
                                                       cAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATaaa                                                                                             <  1:145654/149‑1 (MQ=255)
                                                        aGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAat                                                                                            >  3:401943/1‑149 (MQ=255)
                                                        aGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAat                                                                                            <  1:645057/149‑1 (MQ=255)
                                                            gTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATcc                                                                                        >  2:226040/1‑149 (MQ=255)
                                                             tAGAACACGCTCTGGGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt                                                                                                                   <  4:328340/121‑1 (MQ=255)
                                                             tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt                                                                                                                   >  3:328340/1‑121 (MQ=255)
                                                             tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCcg                                                                                       >  3:122992/1‑149 (MQ=255)
                                                                   acGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTgg                                                                                 >  1:129659/1‑149 (MQ=255)
                                                                     gcgcTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTg                                                                                                           <  1:203278/121‑1 (MQ=255)
                                                                     gcgcTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTg                                                                                                           >  2:203278/1‑121 (MQ=255)
                                                                         tGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGc                                                                            >  3:214024/1‑149 (MQ=255)
                                                                            ggAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCatat                                                                        <  1:354078/149‑1 (MQ=255)
                                                                               aaTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGc                                                                     <  2:344101/149‑1 (MQ=255)
                                                                                aTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGcc                                                                    <  1:611109/149‑1 (MQ=255)
                                                                                aTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGcc                                                                    <  1:651487/149‑1 (MQ=255)
                                                                                      ttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACt                                                              >  2:433421/1‑149 (MQ=255)
                                                                                          tCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTcccc                                                          >  1:659643/1‑149 (MQ=255)
                                                                                            gggAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCgg                                                        >  3:512823/1‑149 (MQ=255)
                                                                                                 aTAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGGCGTTACTCCCCGGTTccc                                                   <  4:527845/149‑1 (MQ=255)
                                                                                                 aTAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTccc                                                   <  1:67826/149‑1 (MQ=255)
                                                                                                    aCTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCa                                                >  3:412132/1‑149 (MQ=255)
                                                                                                    aCTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAACCCCCTTCCAAGGGGGGCAAGGAAAACT‑TTGCAATGGTTCGTTGGTTAAACCCTGTGGTTACCATAAATTTTCGCCGTGGCCTTCATATTTCCCTTTATCCCCGGTTtcccga                                                >  3:613902/1‑146 (MQ=255)
                                                                                                     cTGGTCTAATTCCGCAGACGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCAt                                               <  4:224375/149‑1 (MQ=255)
                                                                                                      tGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTcc                                                    <  2:180005/143‑1 (MQ=255)
                                                                                                      tGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTcc                                                    >  1:180005/1‑143 (MQ=255)
                                                                                                      tGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATc                                              <  3:683187/149‑1 (MQ=255)
                                                                                                       ggTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCACGCATCa                                             >  3:347479/1‑149 (MQ=255)
                                                                                                         tCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGa                                           >  3:873645/1‑149 (MQ=255)
                                                                                                         tCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGa                                           >  3:206540/1‑149 (MQ=255)
                                                                                                         tCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGa                                           >  2:795230/1‑149 (MQ=255)
                                                                                                           tAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAAt                                         >  1:846078/1‑149 (MQ=255)
                                                                                                              ttCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTAc                                      <  3:618507/149‑1 (MQ=255)
                                                                                                               tCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                     <  1:844957/149‑1 (MQ=255)
                                                                                                                ccGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt                                                                                                                   >  4:657556/1‑70 (MQ=255)
                                                                                                                   cAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGgt                                 <  4:873645/149‑1 (MQ=255)
                                                                                                                      ccGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGgtagt                              >  1:7560/1‑149 (MQ=255)
                                                                                                                           tttC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGtt                         <  3:554944/149‑1 (MQ=255)
                                                                                                                             tC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTgc                       >  4:308442/1‑149 (MQ=255)
                                                                                                                              c‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTgcg                      >  1:574014/1‑149 (MQ=255)
                                                                                                                                aaaaaaaaGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGt                                                                  <  3:252091/104‑1 (MQ=255)
                                                                                                                                aaaaaaaaGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTgcgc                     >  1:572378/1‑149 (MQ=255)
                                                                                                                                 aaaaaaaGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGcc                    >  4:710770/1‑149 (MQ=255)
                                                                                                                                  aaaaaaGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCt                   >  2:199552/1‑149 (MQ=255)
                                                                                                                                     tttgCCCCCTTCGCA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgc                >  3:311887/4‑149 (MQ=255)
                                                                                                                                     tttgCCCCCTTCGCA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgc                >  2:199377/4‑149 (MQ=255)
                                                                                                                                     aaaGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgc                <  3:396492/149‑1 (MQ=255)
                                                                                                                                       aGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGt                                    >  1:384381/1‑127 (MQ=255)
                                                                                                                                          ccccTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATcc                                                                                        <  1:871896/72‑1 (MQ=255)
                                                                                                                                          ccccTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATcc                                                                                        >  2:871896/1‑72 (MQ=255)
                                                                                                                                          ccccTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGTCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCTTc           >  3:1709/1‑149 (MQ=255)
                                                                                                                                          ccccTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATc           >  4:860532/1‑149 (MQ=255)
                                                                                                                                          ccccTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATc           >  4:38418/1‑149 (MQ=255)
                                                                                                                                          ccccTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATc           >  1:416409/1‑149 (MQ=255)
                                                                                                                                          ccccTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATc           >  2:314906/1‑149 (MQ=255)
                                                                                                                                             atGCGAA‑GGGGGCAAAGCATACTATGGCAATGTTTCGTTGGTTATACCTAGTGCTAGCGATAAACATCCGCGCTGGAGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCATAACTACGTGGTAGTACGTTGCGCCTTGCGCATCTGt        <  2:640625/148‑1 (MQ=255)
                                                                                                                                                  aa‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttctt   >  1:232101/1‑149 (MQ=255)
                                                                                                                                                  aa‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGAT‑AATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc  <  4:59089/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAA‑GGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTC  >  NC_000913/4424323‑4424617

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑