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,8930GA93.3% 121.0 / 2.2 45Q455* (CAA→TAA) abgBp‑aminobenzoyl‑glutamate hydrolase subunit B
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (23/19);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

TACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGC  >  NC_000913/1401657‑1402127
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          
tACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAAt                                                                                                                                                                                                                               >  1:40884/1‑250 (MQ=255)
  cGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCGCTTGTTAGCATTCTTGCAATAGc                                                                                                                                                                                                                            >  1:198081/1‑251 (MQ=255)
     tAATCGAATAGATAAAGAGCAGAAAAGGATGGGGAACATTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCt                                                                                                                                                                                                                         >  2:561067/1‑251 (MQ=255)
         cAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACCTACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGTTGCTCTTGTTAGCATTCTTGCAATAGCCCTGAAt                                                                                                                                                                                                                     >  1:573885/1‑251 (MQ=255)
         cAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATTCTCTTGTTAGCAGTCTTGCAATAGCCCTGAAt                                                                                                                                                                                                                     <  1:128092/251‑1 (MQ=255)
                   aaGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAag                                                                                                                                                                                                             >  1:35793/1‑249 (MQ=255)
                         aGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTCTGGATGACATACTCATGGGAATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTgcc                                                                                                                                                                                                     >  1:565848/1‑249 (MQ=255)
                           aaaaGGATGGGGCACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCACGGAAGAGATTGACtg                                                                                                                                                                                                   >  2:70056/1‑251 (MQ=255)
                             aaGGATGGGGAGCCTTGTTACCAATTCTTTCGCCCCATCCATAGAGCTTCCCGGATTTGGAGGACGACGGTATGGATGACATTCTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGTCGGTGTCACGTTTTTCGGGATAGGTCAGTGTTACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATGCGTGCAATAGCCCTGAATCAAGGAAGAGATTGACtgtg                                                                                                                                                                                                 <  2:495059/251‑1 (MQ=255)
                                   ggggAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTtgct                                                                                                                                                                                           >  1:132430/1‑251 (MQ=255)
                                                         ttCGACCCAGCCATAGAGCTGCCCGGATTGGGAGGACGACGGTATTGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAAACTAGGAAGAGATTAACTGGGGTTGCTGCCATAGTGTGCGCCGCcag                                                                                                                                                                       <  2:59310/249‑1 (MQ=255)
                                                                       agagCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGa                                                                                                                                                       <  1:208512/251‑1 (MQ=255)
                                                                        gagCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCAGAGTTTTCGCCGCCAGCAGCATTCCTTTATGAg                                                                                                                                                      <  1:199536/251‑1 (MQ=255)
                                                                             tCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCaa                                                                                                                                                   <  1:416125/249‑1 (MQ=255)
                                                                                cGGATTGGGAGGACGACGGTATGGATGCCATACTCATGGACATTACTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGCGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAAt                                                                                                                                                  >  2:28340/1‑247 (MQ=255)
                                                                                  gATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGt                                                                                                                                            <  1:70056/251‑1 (MQ=255)
                                                                                          ggaCGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCt                                                                                                                                    <  2:122360/251‑1 (MQ=255)
                                                                                                      ggATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGc                                                                                                                                 >  1:319391/1‑242 (MQ=255)
                                                                                                      ggATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGc                                                                                                                                 <  2:319391/242‑1 (MQ=255)
                                                                                                      ggATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGtt                                                                                                                        <  2:35793/251‑1 (MQ=255)
                                                                                                             aTACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTa                                                                                                                                                          <  1:110158/210‑1 (MQ=255)
                                                                                                             aTACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTa                                                                                                                                                          >  2:110158/1‑210 (MQ=255)
                                                                                                              tACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGt                                                                                                                 >  2:461052/1‑250 (MQ=255)
                                                                                                                       gCATTCCTCGGTTGTTGTGTTTGCGGTGTTGTTGTTATTTTAAACGTGACGGCGTCACGTTTTTCGGGATAGGGCAGAGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGgt                                                                                                       >  1:49329/1‑251 (MQ=255)
                                                                                                                           tCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTa                                                                                                    <  1:266851/250‑1 (MQ=255)
                                                                                                                                        tgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCaacac                                                                                       >  2:227905/1‑247 (MQ=255)
                                                                                                                                             gCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGc                                                                                 >  2:369243/1‑251 (MQ=255)
                                                                                                                                                                aaatgTGACGGAGTCACGTTTTTCAGGATAGGGCAGTGCTACGGGTGCGTGTCCGTTACTTGCTGATGCTATTGTTAGCATTCTTGAAATCGCCCTGAGTCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTAGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGAGTACCGACGTCAACACAGGGGCGGAAACACTGCGCAACAgg                                                              >  2:401245/5‑251 (MQ=255)
                                                                                                                                                                aaaGGTGACGGGGTCACGTGTTTCGGGATAGGGCAGGGAGACGGTTGCGTGTCCGGTTCTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAAGCAGGGAAGAGATTGTCTGTGGTTGCTGCCCTAGTGTTCGCCGCCAGCAGGATTCCTTTATGAGCAATAGGTGTTCGCCCCTGGCTAACCAGTTGCCACGCATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAg                                                               <  1:361090/250‑1 (MQ=255)
                                                                                                                                                                 aaGGTGACGGTATCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGCTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGc                                                             >  1:490083/1‑251 (MQ=255)
                                                                                                                                                                  aGGTGACGGGGTCCCGTTTGTCGGAATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGct                                                            <  2:327126/251‑2 (MQ=255)
                                                                                                                                                                     tGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCa                                                                                                                                                    <  1:441887/160‑1 (MQ=255)
                                                                                                                                                                     tGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCa                                                                                                                                                    >  2:441887/1‑160 (MQ=255)
                                                                                                                                                                       aCGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTc                                                       >  1:211133/1‑251 (MQ=255)
                                                                                                                                                                                      tCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACGGGGCAACAGGCAGTTTCCAACTGACGTCGcc                                         <  2:198081/250‑1 (MQ=255)
                                                                                                                                                                                              gggctgtgcgaCGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGATAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATGGACTGTGTTTGCTGCCAGAGTTTTCGCCGCCAGCAGCATGCCTTTATGAGCAATAGAGTTTCTCCCCTGGCTAACCAGTTGCCACGGATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCGGATATCAGt                                <  2:264999/241‑1 (MQ=255)
                                                                                                                                                                                                  aGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCCGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCCCGTATGTCGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCACCAGGCAGTTTCCAACTGACGTCGCCGCCATCAGTCGAt                            >  1:278222/1‑251 (MQ=255)
                                                                                                                                                                                                     gATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGc                          <  2:514350/250‑1 (MQ=255)
                                                                                                                                                                                                        aCGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGc                       >  1:564303/1‑250 (MQ=255)
                                                                                                                                                                                                             tGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGCCGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGCCATCAGTCGATGCCgcacgcac                 >  1:511531/1‑251 (MQ=255)
                                                                                                                                                                                                                 tgtCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCACCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGAGCTGAAACACTGGGCAACAGGCAGTTTCCACCTGACGTCGCCGACACCAGTCGATGCCGCAAGCACGttt             >  2:223853/1‑250 (MQ=255)
                                                                                                                                                                                                                 tgtCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCACGGAAGAGATTGACTGCGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCCTTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGACCCGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGAAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGtt              >  2:522361/1‑250 (MQ=255)
                                                                                                                                                                                                                        tACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTg       <  1:192528/250‑1 (MQ=255)
                                                                                                                                                                                                                          cTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGCCGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTggcg    >  1:87832/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCTGATGCTCTTTGTAGCATGCTTGCAATAGCCCTGAAACAAGGAAGAGATTGACTGGGGTAGCTGCCATAGTTTTCGCCGCCAGCAGCAATCCTTTATGATCAATAGATGTTCGCCCCTGGCTAACCAGTGGCCACGTATGTAGCGGTGTACCGCCGGCAAAACTGGGGCTGAAACACTGGGCAACAGGCAGTTTCCACCTGACGGCGCCGACATCAGTCGATGCCGCAAGCACGTTATCggcggcggc  <  2:513513/250‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          
TACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGC  >  NC_000913/1401657‑1402127

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