Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,904,420 A→T 6.6% I176F (ATC→TTC)  manZ → mannose‑specific enzyme IID component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,904,4200AT 6.6% 202.9 / 5.2 61I176F (ATC→TTC) manZmannose‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base T (2/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

GCGTAACCCTGGCGCTGGAAGAACAGCGTGCTAATGGCGCAGAGATCGACGACGGTGCTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCT  >  NC_000913/1904145‑1904696
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                     
gCGTAACCCTGGCGCTGGAAGAACAGCGTGCTAATGGCGCAGAGATCGACGACGGTGCTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGata                                                                                                                                                                                                                                                                           <  1:29516/288‑1 (MQ=255)
    aaCCCTGGCGCTGGAAGAACAGCGTGCTAATGGCGCAGAGATCGACGACGGTGCTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTATGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGCAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATggg                                                                                                                                                                                                                                                                       <  1:484059/288‑1 (MQ=255)
           gcgcTGGAAGAACAGCGTGCTAATGGCGCAGAGATCGACGACGGTGCTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCt                                                                                                                                                                                                                                                               >  1:418998/1‑289 (MQ=255)
                         gCGTGCTAATGGCGCAGAGATCGACGACGGTGCTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAACAACTg                                                                                                                                                                                                                                                 >  1:115354/1‑289 (MQ=255)
                           gTGCTAATGGCGCAGAGATCGACGACGGTGCTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACt                                                                                                                                                                                                                                                  <  1:90649/286‑1 (MQ=255)
                                        agagATCGACGACGGTGCTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCCTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTACTGCAAAACCTGACGGAAGGGGCGTca                                                                                                                                                                                                                                  >  1:484047/1‑288 (MQ=255)
                                               gacgacGGTGCTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATTTCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCg                                                                                                                                                                                                                           >  2:115646/1‑289 (MQ=255)
                                                 cgacgGTGCTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGc                                                                                                                                                                                                                         <  2:146132/289‑1 (MQ=255)
                                                         cTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTg                                                                                                                                                                                                                  >  1:235924/1‑288 (MQ=255)
                                                            tCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCAt                                                                                                                                                                                                              <  2:418998/289‑1 (MQ=255)
                                                             cAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCTTACGTCCGGTATGTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATg                                                                                                                                                                                                             <  2:465939/289‑1 (MQ=255)
                                                                   tATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGGCTATATCCTTGGCCTGGTTGTCCTGGGGGCa                                                                                                                                                                                                       >  1:357320/1‑289 (MQ=255)
                                                                    aTCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCa                                                                                                                                                                                                       <  1:475518/288‑1 (MQ=255)
                                                                            cGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTaa                                                                                                                                                                                               <  1:127499/288‑1 (MQ=255)
                                                                            cGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAAc                                                                                                                                                                                              <  2:127522/289‑1 (MQ=255)
                                                                               tttGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAg                                                                                                                                                                                           >  2:269241/1‑289 (MQ=255)
                                                                                       gtgcACTGGCTGGCGTAGGCGGCCCGATCTTCTGCGGAACCGTACGTCCGGTATTTGCAGCACTGGTTGCAGGTATCGCGATGAGCGGTAGCCTGTTATGTCCGCTGCTGTTCTGCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTAGCGATATCGTTAAAGATATGGGTGGTGGCTTCCTACAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGacac                                                                                                                                                                                    <  2:352788/285‑1 (MQ=255)
                                                                                           aCTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGc                                                                                                                                                                                                                         <  1:345965/247‑1 (MQ=255)
                                                                                           aCTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGc                                                                                                                                                                                                                         >  2:345965/1‑247 (MQ=255)
                                                                                             tggctggTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCaa                                                                                                                                                                             >  1:147347/1‑289 (MQ=255)
                                                                                             tggctggTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCa                                                                                                                                                                              >  1:514237/1‑288 (MQ=255)
                                                                                                     gtAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCg                                                                                                                                                                      <  2:115354/288‑1 (MQ=255)
                                                                                                         gCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATTTCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGt                                                                                                                                                                >  1:80918/1‑289 (MQ=255)
                                                                                                         gCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGt                                                                                                                                                                >  1:350881/1‑289 (MQ=255)
                                                                                                             cccGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCTGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTg                                                                                                                                                              >  2:47759/1‑287 (MQ=255)
                                                                                                             cccGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTg                                                                                                                                                              >  1:235949/1‑287 (MQ=255)
                                                                                                              ccGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGtc                                                                                                                                                            >  2:366601/1‑288 (MQ=255)
                                                                                                                  tcttctGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATTTCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCg                                                                                                                                                       <  2:80918/289‑1 (MQ=255)
                                                                                                                  tcttctGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCg                                                                                                                                                       <  2:506999/289‑1 (MQ=255)
                                                                                                                            aaCCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGa                                                                                                                                              >  2:102538/1‑288 (MQ=255)
                                                                                                                              ccGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACCAATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTTACCa                                                                                                                                           >  1:85979/1‑289 (MQ=255)
                                                                                                                                 tACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGCAATTGGTAACAAAGTGGCAACTTGTCACACTCCGGTGGTTTGTCTCTCGCTTTCTGAcccggcg                                                                                                                                       >  2:126488/1‑284 (MQ=255)
                                                                                                                                      ccGGTATTTGCAGCACTGGGTGCTGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCa                                                                                                                                   <  1:47759/289‑1 (MQ=255)
                                                                                                                                           aTTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGaa                                                                                                                              <  1:366601/289‑1 (MQ=255)
                                                                                                                                           aTTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGaa                                                                                                                              <  2:379727/289‑1 (MQ=255)
                                                                                                                                           aTTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGaa                                                                                                                              <  2:235949/289‑1 (MQ=255)
                                                                                                                                             ttGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAaca                                                                                                                            <  2:107965/289‑1 (MQ=255)
                                                                                                                                                    aCTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGACAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGCCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTtac                                                                                                                      >  2:102568/1‑288 (MQ=255)
                                                                                                                                                      tGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCGGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGGATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTtacta                                                                                                                    <  1:352353/288‑1 (MQ=255)
                                                                                                                                                          tGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATTTCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGt                                                                                                                <  1:115646/288‑1 (MQ=255)
                                                                                                                                                             cGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTaa                                                                                                                                                                                               <  2:176910/207‑1 (MQ=255)
                                                                                                                                                             cGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTaa                                                                                                                                                                                               >  1:176910/1‑207 (MQ=255)
                                                                                                                                                              ggTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGa                                                                                                           <  2:147347/289‑1 (MQ=255)
                                                                                                                                                                    gcgATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTAt                                                                                                       >  1:152179/1‑287 (MQ=255)
                                                                                                                                                                                    tGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGcc                                                                                     <  2:235924/289‑1 (MQ=255)
                                                                                                                                                                                       tAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAgg                                                                                  <  2:484047/289‑1 (MQ=255)
                                                                                                                                                                                                   tgTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAATGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCAc                                                                       <  2:148057/288‑1 (MQ=255)
                                                                                                                                                                                                            tCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAc                                                             <  2:85979/289‑1 (MQ=255)
                                                                                                                                                                                                              cTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCt                                                           <  2:350881/289‑1 (MQ=255)
                                                                                                                                                                                                                     aCCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                                    <  2:357320/289‑1 (MQ=255)
                                                                                                                                                                                                                       cTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTa                                                  <  2:152179/289‑1 (MQ=255)
                                                                                                                                                                                                                       cTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTa                                                  <  2:514237/289‑1 (MQ=255)
                                                                                                                                                                                                                          gTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtgt                                               >  1:235444/1‑289 (MQ=255)
                                                                                                                                                                                                                           tGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACACGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGCACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtttg                                              >  2:447254/1‑289 (MQ=255)
                                                                                                                                                                                                                                        ccGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaaa                                 <  1:414139/289‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTaa                         >  2:342677/1‑288 (MQ=255)
                                                                                                                                                                                                                                                             tGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGaa                                                                                                                              >  1:482483/1‑175 (MQ=255)
                                                                                                                                                                                                                                                             tGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGaa                                                                                                                              <  2:482483/175‑1 (MQ=255)
                                                                                                                                                                                                                                                              ggTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatca           >  2:285398/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                   cTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTATCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGtt      >  2:168647/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                        aaGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGct  >  2:108037/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                     
GCGTAACCCTGGCGCTGGAAGAACAGCGTGCTAATGGCGCAGAGATCGACGACGGTGCTATCAACGGTATCAAAGTCGGTTTGATGGGGCCACTGGCTGGTGTAGGCGACCCGATCTTCTGGGGAACCGTACGTCCGGTATTTGCAGCACTGGGTGCCGGTATCGCGATGAGCGGCAGCCTGTTAGGTCCGCTGCTGTTCTTCATCCTGTTTAACCTGGTGCGTCTGGCAACCCGTTACTACGGCGTAGCGTATGGTTACTCCAAAGGTATCGATATCGTTAAAGATATGGGTGGTGGCTTCCTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCT  >  NC_000913/1904145‑1904696

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