Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,035 A→C 6.5% intergenic (+66/+264) 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,0350AC 6.5% 872.8 / 40.0 274intergenic (+66/+264)yiiF/fdhEprotein YiiF/formate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base A (112/144);  new base C (10/8);  total (122/152)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

GCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCC  >  NC_000913/4079886‑4080184
                                                                                                                                                     |                                                                                                                                                     
gCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGa                                                                                                                                                       >  1:1115302/1‑150 (MQ=255)
gCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGa                                                                                                                                                       >  1:167514/1‑150 (MQ=255)
 cAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                                                      >  2:313997/1‑150 (MQ=255)
  aaGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                                                      <  1:470177/149‑1 (MQ=255)
      aCCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                 >  1:826124/1‑150 (MQ=255)
         aGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAACATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                              >  1:965425/1‑150 (MQ=255)
          gtgtTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                             >  1:1186812/1‑150 (MQ=255)
             ttCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                          <  1:899279/150‑1 (MQ=255)
             ttCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                          <  2:1336907/150‑1 (MQ=255)
             ttCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                          <  1:1931085/150‑1 (MQ=39)
              tCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                         <  1:1969589/150‑1 (MQ=255)
                cTGGCATCTGGCAAGGGGGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                       <  2:856606/150‑1 (MQ=255)
                 tGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                      <  1:723486/150‑1 (MQ=255)
                  ggCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                     >  1:589558/1‑150 (MQ=255)
                   gCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                                                                                                    <  2:1779434/150‑1 (MQ=255)
                    cATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaa                                                                                                                                   >  2:376561/1‑150 (MQ=255)
                      tCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAg                                                                                                                                 >  2:1305774/1‑150 (MQ=255)
                       cTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                 <  1:1987394/133‑1 (MQ=255)
                        tGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTc                                                                                                                               >  1:1959613/1‑150 (MQ=255)
                         ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                 <  2:1903265/131‑1 (MQ=255)
                         ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                 >  1:1903265/1‑131 (MQ=255)
                         ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCg                                                                                                                              <  2:1929167/150‑1 (MQ=255)
                          gCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGt                                                                                                                             >  1:795726/1‑150 (MQ=255)
                                  gtgAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTc                                                                                                                     <  2:157072/150‑1 (MQ=255)
                                  gtgAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTc                                                                                                                     >  2:181613/1‑150 (MQ=255)
                                    gAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCaa                                                                                                                   >  2:476332/1‑150 (MQ=255)
                                    gAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCaa                                                                                                                   >  1:1312610/1‑150 (MQ=255)
                                    gAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCaa                                                                                                                   >  2:279666/1‑150 (MQ=255)
                                     aGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAAt                                                                                                                  <  1:1003570/150‑1 (MQ=255)
                                     aGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAAt                                                                                                                  <  2:550412/150‑1 (MQ=255)
                                        aaGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAAGTCGTATAAATTCAATGTa                                                                                                               <  2:1487367/150‑1 (MQ=255)
                                           aTGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTg                                                                                                            >  2:439125/1‑150 (MQ=255)
                                            tGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGc                                                                                                           <  1:1329368/150‑1 (MQ=255)
                                              gtgtCGAATATCAGCGTAAGCTGCGCGAGGATTGGTAATGGGAAAAATTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGGAg                                                                                                         <  2:1046379/150‑1 (MQ=255)
                                              gtgtCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTg                                                                                                            <  2:136064/147‑1 (MQ=255)
                                              gtgtCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAg                                                                                                         <  2:1036844/150‑1 (MQ=255)
                                                gtCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGtt                                                                                                       <  2:1104014/150‑1 (MQ=255)
                                                gtCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGtt                                                                                                       >  2:1717251/1‑150 (MQ=255)
                                                 tCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGttt                                                                                                      >  1:156846/1‑150 (MQ=255)
                                                  cGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTa                                                                                                     >  1:1144099/1‑150 (MQ=255)
                                                     atatCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTg                                                                                                  <  1:261716/150‑1 (MQ=255)
                                                      tatCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTgg                                                                                                 <  1:1019461/150‑1 (MQ=255)
                                                      tatCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTgg                                                                                                 <  1:696202/150‑1 (MQ=255)
                                                          aGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGG                                                                                             >  1:645225‑M1/1‑150 (MQ=255)
                                                              tAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTG                                                                                         <  1:1695181‑M1/150‑1 (MQ=255)
                                                                aGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGAT                                                                                       <  2:1357172‑M1/150‑1 (MQ=255)
                                                                aGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGAT                                                                                       >  1:1166181‑M1/1‑150 (MQ=255)
                                                                aGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGAT                                                                                       <  2:1199862‑M1/150‑1 (MQ=255)
                                                                aGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATGCAGTTTATTGGTAGGCCTGATa                                                                                      <  1:629152/150‑1 (MQ=255)
                                                                 gCTGGGGGAGGAATGGTTATGGGAAAAAATGATGTGAATCAGATTGCTGTCAACGTGCGCGTTGTTCATGCGGGATGCTGCGTGAACGCCTTATCCGGCCTAGAAAGTCGTATAAATTCTATGTATTGCAGTTTATTGGTAGGGCTGATa                                                                                      <  2:1744590/150‑1 (MQ=25)
                                                                  cTGCGCGAGGAATGGTAATGGTAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTg                                                                                                  <  1:1298762/137‑1 (MQ=255)
                                                                  cTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTg                                                                                                  >  2:1298762/1‑137 (MQ=255)
                                                                   tGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAag                                                                                    <  2:1545336/150‑1 (MQ=255)
                                                                       cgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgca                                                                                >  1:1075774/1‑150 (MQ=255)
                                                                       cgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgca                                                                                >  1:864430/1‑150 (MQ=255)
                                                                        gAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcag                                                                               <  1:313997/150‑1 (MQ=255)
                                                                            aaTGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCa                                                                           <  2:987925‑M2/150‑2 (MQ=255)
                                                                             aTGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCat                                                                          <  1:527316‑M2/150‑3 (MQ=255)
                                                                               ggTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCatca                                                                        >  2:1147240‑M2/1‑146 (MQ=255)
                                                                                gTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCatcag                                                                       <  2:1076794‑M2/150‑6 (MQ=255)
                                                                                 tAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCatcagg                                                                      >  2:805983‑M2/1‑144 (MQ=255)
                                                                                 tAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCatcagg                                                                      >  1:1813249‑M2/1‑144 (MQ=255)
                                                                                  aaTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCatcaggc                                                                     >  1:143784‑M2/1‑143 (MQ=255)
                                                                                   aTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgcttatcaggcc                                                                    <  2:795726‑M1/150‑20 (MQ=255)
                                                                                   aTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgcttatcaggcc                                                                    <  1:1237132‑M1/150‑20 (MQ=255)
                                                                                   aTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                    <  2:167514/150‑1 (MQ=255)
                                                                                     gggAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAAt                                                                  >  2:727664/1‑150 (MQ=255)
                                                                                      ggAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATc                                                                 >  2:1234060/1‑150 (MQ=255)
                                                                                      ggAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCACGCCGCATCCGGCATGa                                                                                                                                                       <  2:1015560/64‑1 (MQ=18)
                                                                                       gAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCt                                                                <  2:674644/150‑1 (MQ=255)
                                                                                       gAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCt                                                                <  2:909413/150‑1 (MQ=255)
                                                                                       gAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCt                                                                <  2:960418/150‑1 (MQ=255)
                                                                                         aaaaaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTa                                                              <  1:1384521/150‑1 (MQ=255)
                                                                                           aaaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgcttatcagg                                                                      <  2:1562541‑M1/140‑18 (MQ=255)
                                                                                           aaaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATgcgctgcgcttatcagg                                                                      >  1:1562541‑M1/1‑123 (MQ=255)
                                                                                             atgatgTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCTTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCt                                                          <  2:622358/150‑1 (MQ=255)
                                                                                             atgatgTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCt                                                          <  2:1443749/150‑1 (MQ=255)
                                                                                                    cAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtcatc                                                   >  2:1957453/1‑150 (MQ=255)
                                                                                                     aaTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtcatca                                                  >  2:757863/1‑150 (MQ=255)
                                                                                                            ttGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCAt                                           >  2:1687451/1‑150 (MQ=255)
                                                                                                                tGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGccc                                       >  1:1772718/1‑150 (MQ=255)
                                                                                                                gaaCAACGCGCACGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagcg                                                                             <  1:38977/110‑1 (MQ=255)
                                                                                                                  aCAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCAt                                     <  1:1360331/150‑1 (MQ=255)
                                                                                                                    aaCGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATaa                                   <  2:826124/150‑1 (MQ=255)
                                                                                                                      cGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAAcc                                 <  1:757863/150‑1 (MQ=255)
                                                                                                                         gcgcGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTa                              <  1:351635/150‑1 (MQ=255)
                                                                                                                         gcgcGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCAATAACCTTa                              >  2:498630/1‑150 (MQ=255)
                                                                                                                            cgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTg                           >  1:1265347/1‑150 (MQ=255)
                                                                                                                            cgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTg                           >  1:500038/1‑150 (MQ=255)
                                                                                                                                     ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggt                    >  2:1733435/2‑148 (MQ=38)
                                                                                                                                      gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtggt                 >  1:68235/1‑150 (MQ=255)
                                                                                                                                      gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtggt                 >  1:83778/1‑150 (MQ=38)
                                                                                                                                        cGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGtggtggtgg               >  1:1051632/1‑150 (MQ=38)
                                                                                                                                         ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTggg              <  1:207559/150‑1 (MQ=255)
                                                                                                                                             gcggcgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTaaa          <  2:1990545/150‑1 (MQ=255)
                                                                                                                                              cggcgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAAt         >  1:1516413/1‑150 (MQ=255)
                                                                                                                                               ggcgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAAtt        <  1:439125/150‑1 (MQ=255)
                                                                                                                                                gcgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTg       >  1:840479/1‑150 (MQ=255)
                                                                                                                                                 cgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGa      >  2:652472/1‑150 (MQ=255)
                                                                                                                                                    gAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGc   <  1:987925/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCC  >  NC_000913/4079886‑4080184

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: