Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 169,074 G→T E531* (GAA→TAA)  fhuA → ferrichrome outer membrane transporter/phage receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913169,0740GT98.4% 202.9 / NA 63E531* (GAA→TAA) fhuAferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (37/25);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGC  >  NC_000913/168831‑169309
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           
ccagCAACCTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGt                                                                                                                                                                                                                                       <  1:434783/247‑1 (MQ=255)
 cGGCATACTACGGCCCTTAACGCATTCTGAATAAACAGAAACAAAAGGGCGTTTATGTTAAGGATCAGGCGCAGTGGGATAAAGAAATGGTCACCCTAGGCGGACGTTATGACTGGGCCGATCAAGAATCTCTTAACCGCGTTGCCGGGACGCCCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAAAGGTGTCACACCTTACTTCAGCTATAGATAATCGttt                                                                                                                                                                                                                                     >  2:376322/1‑251 (MQ=255)
       aCTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCCCCGCGGCCGCTATGACTGGGCAGATCAAGACTCTCTTACCCGCGTTGCCGGGCCGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGCGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACct                                                                                                                                                                                                                               >  1:188893/1‑251 (MQ=255)
        cTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACctt                                                                                                                                                                                                                              >  1:309535/1‑251 (MQ=255)
         cccGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTTTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACcttc                                                                                                                                                                                                                             <  2:284352/250‑1 (MQ=255)
            ggCCCTTACCGCATTCTGAATAAACAGAAACACACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTAGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGATGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACAACTTACTTCAGCTATAGCTAATCGTTTGAACCttcttc                                                                                                                                                                                                                          >  2:267495/1‑251 (MQ=255)
            ggCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTGGAACCttcttc                                                                                                                                                                                                                          <  2:82358/251‑1 (MQ=255)
                            tGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGAt                                                                                                                                                                                                          >  2:12322/1‑251 (MQ=255)
                            tGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGAt                                                                                                                                                                                                          >  1:393011/1‑251 (MQ=255)
                                  aaCAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAat                                                                                                                                                                                                    <  1:401052/251‑1 (MQ=255)
                                  aaCAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGAGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGGGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGAAat                                                                                                                                                                                                    >  2:122290/1‑251 (MQ=255)
                                   aCAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAat                                                                                                                                                                                                    >  1:321205/1‑250 (MQ=255)
                                     aGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATAtt                                                                                                                                                                                                 >  2:184626/1‑251 (MQ=255)
                                         acaaacGGGTGTTTATGTTCCGGATCAGGCTCCGTGGGGTAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGCTAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGCGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCg                                                                                                                                                                                             >  2:98426/1‑251 (MQ=255)
                                          caaacGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGc                                                                                                                                                                                            >  1:329317/1‑251 (MQ=255)
                                            aacGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCAc                                                                                                                                                                                          >  1:105650/1‑251 (MQ=255)
                                                 gCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCt                                                                                                                                                                                     <  1:252103/251‑1 (MQ=255)
                                                                aTCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAgaat                                                                                                                                                                      >  2:427658/1‑248 (MQ=255)
                                                                    ggCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTGCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAg                                                                                                                                                                  <  2:191863/251‑1 (MQ=255)
                                                                                    gTGCTGGTCACCCTAGGCGGTCGGTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGt                                                                                                                                                  <  2:104109/251‑1 (MQ=255)
                                                                                                          gatATGACTGGGCAGATCAAGCATACCTCAACCCCGTTGCCGGGACGACAGATAACCGTGATGACAAACTGTTTACCTGGCGTGGGGGTGTTACCTACCTTTTAGATAAAGGTGTAACACCTTACTTCAGCTCTAGCTAATCGTTTGCACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCCCCGTCTCAAGGTACGCCGTATGACGGCGGCGTGAAATATGCACCGGAAGATCGTCCGATTGt                                                                                                                             >  2:558813/3‑250 (MQ=255)
                                                                                                          gTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGt                                                                                                                             >  1:34069/1‑250 (MQ=255)
                                                                                                             aTGACTGGGCAGATCAAGAATCTCTTATCCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGGCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGGACCGGAAGATCGTCCGATTGTAGtt                                                                                                                         <  2:250592/251‑1 (MQ=255)
                                                                                                             aTGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGtt                                                                                                                         <  1:327043/251‑1 (MQ=255)
                                                                                                                   gggCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGGAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTg                                                                                                                     <  1:458984/249‑1 (MQ=255)
                                                                                                                    ggCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTg                                                                                                                  >  1:131047/1‑251 (MQ=255)
                                                                                                                     gCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGc                                                                                                                 <  1:91794/251‑1 (MQ=255)
                                                                                                                      cAGATCAAGAATCTCTTAACCGCGTTGCCGGGGCGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGcc                                                                                                                >  1:450655/1‑251 (MQ=255)
                                                                                                                      cAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGcc                                                                                                                >  1:5815/1‑251 (MQ=255)
                                                                                                                      cAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGcc                                                                                                                <  2:393011/251‑1 (MQ=255)
                                                                                                                       aGATCAAGCATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGAGCCg                                                                                                               >  2:408162/1‑251 (MQ=255)
                                                                                                                         aTCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGGTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCgtg                                                                                                             <  2:530817/251‑1 (MQ=255)
                                                                                                                           cAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGCGTTAACTACCTGTTCGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATCTTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTCCTGGTGCCgtgt                                                                                                            >  1:382194/1‑250 (MQ=255)
                                                                                                                                 tctcTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAAtct                                                                                                     >  2:115846/1‑251 (MQ=255)
                                                                                                                                    ctTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAAtc                                                                                                      <  1:334217/247‑1 (MQ=255)
                                                                                                                                      tAACCGCGTTGCCGGCACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCACGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACt                                                                                                 >  2:561459/1‑250 (MQ=255)
                                                                                                                                                  cGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAAcc                                                                                    >  1:195698/1‑251 (MQ=255)
                                                                                                                                                  cGGGACGACCGATAAACGTGATGACAAACAGGTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGGTTGAACCTTCGTCGCAAGTGGGGAAGGAGGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCACCCTAAATCCAACAAcc                                                                                    <  2:44934/251‑1 (MQ=255)
                                                                                                                                                   gggACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCt                                                                                   >  1:282429/1‑251 (MQ=255)
                                                                                                                                                         acCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGc                                                                             <  1:427658/251‑1 (MQ=255)
                                                                                                                                                          ccGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCg                                                                            <  2:131047/251‑1 (MQ=255)
                                                                                                                                                            gATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTTAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGa                                                                          <  2:195698/251‑1 (MQ=255)
                                                                                                                                                              tAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTAACGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGAcc                                                                        >  1:357593/1‑251 (MQ=255)
                                                                                                                                                                aaCGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCt                                                                      <  1:71454/251‑1 (MQ=255)
                                                                                                                                                                  cGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACGTCGGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGGAATATTTTCGCACCGTCTCAAGGTAAGCAGTATGAAGTCGGCGTGAAATATTTACCGGAAGATAGTCCGATTTTAGTTACGGGTGCCGTGAATAATCTCACTACAACCAACAACCTGATGGCGGACACTg                                                                     >  1:449625/1‑250 (MQ=255)
                                                                                                                                                                     gatgaCAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACGTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTGCGCAAGTTGGGAAGGATTGTAATATTTTCGCACCGGCTAAAGGTAAGCAGTAGGAAGTCGGCGTGAAATATGTACCGGAAGATCTGCCGATGGTAGTTAGTGGTGCCGTGTAGAATCTCTCTAAAACGAACAACCTCGTGGCGGATCCTGAgg                                                                  <  2:326110/250‑1 (MQ=255)
                                                                                                                                                                     gatgaCAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGGTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAgg                                                                  <  2:282429/250‑1 (MQ=255)
                                                                                                                                                                         aCAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTGCAGCTATAGCTAATCGGTTGAACCGTCTGCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTGACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTcc                                                             <  2:127190/251‑1 (MQ=255)
                                                                                                                                                                             aCAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCcttct                                                         >  1:200335/1‑251 (MQ=255)
                                                                                                                                                                               aGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCttcttc                                                       >  1:474627/1‑251 (MQ=255)
                                                                                                                                                                                            tggtggtgTTAACTAACTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAAACTACTTCGCAAGTTGGGCAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTTATGGTTCCTTCTTCTCGGTTGAAGGTg                                          >  2:504078/1‑251 (MQ=255)
                                                                                                                                                                                            gggtggtgTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCGTCTTCTCGGTTGAAGGTg                                          <  2:105650/250‑1 (MQ=255)
                                                                                                                                                                                                       aCTACCTGTTTGAGAATGGTGTAACGCGTTACGTCAGCTATAGCTAATCGGTGGAACCTTCTGCGCAAGGTGGGAAGGAGGGTAATATTTTCGCACCGGCTAAAGGTAAGCAGTATGAGGTCGGCGTGAAATATGTACCGGAAGATCGTCCGAGTGTAGGTGCTGGTGCCGTGTATAAGCTCGCTAAAACCAACAGCCTGATGGCGGGCCCTGAGGGTGCCTTCTTCGCGGTTGAAGGTGGCGAGAGCcgg                               <  2:188893/251‑2 (MQ=255)
                                                                                                                                                                                                        cTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTATTCGCAAGTTGGGAAGGATGGTAATATTTTAGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGAACCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCcgcg                              >  2:157318/1‑251 (MQ=255)
                                                                                                                                                                                                              gTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACg                          <  2:200335/249‑1 (MQ=255)
                                                                                                                                                                                                                 tGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCg                     >  1:26240/1‑251 (MQ=255)
                                                                                                                                                                                                                   aTAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTAGCACCGTCTAAAGGTAAGCAGTATGCAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGt                    >  1:29546/1‑250 (MQ=255)
                                                                                                                                                                                                                      aTGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGa                 >  1:311286/1‑250 (MQ=255)
                                                                                                                                                                                                                      aTGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAAACTGATGGCGGACCCTGCGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGa                 >  2:271672/1‑250 (MQ=255)
                                                                                                                                                                                                                        ggTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGCAGGATGGTAATATTTTCGCACCGTCTAAAGGAAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAAt              >  2:278488/1‑251 (MQ=255)
                                                                                                                                                                                                                             aaCACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGGACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGAACGTCCGATTGTAGGTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCATCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAg         <  2:468620/251‑1 (MQ=255)
                                                                                                                                                                                                                             aaCACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGCTCCGCGCACGTGGCGTAGAAATCGa           >  2:127099/1‑249 (MQ=255)
                                                                                                                                                                                                                                    tACTTCAGCAATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGGGGCGTAGAAATCGAAGCGAAAgg  >  2:263785/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           
CCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGC  >  NC_000913/168831‑169309

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