Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,496,485 G→A E825E (GAG→GAA nirB → nitrite reductase (NADH) large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,496,4850GA98.5% 228.8 / NA 67E825E (GAG→GAAnirBnitrite reductase (NADH) large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/34);  minor base T (0/1);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGACAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGG  >  NC_000913/3496239‑3496731
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
tGCCGACAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAAcgcg                                                                                                                                                                                                                                                    >  1:192644/1‑251 (MQ=255)
 gCCGACAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAacgc                                                                                                                                                                                                                                                     <  1:67587/249‑1 (MQ=255)
      cAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGc                                                                                                                                                                                                                                              >  1:285014/1‑251 (MQ=255)
       aaaCTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCATTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACGCTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCa                                                                                                                                                                                                                                             <  1:247723/251‑1 (MQ=255)
            gACGCGTACCGCACCGTGGTTAGAAAGCCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGTGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGt                                                                                                                                                                                                                                         <  2:242834/250‑1 (MQ=255)
                 gTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCGCCGTCCGGCa                                                                                                                                                                                                                                   <  2:408626/251‑1 (MQ=255)
                      gCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGcc                                                                                                                                                                                                                              <  2:153397/251‑1 (MQ=255)
                           gTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATg                                                                                                                                                                                                                         >  2:408373/1‑251 (MQ=255)
                                       ccTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTGCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAg                                                                                                                                                                                                             <  2:138880/251‑1 (MQ=255)
                                             aGGCGGGATCGATTACCTGAAAGCAGTGATCAGTGACGACAAACTGGTGCTGAACGCACATCTTGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCCATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCATATCGCGAACAGGGCCTTCCGTCAACGCCGTATGAACGTATCCCTGTAActc                                                                                                                                                                                                       <  2:54055/251‑1 (MQ=255)
                                              ggcggcATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAActct                                                                                                                                                                                                      <  2:206796/251‑1 (MQ=255)
                                                cggcATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTCACTCtgg                                                                                                                                                                                                    >  1:268655/1‑251 (MQ=255)
                                                 ggcATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCtggg                                                                                                                                                                                                   >  1:323524/1‑250 (MQ=255)
                                                     tCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCtggtgg                                                                                                                                                                                                 >  2:154965/1‑249 (MQ=255)
                                                           cccTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGAATGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGGCCGGCAACTCCGTATGAACGCATCCCAGTAACTCTGGTGGAGGAcccc                                                                                                                                                                                         <  2:238886/250‑4 (MQ=255)
                                                           aCCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAAc                                                                                                                                                                                         >  1:81140/1‑251 (MQ=255)
                                                                    cAGTGATCATTGACGACAAACTGGGGCTGAACGCCCATCTGGAAGAAGAGCTGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGAGGAGGACAACGCATGAGcc                                                                                                                                                                                >  2:70741/1‑251 (MQ=255)
                                                                     agttaTCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCa                                                                                                                                                                               <  1:228509/247‑1 (MQ=255)
                                                                     aGTGATCATTGACGACAAACTGGGGCTGAACGCACATCGGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCTCAGACGCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAGCGCTGGTGGAGGACAACGCATGAGCCa                                                                                                                                                                               <  2:412444/251‑1 (MQ=255)
                                                                     aGTGATCATTGACGACAAACGGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCa                                                                                                                                                                               <  2:390499/251‑1 (MQ=255)
                                                                        gATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAAATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTg                                                                                                                                                                            >  1:133769/1‑251 (MQ=255)
                                                                                acgacAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACa                                                                                                                                                                    <  2:257442/251‑1 (MQ=255)
                                                                                    cAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTg                                                                                                                                                                <  1:173953/251‑1 (MQ=255)
                                                                                     aaaCTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACTCCGTATGAACGTATCCCAGTACATCTGGTGGTGGACAACGCATGAGCCAGTGGAACGACATcagc                                                                                                                                                               >  2:330162/1‑248 (MQ=255)
                                                                                                gAACGCACCTCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTGCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCCGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGCGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGAt                                                                                                                                                      <  2:3597/249‑1 (MQ=255)
                                                                                                gAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATg                                                                                                                                                     <  2:351309/250‑1 (MQ=255)
                                                                                                     cacaTCTGGAATAAGAGATGGTGCGCCGGCGTGAAGCGGTACTGTGTGAGTGGACTGAATCGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGGACGTATCCCAGTAACTCTGGTGGTGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAAGCGATGACATcc                                                                                                                                               <  2:86449/251‑1 (MQ=255)
                                                                                                           tGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATcctgcctg                                                                                                                                         >  1:387209/1‑251 (MQ=255)
                                                                                                            ggAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGa                                                                                                                                        <  1:216591/251‑1 (MQ=255)
                                                                                                              aagnagAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCAAGGGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGaaa                                                                                                                                      <  2:369661/251‑1 (MQ=255)
                                                                                                                gaagaGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAAcc                                                                                                                                    <  1:148943/251‑1 (MQ=255)
                                                                                                                           gcgcCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTg                                                                                                                                                                            >  2:176271/1‑200 (MQ=255)
                                                                                                                           gcgcCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTg                                                                                                                                                                            <  1:176271/200‑1 (MQ=255)
                                                                                                                                      aaGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTg                                                                                                              <  1:408373/251‑1 (MQ=255)
                                                                                                                                            tACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGc                                                                                                        >  1:5796/1‑251 (MQ=255)
                                                                                                                                               tgtgtgAGTGGACTGAAACGGTCAATACTCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGCACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATAGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGAc                                                                                                            >  1:189460/1‑244 (MQ=255)
                                                                                                                                               tgtgtgAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCAGCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTTAAACCGGCGTCTGCGCGCTGTTAGGTGAc                                                                                                            <  2:189460/244‑1 (MQ=255)
                                                                                                                                                gtgtgAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAggg                                                                                                    >  2:8501/1‑250 (MQ=255)
                                                                                                                                                      gTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGTCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGAt                                                                                              <  1:367117/251‑1 (MQ=255)
                                                                                                                                                         gACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGAtttt                                                                                           <  2:218048/251‑1 (MQ=255)
                                                                                                                                                             gAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGc                                                                                       >  2:287471/1‑251 (MQ=255)
                                                                                                                                                                aCGGTCAATACGCCGTATGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCTCAACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCATCCGGCACCGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACACCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGAGACGAGCAGGTCGCGATTTTCCGCCCg                                                                                    >  1:195145/1‑251 (MQ=255)
                                                                                                                                                                     caacaCGCAGTCTGCGCAGACTCGCTTCAAACACATCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCAGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAAAGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGTGCTGTTAGGTGACGAGCAGGTCGCGATTCTCCGCCCGTATc                                                                                <  2:47589/246‑1 (MQ=255)
                                                                                                                                                                     cAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATc                                                                                >  1:44034/1‑250 (MQ=255)
                                                                                                                                                                       aTACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATcaca                                                                             <  2:156069/251‑1 (MQ=255)
                                                                                                                                                                            ccGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGAt                                                                        >  1:179738/1‑251 (MQ=255)
                                                                                                                                                                              gTCTGCGCAGCCTGGCTTCAAGCACTTCATCAACAGCGACAAGCGTGACCCGACCGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTTTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGACTGACACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCACGTATCACAGCGATc                                                                       >  2:185976/1‑250 (MQ=255)
                                                                                                                                                                              gTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATc                                                                       >  1:15178/1‑250 (MQ=255)
                                                                                                                                                                                  gcgcAGACTCGCTTCAAACACTTAAGCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAg                                                                     <  2:282033/248‑1 (MQ=255)
                                                                                                                                                                                        aCTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCg                                                            >  1:117801/1‑251 (MQ=255)
                                                                                                                                                                                           cGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATc                                                         <  2:5796/251‑1 (MQ=255)
                                                                                                                                                                                            gCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATc                                                         >  1:161201/1‑250 (MQ=255)
                                                                                                                                                                                                  aaCACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACa                                                  <  2:192644/251‑1 (MQ=255)
                                                                                                                                                                                                        tcataaACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCCGAACGCGAACCGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTCACTCTTGTGGAGGACAACGCATGCGCCAGTGGAAAGACATCTGCAAAATCGATGCCATCCTGCCTGCAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGCGTTTGCGATCAGCAACATCGAcc                                            >  2:38199/1‑251 (MQ=255)
                                                                                                                                                                                                         catcaACAGCTACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGAccc                                           <  2:117801/251‑1 (MQ=255)
                                                                                                                                                                                                            caACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGtt                                        >  1:209501/1‑251 (MQ=255)
                                                                                                                                                                                                             aaCAGCGACAAGCGGGACCCGAACGTGCAGATGGTGCTCGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGAAGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGttc                                       <  2:133769/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaGCGTGACCCGAACGTGCAGATGGTGGCATAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTcc                              <  2:331300/251‑1 (MQ=255)
                                                                                                                                                                                                                          gTGACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGGTTGCGATCAGCAACATCGACCCGTTCTTCGATTCCAGCg                          <  1:108974/251‑1 (MQ=255)
                                                                                                                                                                                                                            gACCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGt                         <  1:72779/250‑1 (MQ=255)
                                                                                                                                                                                                                             aCCCGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGc                       >  2:21960/1‑251 (MQ=255)
                                                                                                                                                                                                                                cGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGt                    >  1:406621/1‑251 (MQ=255)
                                                                                                                                                                                                                                cGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGt                    >  1:228743/1‑251 (MQ=255)
                                                                                                                                                                                                                                cGAACGTGCAGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGt                    <  1:94070/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aGATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGAc           >  1:85649/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gATGGTGCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACt          >  2:340468/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gCCAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGTAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGc    >  1:134096/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  cAGAACGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCgg  >  1:238071/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
TGCCGACAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGAACAGCACCGTCCGGCAACGCCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGG  >  NC_000913/3496239‑3496731

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