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

CCATGAGTACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑GGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAG  >  NC_000913/1401650‑1402135
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ccATGAGTACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCt                                                                                                                                                                                                                                             <  1:317468/251‑1 (MQ=255)
              aTCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGa                                                                                                                                                                                                                               <  2:455331/251‑1 (MQ=255)
                  aaTAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCa                                                                                                                                                                                                                           <  2:265114/251‑1 (MQ=255)
                   aTAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCaa                                                                                                                                                                                                                          <  1:631413/251‑1 (MQ=255)
                   aTAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCa                                                                                                                                                                                                                           >  1:249088/1‑250 (MQ=255)
                       aTAAAGAGCAGAAAAGGATGGGGAACCTTGTTACAAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGCCGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGCTACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGa                                                                                                                                                                                                                       >  2:293031/1‑250 (MQ=255)
                          aaGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAagag                                                                                                                                                                                                                   <  1:348267/251‑1 (MQ=255)
                           agagCAGAAACGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTACTCGGTTGTTGTGTTTGCGTTGTTGTTGTT‑TTTTAAAGGTGACGGGGTCAAGTTTTTCGGGAGAGGGGCGTGGTACGGTTGCGTGTCCGTTACTTTCTGATGCTCTTGTTAGGC‑TTTTTGCAATAGCCCAGAATCAAGGAagagtt                                                                                                                                                                                                                 >  2:446962/1‑249 (MQ=255)
                            gagCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCCGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCACGGAag                                                                                                                                                                                                                     >  1:144329/1‑247 (MQ=255)
                                     aGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTgg                                                                                                                                                                                                        <  2:320885/251‑1 (MQ=255)
                                        aTGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTtg                                                                                                                                                                                                     <  1:634511/251‑1 (MQ=255)
                                              aaCCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCa                                                                                                                                                                                               <  2:206172/251‑1 (MQ=255)
                                                                 tCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCcagca                                                                                                                                                                             <  1:563050/250‑1 (MQ=255)
                                                                 tCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATGCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGGCCGTTACGTGCTGATGCTCTTGGT‑AGCATGCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGGGGTGGCTGCCATAGTTTTCGCCGCcagca                                                                                                                                                                             <  2:699542/250‑1 (MQ=255)
                                                                  cGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTGGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCagcagc                                                                                                                                                                           <  1:483775/251‑1 (MQ=255)
                                                                  cGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTGGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGGGGTTGCTGCCATAGTTTTCGCCGCCagcagc                                                                                                                                                                           <  2:643277/251‑1 (MQ=255)
                                                                   gACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCagcagc                                                                                                                                                                           >  1:229812/1‑250 (MQ=255)
                                                                      ccAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTc                                                                                                                                                                       >  1:18822/1‑251 (MQ=255)
                                                                         gCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCtt                                                                                                                                                                    >  1:112594/1‑251 (MQ=255)
                                                                              agagCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGa                                                                                                                                                               <  2:18822/251‑1 (MQ=255)
                                                                                          aTTGGGAGGACGTCGGGATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTTTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGtt                                                                                                                                                   <  2:644607/251‑1 (MQ=255)
                                                                                          aTTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGtt                                                                                                                                                   <  2:123088/251‑1 (MQ=255)
                                                                                                     gacgGTCTGGATGACATACTCATGGGCATTCCTCGGTTGTTGCGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCCCGTTTTTCGGGCTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTATTGCAATAGCCCTGAATCAAGGACGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCt                                                                                                                                        >  2:329826/1‑251 (MQ=255)
                                                                                                       cgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGGGGCCGTTACGTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTGATGAGCAATAGATGTACGCCCCTGGCTaa                                                                                                                                      <  2:717011/251‑1 (MQ=255)
                                                                                                              gATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTg                                                                                                                               >  1:144636/1‑251 (MQ=255)
                                                                                                               aTGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGc                                                                                                                              >  2:371110/1‑251 (MQ=255)
                                                                                                                  aCATACTCATGGGCATTCCTCGGTTGTTTTGTTGGCGTTGTTGTTGTTATTTGAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTGGCGGATGCTCTTGTT‑AGCATTCTGGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATGGTTTTCGCCGCCAGCAGCAGTCCTTTATGAGCAAAAGATGTTCGCCGCCGACTAACCAGTGGCCAc                                                                                                                           <  2:662565/251‑1 (MQ=255)
                                                                                                                  aCATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCAc                                                                                                                           <  2:483045/251‑1 (MQ=255)
                                                                                                                    aTACTCATGGGAATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAACGGTGCCGGTGTCACGTTTTTCGGGATCGGGCCGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGCTGTTCGCCCCTGGCTAACCAGTTGCAACGt                                                                                                                         >  2:497088/1‑251 (MQ=255)
                                                                                                                                        gttgttGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCCCGTTTTTCGGTATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACgg                                                                                                     >  1:325040/1‑251 (MQ=255)
                                                                                                                                              gtgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAAc                                                                                               >  2:320854/1‑251 (MQ=255)
                                                                                                                                               tgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAAc                                                                                               >  2:238095/1‑250 (MQ=255)
                                                                                                                                                gtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTATCCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACa                                                                                              <  2:249088/250‑1 (MQ=255)
                                                                                                                                                  ttGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAggg                                                                                           >  2:484565/1‑251 (MQ=255)
                                                                                                                                                            gttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCg                                                                                                                                                 >  1:331897/1‑187 (MQ=255)
                                                                                                                                                            gttgttATTTGAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCg                                                                                                                                                 <  2:331897/187‑1 (MQ=255)
                                                                                                                                                              tgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTg                                                                               >  1:166545/1‑251 (MQ=255)
                                                                                                                                                                               cGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTcc                                                              >  2:176792/1‑251 (MQ=255)
                                                                                                                                                                               cGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACGGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCGTTATGAGCAATAGATGTTCGCCCCGGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTcc                                                              <  2:480735/251‑1 (MQ=255)
                                                                                                                                                                                   gtCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAAc                                                           <  2:397005/250‑1 (MQ=255)
                                                                                                                                                                                     cACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTg                                                         >  1:562376/1‑250 (MQ=255)
                                                                                                                                                                                          ttttCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTc                                                    <  1:503624/250‑1 (MQ=255)
                                                                                                                                                                                           tttCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAACAGATTGACGGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACGGACGTCg                                                   <  2:261487/250‑1 (MQ=255)
                                                                                                                                                                                            ttCGGGATAGGGCAGGGATACGGTTGCGTGTCCGTTACTTGCTGGTGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCAGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACTTATGTAGCGGTGTACCGCCGGCAAACCAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTctcc                                                 >  2:445579/1‑248 (MQ=255)
                                                                                                                                                                                              cGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGa                                               <  1:320854/251‑1 (MQ=255)
                                                                                                                                                                                              cGGGATAGGGCAGTGATACGGTTGCGTGGCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGGGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACATGCAGTTTCCAACTGACGTCGCCg                                                <  2:539218/250‑1 (MQ=255)
                                                                                                                                                                                                  aTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATc                                           <  2:622973/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCGTGTCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACg                        >  2:583590/1‑251 (MQ=255)
                                                                                                                                                                                                                          tCCGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGACGACATCAGTCGATGCCGCAAGCACGTTATc                   >  1:267511/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATc                   <  1:484565/250‑1 (MQ=255)
                                                                                                                                                                                                                            cGTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTGGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCgg                 <  1:121213/251‑1 (MQ=255)
                                                                                                                                                                                                                             gTTACTTGCTGATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGt                >  2:524623/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gATGCTCTTGTT‑AGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGACGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCATGGCTAACCAGTTGCCACGTATGTAGCGGTGTAACGACGGCAAAACAGGGGATGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGCTGCCGCAAGCACGTTATCGGTGGCGGCatat      >  1:50009/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ctTGTT‑AGCTTTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACAGCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAg  <  2:196665/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
CCATGAGTACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTT‑GGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAG  >  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: ‑