Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,278,959 C→T L335L (CTG→TTG)  narK → nitrate:nitrite antiporter NarK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,278,9590CT98.4% 212.3 / NA 61L335L (CTG→TTG) narKnitrate:nitrite antiporter NarK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/31);  minor base G (0/1);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGGGTCATCTGTGGATTATGAGCCTGCTGTATCTGGCAACCTTCGGCTCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTGGGTTTCATCTCTGCGATTG  >  NC_000913/1278714‑1279196
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
ggggTCATCTGTGGATTATGAGCCTGCTGTATCTGGCAACCTTCGGCTCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTc                                                                                                                                                                                                                                          >  1:456062/1‑251 (MQ=255)
  ggTCATCTGTGGAGTATGAGCCTGCTGTATCTGGCAACCTTCGGCTCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGGCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCt                                                                                                                                                                                                                                        <  2:28894/251‑1 (MQ=255)
               ccccacacgacgcccctccgatcaCCTTCGGCTCCTTCAACCGCTTCTCCGCGGGCTTTGCAAAGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGGAATACGCTTTCGTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGGTTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGGTGTTCCTGACCTTACCGACt                                                                                                                                                                                                                           <  2:307620/228‑1 (MQ=255)
                tATGAGCCTGCTGTAGCTGGCAACTGTCGGCGCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTTCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTTCATTATCTGACCGGCTGGGCTGAACTCGTGTCACGCTGTTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTGCCTGACCTTGCCGACt                                                                                                                                                                                                                           <  2:357609/250‑1 (MQ=255)
                 aTGAGCCTGCTGTATCTTGCACCCTTCGGCTCCTTCATCGGCGTCTCCGCGGGCGTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCGGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGTTGGCGCGTTCTGCAGGTGGTGCAGTATCTGACCGTCTGGTCGGAACGCGTGTCACGCTGGTGAACTTTATTCGTATGGCGAGTTTCAGCGGCTTGTTGTGCCTGACCTTGCCGACTg                                                                                                                                                                                                                          <  2:267459/250‑1 (MQ=255)
                      ccTGCTGTATCTGGCAACCTTCGGCTCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGTGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGc                                                                                                                                                                                                                    <  1:526900/251‑1 (MQ=255)
                      ccTGCTGTATCTGGCAACCTTCGGCTCCTTCATCGGCTTCTCCCCGGGCTTTGCAATGCTGTTAACAACGCAGTTCCCGGATGTTCCGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGCCGGGc                                                                                                                                                                                                                    >  2:231044/1‑251 (MQ=255)
                       ctgctgTATCTGGCAACCTTCGGCTCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAACACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGcc                                                                                                                                                                                                                   >  2:317066/1‑250 (MQ=255)
                                  ggCAACCTTCGGCTCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAg                                                                                                                                                                                                          <  2:181629/249‑1 (MQ=255)
                                      accgtCGGCTCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTTACCGTCTGGGCGGAACGCGTGTCACGCTGGTGAACTTGATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCAt                                                                                                                                                                                                    <  1:40043/247‑1 (MQ=255)
                                       ccTTCGGCTCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTCCCGACt                                                                                                                                                                                                                           >  1:307620/1‑227 (MQ=255)
                                        cTTCGGCTCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATgg                                                                                                                                                                                                  >  1:6123/1‑251 (MQ=255)
                                                  ttCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTcgc                                                                                                                                                                                        >  1:536066/1‑251 (MQ=255)
                                                       cGGCTTCTCCCCGTGCTTGGCAATGCTGTCAAAAACGCAGTTCCCAGATTTTCAGAGTCTGCAATACGCTTGCTTCGGGCCGTTTATTAGTGCGCTGTCGCTTGCTGCAGGTGGTGCATTATCTGTCCGTCTGGGCGGATCGCGGGTCACGCTTGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGGTCCTAATCTTTCCTTCTTGCGTGCAGGGCAAAAACTTCATGGCGTTCTTCGCGGTCt                                                                                                                                                                                   <  1:230243/251‑1 (MQ=255)
                                                        ggCTTCTCCGCGGGCTTTGCAATGCTGTCAAACACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCtt                                                                                                                                                                                  >  2:238236/1‑251 (MQ=255)
                                                        ggCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGTCGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCtt                                                                                                                                                                                  <  1:417407/251‑1 (MQ=255)
                                                        ggCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCtt                                                                                                                                                                                  >  2:190198/1‑251 (MQ=255)
                                                        ggCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCACGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCtt                                                                                                                                                                                  >  1:433776/1‑251 (MQ=255)
                                                         gCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTc                                                                                                                                                                                 <  2:480791/251‑1 (MQ=255)
                                                               ccGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATGCGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTg                                                                                                                                                                              <  2:291009/248‑1 (MQ=255)
                                                               ccGCGGGCTTTGCAATGCCGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGgcg                                                                                                                                                                           <  1:9164/251‑1 (MQ=255)
                                                                      cTTTGCAATGCTGTCAAAAACGCAGTTCCAGGCTGTTCAGCTTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGACGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTcc                                                                                                                                                                    >  2:98460/1‑251 (MQ=255)
                                                                                  gTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAgctgggc                                                                                                                                                        >  2:438359/1‑251 (MQ=255)
                                                                                     aaaaaCGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTCTTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAgctgggctgg                                                                                                                                                     >  1:182332/1‑251 (MQ=255)
                                                                                          cGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGt                                                                                                                                                >  2:270944/1‑251 (MQ=255)
                                                                                             aGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGt                                                                                                                                             <  1:438359/251‑1 (MQ=255)
                                                                                                  ccGGATGTTCAGATTCTGCATGGCGCTTTCTTCGGGCCGTTTATTGGGGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCAc                                                                                                                                        <  2:383268/251‑1 (MQ=255)
                                                                                                   cGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACt                                                                                                                                       <  1:190198/251‑1 (MQ=255)
                                                                                                    ggatgttacgaatcggccaTACGCTTTCTGGGGGCGGTTTATTGGGGCGGTGGCGCGTGCTGCAGGTGGTGCATTATCTGACCTTCTGGGCGGAACTCGTGTGTCGCTGGTGAGCTTTATTCTGATGGCTGTTTTCAGCGGCCTGTTGTTCCTGGCCTTGCCGTCTGGCGTGCGTGGCGTAAGCTTCATGGGGTTCTTCGCGGTCTTCCTGTCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCCCtt                                                                                                                                      <  2:307133/233‑1 (MQ=255)
                                                                                                        gTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCTCTCGAGCGTTCTGCAGGTGGTGCCTTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGCAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGAAAGCTGGGCTGGGTAGTGGTTCCACTTTca                                                                                                                                   >  2:534215/1‑249 (MQ=255)
                                                                                                          tCAGATTCTGCCATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGTGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGCCTGACGGGCAGGGCGTAAGATTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTCCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTcc                                                                                                                                   >  2:257006/1‑248 (MQ=255)
                                                                                                              aTTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAgatgct                                                                                                                            >  1:61586/1‑249 (MQ=255)
                                                                                                                      aTACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTgt                                                                                                                    >  2:5681/1‑250 (MQ=255)
                                                                                                                                 tcGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGt                                                                                                            <  1:521054/248‑1 (MQ=255)
                                                                                                                                 tcGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGt                                                                                                            <  2:6123/248‑1 (MQ=255)
                                                                                                                                       cGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCCTTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGAc                                                                                                    >  1:68399/1‑250 (MQ=255)
                                                                                                                                                       ggCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGTTTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGaaa                                                                                    <  1:466581/250‑1 (MQ=255)
                                                                                                                                                            gTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGaa                                                                              <  1:317066/251‑1 (MQ=255)
                                                                                                                                                               cTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACg                                                                                                                                                                                                                       <  2:85076/111‑1 (MQ=255)
                                                                                                                                                               cTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACg                                                                                                                                                                                                                       >  1:85076/1‑111 (MQ=255)
                                                                                                                                                                    ggtggtGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTc                                                                      >  1:227648/1‑251 (MQ=255)
                                                                                                                                                                             ttATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACg                                                             >  2:39190/1‑251 (MQ=255)
                                                                                                                                                                              tATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGt                                                            >  1:362939/1‑251 (MQ=255)
                                                                                                                                                                                   gACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCg                                                         <  1:59139/249‑1 (MQ=255)
                                                                                                                                                                                   gACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGAGTTCAGTGATCTTCCGTAAACTGACAATGGAACGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGc                                                          <  2:68399/248‑1 (MQ=255)
                                                                                                                                                                                         cTGGGCGGAACTCGTGTAACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTCGTGGTTCCACTTTCACGATGATTTCAGTGATCTTCCGTAACCTGACAATGGATAGCGTGTAAGCAGAAGGGGGGTCTGAAGAACGTGAGATGCGTg                                                  >  2:517695/1‑250 (MQ=255)
                                                                                                                                                                                           gggCGGAACTCGCGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGCCCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGCAGTGGTTCCACTTTCCAGATGATTACATTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGa                                                 >  2:202566/1‑249 (MQ=255)
                                                                                                                                                                                                     cGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAAcc                                       <  1:238236/249‑1 (MQ=255)
                                                                                                                                                                                                           aCGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTCCCTGGCGCTGTTCCTGACAGCTGGGCTGTGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACAcggc                               <  1:98460/251‑1 (MQ=255)
                                                                                                                                                                                                                ggTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACAcggcggcgg                          >  1:64599/1‑251 (MQ=255)
                                                                                                                                                                                                                ggTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACAGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACACTGGATCGCGTGa                                                                                      >  1:80433/1‑191 (MQ=255)
                                                                                                                                                                                                                     aCTTTATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCCGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTg                     >  2:502342/1‑251 (MQ=255)
                                                                                                                                                                                                                        ttATTCGGATGGCGATTTTCAGCGGCCGGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAATCTTCATGGCGTTCTTCGCGGTTTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGGGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCGGACGAACGTGCGATGCGTGAAGCGTCAACCTACACGGCGGCGGCGCTggg                   <  2:61586/250‑1 (MQ=255)
                                                                                                                                                                                                                         tATTCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGGTAGCTTCCGTAAACGGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTGGGt                  <  1:534215/250‑1 (MQ=255)
                                                                                                                                                                                                                           ttCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTggggttg               >  2:399830/1‑247 (MQ=255)
                                                                                                                                                                                                                           ttCTGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTGGGttt                <  1:399830/250‑1 (MQ=255)
                                                                                                                                                                                                                              tGATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTGGGTTTCAtc            >  1:280862/1‑251 (MQ=255)
                                                                                                                                                                                                                               gATGGCGATTTTCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTGGGTTTCAt             <  1:175143/249‑1 (MQ=255)
                                                                                                                                                                                                                               gATGGCGATTTGCAGCGGCCTGTTGTTCCTGACCTTACCGTCTGACGTGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTGTCCAGATGATTTCAGTGATCTTCCGTAAGCTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTGGGTTTCAt             <  1:40649/249‑1 (MQ=255)
                                                                                                                                                                                                                                    cGATTTTCAGCGGCCTGTTGTTCCTGATCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTGGGTTTCATCTCTGc       <  1:11879/250‑1 (MQ=255)
                                                                                                                                                                                                                                         ttCAGCGGCCTGTTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTGGGTTTCATCTCTGCGATTg  <  1:39190/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
GGGGTCATCTGTGGATTATGAGCCTGCTGTATCTGGCAACCTTCGGCTCCTTCATCGGCTTCTCCGCGGGCTTTGCAATGCTGTCAAAAACGCAGTTCCCGGATGTTCAGATTCTGCAATACGCTTTCTTCGGGCCGTTTATTGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAAGGGGGTTCTGACGAACGTGCGATGCGTGAAGCGGCAACCGACACGGCGGCGGCGCTGGGTTTCATCTCTGCGATTG  >  NC_000913/1278714‑1279196

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: