Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,401,893 G→A Q455* (CAA→TAA)  abgB ← p‑aminobenzoyl‑glutamate hydrolase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,401,8930GA98.2% 187.4 / ‑6.7 54Q455* (CAA→TAA) abgBp‑aminobenzoyl‑glutamate hydrolase subunit B
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/30);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGAGTACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAG  >  NC_000913/1401650‑1402135
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ccATGAGTACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCt                                                                                                                                                                                                                                             <  1:317468/251‑1 (MQ=255)
              aTCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGa                                                                                                                                                                                                                               <  2:455331/251‑1 (MQ=255)
                  aaTAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCa                                                                                                                                                                                                                           <  2:265114/251‑1 (MQ=255)
                   aTAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCaa                                                                                                                                                                                                                          <  1:631413/251‑1 (MQ=255)
                   aTAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCa                                                                                                                                                                                                                           >  1:249088/1‑250 (MQ=255)
                       aTAAAGAGCAGAAAAGGATGGGGAACCTTGTTACAAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGCCGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGCTACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGa                                                                                                                                                                                                                       >  2:293031/1‑250 (MQ=255)
                          aaGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAagag                                                                                                                                                                                                                   <  1:348267/251‑1 (MQ=255)
                           agagCAGAAACGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTACTCGGTTGTTGTGTTTGCGTTGTTGTTGTTTTTTAAAGGTGACGGGGTCAAGTTTTTCGGGAGAGGGGCGTGGTACGGTTGCGTGTCCGTTACTTTCTGATGCTCTTGTTAGGCTTTTTGCAATAGCCCAGAATCAAGGAagagtt                                                                                                                                                                                                                 >  2:446962/1‑249 (MQ=255)
                            gagCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCCGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCACGGAag                                                                                                                                                                                                                     >  1:144329/1‑247 (MQ=255)
                                     aGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTgg                                                                                                                                                                                                        <  2:320885/251‑1 (MQ=255)
                                        aTGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTtg                                                                                                                                                                                                     <  1:634511/251‑1 (MQ=255)
                                              aaCCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCa                                                                                                                                                                                               <  2:206172/251‑1 (MQ=255)
                                                                 tCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCcagca                                                                                                                                                                             <  1:563050/250‑1 (MQ=255)
                                                                 tCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATGCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGGCCGTTACGTGCTGATGCTCTTGGTAGCATGCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGGGGTGGCTGCCATAGTTTTCGCCGCcagca                                                                                                                                                                             <  2:699542/250‑1 (MQ=255)
                                                                  cGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTGGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCagcagc                                                                                                                                                                           <  1:483775/251‑1 (MQ=255)
                                                                  cGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTGGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGGGGTTGCTGCCATAGTTTTCGCCGCCagcagc                                                                                                                                                                           <  2:643277/251‑1 (MQ=255)
                                                                   gACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCagcagc                                                                                                                                                                           >  1:229812/1‑250 (MQ=255)
                                                                      ccAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTc                                                                                                                                                                       >  1:18822/1‑251 (MQ=255)
                                                                         gCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCtt                                                                                                                                                                    >  1:112594/1‑251 (MQ=255)
                                                                              agagCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGa                                                                                                                                                               <  2:18822/251‑1 (MQ=255)
                                                                                          aTTGGGAGGACGTCGGGATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTTTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGtt                                                                                                                                                   <  2:644607/251‑1 (MQ=255)
                                                                                          aTTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGtt                                                                                                                                                   <  2:123088/251‑1 (MQ=255)
                                                                                                     gacgGTCTGGATGACATACTCATGGGCATTCCTCGGTTGTTGCGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCCCGTTTTTCGGGCTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTATTGCAATAGCCCTGAATCAAGGACGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCt                                                                                                                                        >  2:329826/1‑251 (MQ=255)
                                                                                                       cgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGGGGCCGTTACGTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTGATGAGCAATAGATGTACGCCCCTGGCTaa                                                                                                                                      <  2:717011/251‑1 (MQ=255)
                                                                                                              gATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTg                                                                                                                               >  1:144636/1‑251 (MQ=255)
                                                                                                               aTGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGc                                                                                                                              >  2:371110/1‑251 (MQ=255)
                                                                                                                  aCATACTCATGGGCATTCCTCGGTTGTTTTGTTGGCGTTGTTGTTGTTATTTGAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTGGCGGATGCTCTTGTTAGCATTCTGGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATGGTTTTCGCCGCCAGCAGCAGTCCTTTATGAGCAAAAGATGTTCGCCGCCGACTAACCAGTGGCCAc                                                                                                                           <  2:662565/251‑1 (MQ=255)
                                                                                                                  aCATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCAc                                                                                                                           <  2:483045/251‑1 (MQ=255)
                                                                                                                    aTACTCATGGGAATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAACGGTGCCGGTGTCACGTTTTTCGGGATCGGGCCGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGCTGTTCGCCCCTGGCTAACCAGTTGCAACGt                                                                                                                         >  2:497088/1‑251 (MQ=255)
                                                                                                                                        gttgttGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCCCGTTTTTCGGTATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACgg                                                                                                     >  1:325040/1‑251 (MQ=255)
                                                                                                                                              gtgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAAc                                                                                               >  2:320854/1‑251 (MQ=255)
                                                                                                                                               tgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAAc                                                                                               >  2:238095/1‑250 (MQ=255)
                                                                                                                                                gtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTATCCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACa                                                                                              <  2:249088/250‑1 (MQ=255)
                                                                                                                                                  ttGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAggg                                                                                           >  2:484565/1‑251 (MQ=255)
                                                                                                                                                            gttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCg                                                                                                                                                 >  1:331897/1‑187 (MQ=255)
                                                                                                                                                            gttgttATTTGAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCg                                                                                                                                                 <  2:331897/187‑1 (MQ=255)
                                                                                                                                                              tgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTg                                                                               >  1:166545/1‑251 (MQ=255)
                                                                                                                                                                               cGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTcc                                                              >  2:176792/1‑251 (MQ=255)
                                                                                                                                                                               cGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACGGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCGTTATGAGCAATAGATGTTCGCCCCGGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTcc                                                              <  2:480735/251‑1 (MQ=255)
                                                                                                                                                                                   gtCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAAc                                                           <  2:397005/250‑1 (MQ=255)
                                                                                                                                                                                     cACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTg                                                         >  1:562376/1‑250 (MQ=255)
                                                                                                                                                                                          ttttCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTc                                                    <  1:503624/250‑1 (MQ=255)
                                                                                                                                                                                           tttCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAACAGATTGACGGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACGGACGTCg                                                   <  2:261487/250‑1 (MQ=255)
                                                                                                                                                                                            ttCGGGATAGGGCAGGGATACGGTTGCGTGTCCGTTACTTGCTGGTGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCAGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACTTATGTAGCGGTGTACCGCCGGCAAACCAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTctcc                                                 >  2:445579/1‑248 (MQ=255)
                                                                                                                                                                                              cGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGa                                               <  1:320854/251‑1 (MQ=255)
                                                                                                                                                                                              cGGGATAGGGCAGTGATACGGTTGCGTGGCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGGGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACATGCAGTTTCCAACTGACGTCGCCg                                                <  2:539218/250‑1 (MQ=255)
                                                                                                                                                                                                  aTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATc                                           <  2:622973/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACg                        >  2:583590/1‑251 (MQ=255)
                                                                                                                                                                                                                          tCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGACGACATCAGTCGATGCCGCAAGCACGTTATc                   >  1:267511/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATc                   <  1:484565/250‑1 (MQ=255)
                                                                                                                                                                                                                            cGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTGGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCgg                 <  1:121213/251‑1 (MQ=255)
                                                                                                                                                                                                                             gTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGt                >  2:524623/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGACGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCATGGCTAACCAGTTGCCACGTATGTAGCGGTGTAACGACGGCAAAACAGGGGATGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGCTGCCGCAAGCACGTTATCGGTGGCGGCatat      >  1:50009/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ctTGTTAGCTTTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACAGCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAg  <  2:196665/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
CCATGAGTACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAG  >  NC_000913/1401650‑1402135

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