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

GCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCTCGCCTCTGACCTCG  >  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)
                                                     gggCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGATGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑t                                                                                                                                                                                              <  1:440796/251‑1 (MQ=255)
                                                       gcgcGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGAAGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGGCGGGCTGGAACTGGCggggggcgggggggggagt                                                                                                                                                                                             >  1:179420/1‑235 (MQ=255)
                                                                  cAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAGCCATGAACCAGGTTCGCCGGCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACGCGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTc                                                                                                                                                                               <  1:252353/251‑1 (MQ=255)
                                                                  cAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑CGTGG‑TTTc                                                                                                                                                                               <  1:84974/251‑1 (MQ=255)
                                                                  cAGCTCAGTAAAGAGGACTGGCAGCAGACTTGTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTc                                                                                                                                                                               <  2:59242/251‑1 (MQ=255)
                                                                              gagGACTGGCAGCAGACTTTTGCGGTTAACGACGGCGGTGCGTTTAACCTGTTCCAGCCAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCCCTGTGGCGTCCGACGCCGCGCACACGCCCCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGGGTTGAGCGTCGGGCTGGAACTGGCGGGGAACGGCGTGCGC‑tcgtatgtgggttttcctcgcgcccc                                                                                                                                                                     >  1:536188/1‑226 (MQ=255)
                                                                                ggACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCcaccg                                                                                                                                                                 <  1:349217/251‑1 (MQ=255)
                                                                                  aCTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCcaccgac                                                                                                                                                               <  1:445188/251‑1 (MQ=255)
                                                                                         gcagACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGata                                                                                                                                                        <  1:367841/251‑1 (MQ=255)
                                                                                           agACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATg                                                                                                                                                      <  2:237556/251‑1 (MQ=255)
                                                                                                         aaCGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCATTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCtgtg                                                                                                                                        <  1:514409/251‑1 (MQ=255)
                                                                                                                             ccTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCgaag                                                                                                                    <  1:288373/251‑1 (MQ=255)
                                                                                                                                gTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCgaagaa                                                                                                                  <  2:116214/250‑1 (MQ=255)
                                                                                                                                                aaCCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑tt                                                                                                <  1:497967/251‑1 (MQ=255)
                                                                                                                                                     gTTTCGCCGTCAGCGGGGCGGGTCGATTGTCACGGTTTCGTCCGACGCCTCGCACACGCCGCGTATTGGCATGAGTGCGTATGGCGCGTCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGGTATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCg                                                                                           <  2:304135/251‑1 (MQ=255)
                                                                                                                                                       ttCGCCGGCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTGATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGa                                                                                          <  2:343077/250‑1 (MQ=255)
                                                                                                                                                       gtCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGa                                                                                          <  2:16946/249‑1 (MQ=255)
                                                                                                                                                        tCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAg                                                                                         <  1:230116/250‑1 (MQ=255)
                                                                                                                                                            cGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑tt                                                                                   <  1:382127/251‑1 (MQ=255)
                                                                                                                                                               cAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTaa                                                                                <  1:179115/251‑1 (MQ=255)
                                                                                                                                                                  cggggcggggCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACt                                                                             <  1:273339/251‑1 (MQ=255)
                                                                                                                                                                                tGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGGGAGCCGGTAAATG‑GG‑TTTCCCTTGGC‑CCACCCAACCCGGAATGCTAAGCACGCCGGGGGGGGGCGGTGAACCCGGAGAACAGCGTATTCGCGGC‑TTTGGCGGGCCG‑GTTAAAATCCGGATTCCGCT‑gg                                                              >  1:337238/1‑251 (MQ=255)
                                                                                                                                                                                 gTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑ggg                                                             <  2:109127/251‑1 (MQ=255)
                                                                                                                                                                                          gCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGc                                                    <  1:189595/251‑1 (MQ=255)
                                                                                                                                                                                                 aCGCCGCGCACACGCCGCTTATTGGCATGAGTGCGTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCca                                             <  2:545324/251‑1 (MQ=255)
                                                                                                                                                                                                    ccgcGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCCACaa                                          <  2:113234/251‑1 (MQ=255)
                                                                                                                                                                                                      gcgcACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGGAGCGGCGGGCGG‑TGGTA‑TTTGG‑GTTCGCCTGGCCCCACCGACCACGATATGCCACCCACGCTGTGGGGGAGGGGTGACGCCGAAGAACAGCGGATTTGCGGCGTTTGG‑GGGGAG‑TTTTAACTTCGCCTTCCGCT‑GGGGgaaaaagcccctcccca                                           >  1:39679/1‑235 (MQ=255)
                                                                                                                                                                                                         cacaCGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCCACAAGAGAt                                     <  2:284279/251‑1 (MQ=255)
                                                                                                                                                                                                          acacGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAAGCGCCCGTCCACAAGAGATc                                    <  1:212087/251‑1 (MQ=255)
                                                                                                                                                                                                                cgcgTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGGGGGGGG‑TGGAA‑AGTGG‑GTTTGCCTGGGTCCCCCGGCACCGATATGGAACGCACGCTGTGGGGGGGGGGTGAAGCCGAAGAACAACGTATTCGCGGC‑TTTGGGGAGCAG‑GTTTAACTCGGCATTCCGCTTGGGGAAAA‑CGCCCGTCCACAAgagaccggccac                              >  1:235656/1‑243 (MQ=255)
                                                                                                                                                                                                                   gTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTTGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAgc                                                                                                                                                                                                         <  2:299176/83‑1 (MQ=255)
                                                                                                                                                                                                                   gTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCggggagg                                                                                                                                                                                                         >  1:299176/1‑79 (MQ=255)
                                                                                                                                                                                                                        ggCATGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGAttt                       <  2:322246/250‑1 (MQ=255)
                                                                                                                                                                                                                           aTGAGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGGGGGC‑TGTTA‑TGTGGGTTTC‑CCCGGGTCCCCCGACACCCATATGCAAAGCCCGCTGTGGGGGGGCGGTGGCGCCGGAGGACAGCGGATTCGCGGC‑TTTGGGGGGGAGGTTTAAA‑TCGGGCTTCCGCT‑GGGGGAAAACGCCCGGCcaccaaaggtcgcccaaccgcatttttt                   >  1:61352/1‑226 (MQ=255)
                                                                                                                                                                                                                              aGTGCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCt                <  1:339491/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCTCg              <  2:160528/250‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCTCGc             <  2:305185/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTTATGGCGCATCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGGATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCCACAAGAGATCGCCAACCCGATTTTGTGCCTCg              <  2:368292/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aTCGAAAGCGACGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCTCGCCTCTGACCTCg  <  1:125441/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGC‑TGTAA‑TGTGG‑TTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTGTGGGTGAGCGATGACGCCGAAGAACAGCGTATTCGCGGC‑TTTGGCGAGCAG‑TTTAAACTCGGCATTCCGCT‑GGGGAAAATCGCCCGTCCACAAGAGATCGCCAACACGATTTTGTTCCTCGCCTCTGACCTCG  >  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: ‑