Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,041 T→G 6.1% intergenic (+72/+258) yiiF → / ← fdhE protein YiiF/formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,0410TG 6.1% 855.9 / 37.0 262intergenic (+72/+258)yiiF/fdhEprotein YiiF/formate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base T (109/137);  new base G (8/8);  total (117/145)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 276 total aligned reads displayed.

ACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCAAAC  >  NC_000913/4079892‑4080188
                                                                                                                                                     |                                                                                                                                                   
aCCAGTGTTCCTGGCATCTGGTAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                     <  1:235207/150‑1 (MQ=255)
 ccAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                    <  2:1115302/150‑1 (MQ=255)
    gtgtTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                                 >  1:1186812/1‑150 (MQ=255)
      gtTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                               >  2:637715/1‑150 (MQ=255)
      gtTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                                <  2:409450/149‑1 (MQ=38)
       ttCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                              >  2:1692397/1‑150 (MQ=255)
       ttCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                              <  2:1336907/150‑1 (MQ=255)
        tCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                             <  2:1558671/150‑1 (MQ=255)
         ccTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                            >  1:1990545/1‑150 (MQ=255)
          cTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                           >  1:909413/1‑150 (MQ=255)
          cTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                           >  2:1360331/1‑150 (MQ=255)
           tGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcatg                                                                                                                                          <  1:838509/150‑4 (MQ=38)
               aTCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaa                                                                                                                                      >  1:30405/1‑150 (MQ=255)
                tCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAg                                                                                                                                     >  2:1305774/1‑150 (MQ=255)
                   ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCg                                                                                                                                  <  2:1740826/150‑1 (MQ=255)
                   ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCg                                                                                                                                  <  1:168788/150‑1 (MQ=255)
                   ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCg                                                                                                                                  <  2:1157875/150‑1 (MQ=255)
                     cAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGta                                                                                                                                <  2:281960/150‑1 (MQ=255)
                       aGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGtata                                                                                                                              >  1:1467751/1‑150 (MQ=255)
                         ggTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACAAAGTCGTATAAAt                                                                                                                           <  2:1773723/150‑1 (MQ=255)
                             tgAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCa                                                                                                                        <  2:440025/150‑1 (MQ=255)
                              gAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCaa                                                                                                                       >  1:1312610/1‑150 (MQ=255)
                               aGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAAt                                                                                                                      <  1:1003570/150‑1 (MQ=255)
                                ggAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATg                                                                                                                     >  2:338355/1‑150 (MQ=255)
                                 gAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGt                                                                                                                    >  1:1563678/1‑150 (MQ=255)
                                    gATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTAtt                                                                                                                 <  2:676804/150‑1 (MQ=255)
                                        gtgtCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTg                                                                                                                >  1:136064/1‑147 (MQ=255)
                                        gtgtCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAg                                                                                                             <  2:1036844/150‑1 (MQ=255)
                                          gtCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGtt                                                                                                           <  2:139555/150‑1 (MQ=255)
                                           tCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGttt                                                                                                          >  1:156846/1‑150 (MQ=255)
                                            cGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTa                                                                                                         >  1:1144099/1‑150 (MQ=255)
                                                 atCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGt                                                                                                    <  1:1984252/150‑1 (MQ=255)
                                                    aGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGG                                                                                                 <  1:1419985‑M1/150‑1 (MQ=255)
                                                    aGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGG                                                                                                 >  1:645225‑M1/1‑150 (MQ=255)
                                                     gCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGC                                                                                                <  1:1830540‑M1/150‑1 (MQ=255)
                                                     gCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGC                                                                                                >  1:700843‑M1/1‑150 (MQ=255)
                                                       gTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCT                                                                                              <  2:589558‑M1/150‑1 (MQ=255)
                                                       gTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCT                                                                                              <  2:484949‑M1/150‑1 (MQ=255)
                                                        tAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTG                                                                                             <  1:1148101‑M1/150‑1 (MQ=255)
                                                          aGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGAT                                                                                           <  2:1357172‑M1/150‑1 (MQ=255)
                                                           gCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATa                                                                                          <  2:657460/150‑1 (MQ=255)
                                                            cTGCGCGAGGAATGGTAATGGTAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTg                                                                                                      <  1:1298762/137‑1 (MQ=255)
                                                             tGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAag                                                                                        <  2:1545336/150‑1 (MQ=255)
                                                              gcgcGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagc                                                                                       <  1:825421/150‑1 (MQ=255)
                                                              gcgcGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagc                                                                                       <  1:1019702/150‑1 (MQ=255)
                                                                  gAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcag                                                                                   <  1:881311/150‑1 (MQ=255)
                                                                   aGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagc                                                                                  >  1:1109833/1‑150 (MQ=255)
                                                                    ggAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagcg                                                                                 >  2:821841/1‑150 (MQ=255)
                                                                      aaTGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgctta                                                                               <  1:597804‑M1/150‑13 (MQ=255)
                                                                      aaTGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCa                                                                               <  2:987925‑M2/150‑2 (MQ=255)
                                                                      aaTGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCa                                                                               <  1:1960330‑M2/150‑2 (MQ=255)
                                                                          gTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCatcag                                                                           <  2:1205569‑M2/150‑6 (MQ=255)
                                                                          gTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATATGTAGGCCTGATAAGCGCAGCGCatcag                                                                           <  2:1828558‑M2/150‑6 (MQ=255)
                                                                          gTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTAATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCatcag                                                                           <  1:542835‑M2/150‑6 (MQ=255)
                                                                            aaTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCatcaggc                                                                         >  1:1333398‑M2/1‑143 (MQ=255)
                                                                            aaTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGAATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCatcaggc                                                                         <  1:1922938‑M2/150‑8 (MQ=255)
                                                                             aTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgcttatcaggcc                                                                        <  1:1237132‑M1/150‑20 (MQ=255)
                                                                             aTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgcttatcaggcc                                                                        <  2:1933739‑M1/150‑20 (MQ=255)
                                                                             aTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgcttatcaggcc                                                                        <  2:795726‑M1/150‑20 (MQ=255)
                                                                                ggAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcag                                                                                   <  1:783241/136‑1 (MQ=255)
                                                                                ggAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATc                                                                     >  2:1234060/1‑150 (MQ=255)
                                                                                 gAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCt                                                                    >  2:105238/1‑150 (MQ=255)
                                                                                 gAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCt                                                                    <  2:909413/150‑1 (MQ=255)
                                                                                 gAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCt                                                                    <  2:960418/150‑1 (MQ=255)
                                                                                  aaaaaaTGATGTCAATCAGATTGCTGAAAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAATCGCAGCGCATCAGGCAATCtt                                                                   >  1:727193/1‑150 (MQ=255)
                                                                                    aaaaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTAc                                                                 <  1:637715/150‑1 (MQ=255)
                                                                                      aaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACgt                                                               <  1:993238/150‑2 (MQ=255)
                                                                                          atgTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGt                                                           <  1:1931763/150‑1 (MQ=255)
                                                                                           tgTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtc                                                          <  2:1959613/150‑1 (MQ=255)
                                                                                            gTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgcttatcaggcctaccaataaactgca                                                         <  1:83202‑M1/150‑35 (MQ=255)
                                                                                             tCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgcttatcaggcctaccaataaactgcaa                                                        <  2:828963‑M1/150‑36 (MQ=255)
                                                                                              cAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtcatc                                                       >  1:1558590/1‑150 (MQ=255)
                                                                                              cAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtcatc                                                       >  2:858798/1‑150 (MQ=255)
                                                                                               aaTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtcatca                                                      >  2:757863/1‑150 (MQ=255)
                                                                                                aTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAg                                                     <  1:984786/150‑1 (MQ=255)
                                                                                                          gaaCAACGCGCACGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagcg                                                                                 <  1:38977/110‑1 (MQ=255)
                                                                                                                   gcgcGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTa                                  >  2:1295335/1‑150 (MQ=255)
                                                                                                                   gcgcGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTa                                  <  1:351635/150‑1 (MQ=255)
                                                                                                                    cgcgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTAt                                 >  2:1683528/1‑150 (MQ=255)
                                                                                                                      cgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTg                               >  1:1265347/1‑150 (MQ=255)
                                                                                                                      cgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTg                               >  1:492914/1‑150 (MQ=255)
                                                                                                                        ttgttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGcc                                            >  1:74691/1‑135 (MQ=255)
                                                                                                                          gttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGAtgt                           >  2:1797614/1‑150 (MQ=255)
                                                                                                                            gtaggCCGGATAAGGCGTTCACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtgg                         >  2:733627/5‑150 (MQ=37)
                                                                                                                            gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtgg                         >  1:1066070/5‑150 (MQ=38)
                                                                                                                             cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgcttatcaggcctaccaataaactgcaatacattgaatttatacgactttgtaggccgga                        <  2:1661661‑M1/150‑68 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggt                        <  1:1733435/147‑1 (MQ=38)
                                                                                                                                gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtggt                     >  1:68235/1‑150 (MQ=255)
                                                                                                                                 ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtggtg                    <  2:1905861/150‑1 (MQ=255)
                                                                                                                                   ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTggg                  <  1:338355/150‑1 (MQ=255)
                                                                                                                                       gcggcgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTaaa              <  2:1990545/150‑1 (MQ=255)
                                                                                                                                         ggcgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAAtt            <  1:439125/150‑1 (MQ=255)
                                                                                                                                          gcgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTg           >  2:1727834/1‑150 (MQ=255)
                                                                                                                                           cgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGa          <  1:1797614/150‑1 (MQ=255)
                                                                                                                                           cgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGa          >  2:146449/1‑150 (MQ=255)
                                                                                                                                               aaCGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGcc      >  1:501217/1‑150 (MQ=255)
                                                                                                                                               aaCGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGcc      >  2:1091725/1‑150 (MQ=255)
                                                                                                                                               aaCGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGcc      <  1:858798/150‑1 (MQ=255)
                                                                                                                                                aCGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCa     <  2:1568055/150‑1 (MQ=255)
                                                                                                                                                   ccTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCAAAc  >  2:534338/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCAAAC  >  NC_000913/4079892‑4080188

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: