Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 628,998 G→A A150T (GCG→ACG)  entA → 2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913628,9980GA95.8% 153.3 / ‑6.2 49A150T (GCG→ACG) entA2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (9/37);  total (9/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCTCGCCTCTGACCTCG  >  NC_000913/628752‑629238
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
gCTGAAACGGAGCGACTGGACGCGCTGGGTAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGTGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGacgcg                                                                                                                                                                                                                                                    <  2:411278/251‑2 (MQ=255)
            cGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGc                                                                                                                                                                                                                                        <  1:103761/251‑1 (MQ=255)
                 ggACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGa                                                                                                                                                                                                                                   <  1:394345/251‑1 (MQ=255)
                    cgggcTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCGCAGTAAAGAGGACTGGCAGCAGATTTTTGCGGTTAACTGCGGCGGTGCGTGTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGTGATTGTCACTGTGGCGTACGACGCCGCGCACACGCGGCGTATTGGCATGAGTGCGTATGGCGCATCGAAAGCGGCGTTGAAAAGCATGGCGTTGAGCg                                                                                                                                                                                                                                <  2:549419/248‑1 (MQ=255)
                           gTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCt                                                                                                                                                                                                                         <  2:180878/251‑1 (MQ=255)
                              aaTGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGa                                                                                                                                                                                                                      <  1:90974/251‑1 (MQ=255)
                                       ggAAGTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCggg                                                                                                                                                                                                             <  2:309803/251‑1 (MQ=255)
                                        gAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGt                                                                                                                                                                                                            <  1:139656/251‑1 (MQ=255)
                                        gAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGGTGAGCGTCGGGCTGGAACTGGCgggg                                                                                                                                                                                                            >  1:250571/1‑250 (MQ=255)
                                         aaTTTTACGCATGGGCGCGACCGAGCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTACGCCGTCAGTGGGGCGGGGCGATTGTCGCTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGTCGCTGAGAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTa                                                                                                                                                                                                           <  1:89855/251‑1 (MQ=255)
                                             ttACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTTTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCCAGAGTGCTTATTGCGCATCGAAAGCGACGCTGAAAAGCCCGGCGTTGAGCGTCGGGCTGGAACTGGCggggggggg                                                                                                                                                                                                       >  1:198589/1‑245 (MQ=255)
                                                     gggCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGATGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCt                                                                                                                                                                                              <  1:440796/251‑1 (MQ=255)
                                                       gcgcGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGAAGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGGCGGGCTGGAACTGGCggggggcgggggggggagt                                                                                                                                                                                             >  1:179420/1‑235 (MQ=255)
                                                                  cAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAGCCATGAACCAGGTTCGCCGGCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACGCGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTc                                                                                                                                                                               <  1:252353/251‑1 (MQ=255)
                                                                  cAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAACGTGGTTTc                                                                                                                                                                               <  1:84974/251‑1 (MQ=255)
                                                                  cAGCTCAGTAAAGAGGACTGGCAGCAGACTTGTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTc                                                                                                                                                                               <  2:59242/251‑1 (MQ=255)
                                                                              gagGACTGGCAGCAGACTTTTGCGGTTAACGACGGCGGTGCGTTTAACCTGTTCCAGCCAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCCCTGTGGCGTCCGACGCCGCGCACACGCCCCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGGGTTGAGCGTCGGGCTGGAACTGGCGGGGAACGGCGTGCGCtcgtatgtgggttttcctcgcgcccc                                                                                                                                                                     >  1:536188/1‑226 (MQ=255)
                                                                                ggACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCcaccg                                                                                                                                                                 <  1:349217/251‑1 (MQ=255)
                                                                                  aCTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCcaccgac                                                                                                                                                               <  1:445188/251‑1 (MQ=255)
                                                                                         gcagACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGata                                                                                                                                                        <  1:367841/251‑1 (MQ=255)
                                                                                           agACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATg                                                                                                                                                      <  2:237556/251‑1 (MQ=255)
                                                                                                         aaCGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCATTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCtgtg                                                                                                                                        <  1:514409/251‑1 (MQ=255)
                                                                                                                             ccTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCgaag                                                                                                                    <  1:288373/251‑1 (MQ=255)
                                                                                                                                gTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCgaagaa                                                                                                                  <  2:116214/250‑1 (MQ=255)
                                                                                                                                                aaCCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCtt                                                                                                <  1:497967/251‑1 (MQ=255)
                                                                                                                                                     gTTTCGCCGTCAGCGGGGCGGGTCGATTGTCACGGTTTCGTCCGACGCCTCGCACACGCCGCGTATTGGCATGAGTGCGTATGGCGCGTCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGGTATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCg                                                                                           <  2:304135/251‑1 (MQ=255)
                                                                                                                                                       ttCGCCGGCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTGATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGa                                                                                          <  2:343077/250‑1 (MQ=255)
                                                                                                                                                       gtCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGa                                                                                          <  2:16946/249‑1 (MQ=255)
                                                                                                                                                        tCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAg                                                                                         <  1:230116/250‑1 (MQ=255)
                                                                                                                                                            cGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGtt                                                                                   <  1:382127/251‑1 (MQ=255)
                                                                                                                                                               cAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTaa                                                                                <  1:179115/251‑1 (MQ=255)
                                                                                                                                                                  cggggcggggCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACt                                                                             <  1:273339/251‑1 (MQ=255)
                                                                                                                                                                                tGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGGGAGCCGGTAAATGGGTTTCCCTTGGCCCACCCAACCCGGAATGCTAAGCACGCCGGGGGGGGGCGGTGAACCCGGAGAACAGCGTATTCGCGGCTTTGGCGGGCCGGTTAAAATCCGGATTCCGCTgg                                                              >  1:337238/1‑251 (MQ=255)
                                                                                                                                                                                 gTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTggg                                                             <  2:109127/251‑1 (MQ=255)
                                                                                                                                                                                          gCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGc                                                    <  1:189595/251‑1 (MQ=255)
                                                                                                                                                                                                 aCGCCGCGCACACGCCGCTTATTGGCATGAGTGCGTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCca                                             <  2:545324/251‑1 (MQ=255)
                                                                                                                                                                                                    ccgcGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCCACaa                                          <  2:113234/251‑1 (MQ=255)
                                                                                                                                                                                                      gcgcACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGGAGCGGCGGGCGGTGGTATTTGGGTTCGCCTGGCCCCACCGACCACGATATGCCACCCACGCTGTGGGGGAGGGGTGACGCCGAAGAACAGCGGATTTGCGGCGTTTGGGGGGAGTTTTAACTTCGCCTTCCGCTGGGGgaaaaagcccctcccca                                           >  1:39679/1‑235 (MQ=255)
                                                                                                                                                                                                         cacaCGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCCACAAGAGAt                                     <  2:284279/251‑1 (MQ=255)
                                                                                                                                                                                                          acacGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAAGCGCCCGTCCACAAGAGATc                                    <  1:212087/251‑1 (MQ=255)
                                                                                                                                                                                                                cgcgTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGGGGGGGGTGGAAAGTGGGTTTGCCTGGGTCCCCCGGCACCGATATGGAACGCACGCTGTGGGGGGGGGGTGAAGCCGAAGAACAACGTATTCGCGGCTTTGGGGAGCAGGTTTAACTCGGCATTCCGCTTGGGGAAAACGCCCGTCCACAAgagaccggccac                              >  1:235656/1‑243 (MQ=255)
                                                                                                                                                                                                                   gTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTTGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAgc                                                                                                                                                                                                         <  2:299176/83‑1 (MQ=255)
                                                                                                                                                                                                                   gTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCggggagg                                                                                                                                                                                                         >  1:299176/1‑79 (MQ=255)
                                                                                                                                                                                                                        ggCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGAttt                       <  2:322246/250‑1 (MQ=255)
                                                                                                                                                                                                                           aTGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGGGGGCTGTTATGTGGGTTTCCCCGGGTCCCCCGACACCCATATGCAAAGCCCGCTGTGGGGGGGCGGTGGCGCCGGAGGACAGCGGATTCGCGGCTTTGGGGGGGAGGTTTAAATCGGGCTTCCGCTGGGGGAAAACGCCCGGCcaccaaaggtcgcccaaccgcatttttt                   >  1:61352/1‑226 (MQ=255)
                                                                                                                                                                                                                              aGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCt                <  1:339491/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCTCg              <  2:160528/250‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCTCGc             <  2:305185/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGGATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCCACAAGAGATCGCCAACCCGATTTTGTGCCTCg              <  2:368292/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aTCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCTCGCCTCTGACCTCg  <  1:125441/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGCTTTGGCGAGCAGTTTAAACTCGGCATTCCGCTGGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCTCGCCTCTGACCTCG  >  NC_000913/628752‑629238

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