Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 924,859 T→C R532R (CGT→CGC clpA → ATP‑dependent Clp protease ATP‑binding subunit ClpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913924,8590TC96.5% 184.0 / ‑4.2 58R532R (CGT→CGCclpAATP‑dependent Clp protease ATP‑binding subunit ClpA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (31/25);  minor base G (0/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

AAAACCTCGGCGATCGCTTGAAAATGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGTCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCG  >  NC_000913/924612‑925104
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
aaaaCCTCGGCGATCGCTTGAAAATGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGCACGCCATACCGTCAGCCGCCtt                                                                                                                                                                                                                                                    >  1:371286/1‑251 (MQ=255)
       cGGCGATCGCTTGAAAATGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGt                                                                                                                                                                                                                                              <  1:72573/250‑1 (MQ=255)
        ggCGATCGCTTGAAAATGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTgg                                                                                                                                                                                                                                            >  2:203467/1‑250 (MQ=255)
           gATCGCTTGAAAATGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGcc                                                                                                                                                                                                                                         >  1:168676/1‑251 (MQ=255)
              cGCTTGAAAATGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTcc                                                                                                                                                                                                                                      <  2:393293/251‑1 (MQ=255)
                 ttGAAAATGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCggg                                                                                                                                                                                                                                   >  1:135247/1‑251 (MQ=255)
                    aaaaTGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGACGCCATTAAGATGGCGCGTGCAGGTTTAGGTCATGACAATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATa                                                                                                                                                                                                                                <  1:390812/251‑1 (MQ=255)
                    aaaaTGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGTCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATa                                                                                                                                                                                                                                <  1:178426/251‑1 (MQ=255)
                    aaaaTGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATa                                                                                                                                                                                                                                <  1:220000/251‑1 (MQ=255)
                            ggacgtCGGTCAGGATAAAGCCATTGAGGCGCTGACTGGAGCCACTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAGTCCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGTCGGCACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTTCGCTTTGATATGTCCGATTGTATGGAACGCCGTACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGtt                                                                                                                                                                                                                        <  1:428471/246‑1 (MQ=255)
                             gTCTTCGGTCAGGATAAAGCCATTGAGTCGCGGACTGAAGCCATTAAGATGGCGCGTTCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCTGGAAAACAGAGTTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATTCCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATTCGTTGGttt                                                                                                                                                                                                                       <  2:188056/251‑1 (MQ=255)
                             gTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGttt                                                                                                                                                                                                                       <  1:326290/251‑1 (MQ=255)
                                 tCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGa                                                                                                                                                                                                                    <  1:167762/250‑1 (MQ=255)
                                  cGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATc                                                                                                                                                                                                                  <  1:334031/251‑1 (MQ=255)
                                                 cATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGctg                                                                                                                                                                                                    >  2:386392/1‑250 (MQ=255)
                                                          gcTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCgg                                                                                                                                                                                          >  1:246812/1‑251 (MQ=255)
                                                                 gaAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGtcatca                                                                                                                                                                                    >  1:22442/1‑250 (MQ=255)
                                                                      cATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATc                                                                                                                                                                              <  2:80217/251‑1 (MQ=255)
                                                                             aTGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGc                                                                                                                                                                       <  1:458763/251‑1 (MQ=255)
                                                                             aTGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGc                                                                                                                                                                       >  2:106340/1‑251 (MQ=255)
                                                                               ggCGCGTGCATGTTTAGGTCACGAACATAATCCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGTGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCATATGCgg                                                                                                                                                                     <  1:362980/251‑1 (MQ=255)
                                                                                gcgcGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGt                                                                                                                                                                    <  1:191350/251‑1 (MQ=255)
                                                                                   cgTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGtgct                                                                                                                                                                 <  1:319769/251‑1 (MQ=255)
                                                                                         ggTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTgctgctgct                                                                                                                                                           >  2:121322/1‑251 (MQ=255)
                                                                                                    cGAACATAAGCCGGTGGGTTCGTTCCTGTTTGTCGGCGCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGTGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATTCCGTCAGCCGGCTTATTGGTGCACCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTTCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCg                                                                                                                                                <  1:413558/251‑1 (MQ=255)
                                                                                                      aaCATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCgag                                                                                                                                              <  2:292203/251‑1 (MQ=255)
                                                                                                               cGGTTGGTTCGTTCCTGTTTGCCGGCACTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCAc                                                                                                                                     <  2:133738/251‑1 (MQ=255)
                                                                                                                      ttcgttcCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGAACAGCTTTCGCAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTACGGGATACGTTGGTTTTGATCAGGGAGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACg                                                                                                                              >  1:52440/1‑251 (MQ=255)
                                                                                                                             cTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCa                                                                                                                        >  1:328597/1‑250 (MQ=255)
                                                                                                                              tGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAat                                                                                                                      >  1:14425/1‑251 (MQ=255)
                                                                                                                                   gCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCt                                                                                                                 <  2:97959/251‑1 (MQ=255)
                                                                                                                                               ggggTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCGTCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGAt                                                                                                     >  1:267906/1‑251 (MQ=255)
                                                                                                                                               ggggTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGAt                                                                                                     >  1:75428/1‑251 (MQ=255)
                                                                                                                                                             agagGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGc                                                                                       >  1:125677/1‑251 (MQ=255)
                                                                                                                                                                   gaaggTACAGCTTTCGAAAGCTTTGGGCATTGAGCTACTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGAGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGAGTTCAATATACTGTTGCAGGTGATGGATAACGGTCCGCTGGCCg                                                                                 >  2:126928/4‑251 (MQ=255)
                                                                                                                                                                   gACGGTACAGCTTTCGAATGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCAGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGAcc                                                                                  <  1:346332/250‑1 (MQ=255)
                                                                                                                                                                   gACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCg                                                                                 >  2:6249/1‑251 (MQ=255)
                                                                                                                                                                       gTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATa                                                                              >  1:144337/1‑250 (MQ=255)
                                                                                                                                                                          cAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTGGATATGTCCGAGTATATGGAAGGCGATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTGTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCTGACGTGTTCAATATTCTGGTGCAGGTGATGGATAACGGTACGCTGACCGATaacaa                                                                          <  1:126928/251‑1 (MQ=255)
                                                                                                                                                                                 cGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCATCCGGGATACGTTGGCTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATCACGGTACGCTGACCGATAACAACGGACGc                                                                   >  2:133973/1‑251 (MQ=255)
                                                                                                                                                                                   aaaGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCa                                                                  >  2:88998/1‑250 (MQ=255)
                                                                                                                                                                                          tGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGa                                                           <  1:106340/250‑1 (MQ=255)
                                                                                                                                                                                             gcatcgtgattaggcgcTATGATATGTCCGAGTTTACTGAACACCCTACCGTCAGCCAGCTGATTGGTGCGGCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCTGTCATCAATCATCTACAGGCGGTGATGCTGTTGGTCGATAGCGAGAAAGCGCACCCGGACGTGTGCAATATTCTGTTGCAGGTGATGGAGAAGGGTATGCTGGCATATACCAACGGTCGTAAAGCAGACtac                                                       <  2:25165/238‑3 (MQ=255)
                                                                                                                                                                                                        ctgcgtcttGATATTTCCGAGGATATGCAACGCCATACCGTCAGCCCCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATCCGGTCATCAAGCATCCACCTGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGCGTTCAATATTCTGTTGCCGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAGCgtggg                                            >  2:55091/8‑250 (MQ=255)
                                                                                                                                                                                                             cTTTGATATGTCCGAGTATATGGCACGCCATACCGTCAGCCCCCTTATTGGCGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAACTCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTgg                                       >  2:448950/1‑251 (MQ=255)
                                                                                                                                                                                                                 gATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGtgat                                   <  2:22442/251‑1 (MQ=255)
                                                                                                                                                                                                                    aggTCCTAGTATACTGAACTCCATACCTTCAGTCGCCCTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCTGTGCTGCTGCTGGCCGAAATCGAGAAAGCGCACCCGGACGTGTTAAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGGCGCAAAGCAGACTTCCGTAACGTGGGGCTGGTGATGac                                <  2:222861/249‑1 (MQ=255)
                                                                                                                                                                                                                       tCCGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGTGATGaccac                             >  2:139597/1‑251 (MQ=255)
                                                                                                                                                                                                                         cGAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGTGATGaccacc                            >  2:144554/1‑250 (MQ=255)
                                                                                                                                                                                                                          gAGTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGATGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCTGGTGATGGATAACGGTACGCTGAACGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGTGATGAccacca                           >  2:103239/1‑250 (MQ=255)
                                                                                                                                                                                                                            gTATATGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGTGATGACCACCAACg                        >  2:433188/1‑251 (MQ=255)
                                                                                                                                                                                                                             tataTGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCACAGCAGACTTCCGTAACGTGGTGCTGGTGATGACCACCAACGc                       >  1:96891/1‑251 (MQ=255)
                                                                                                                                                                                                                              atatGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGTGATGACCACCAACGcc                      >  1:289126/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tGGAACGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGCTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGAAGGTGATGGATAACGGAACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTACCGTGGTGCTGGTGATGACCACCAACGCCggg                   >  2:353487/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aCGCCATACCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCAGCCACATGCGGTGCTGCTGCTTGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGGTGGTGATGACCACCAACGCCGGGGTAc               <  1:448950/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aCCGTCAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGTGACGACCACCAACGCCGGGGTACGGGAAAc        >  2:97352/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cAGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAgc   <  1:386392/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  aGCCGCCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAgcg  >  2:259828/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
AAAACCTCGGCGATCGCTTGAAAATGCTGGTCTTCGGTCAGGATAAAGCCATTGAGGCGCTGACTGAAGCCATTAAGATGGCGCGTGCAGGTTTAGGTCACGAACATAAACCGGTTGGTTCGTTCCTGTTTGCCGGCCCTACCGGGGTCGGGAAAACAGAGGTGACGGTACAGCTTTCGAAAGCTTTGGGCATTGAGCTTCTGCGCTTTGATATGTCCGAGTATATGGAACGCCATACCGTCAGCCGTCTTATTGGTGCGCCTCCGGGATACGTTGGTTTTGATCAGGGCGGTTTGCTGACTGATGCGGTCATCAAGCATCCACATGCGGTGCTGCTGCTGGACGAAATCGAGAAAGCGCACCCGGACGTGTTCAATATTCTGTTGCAGGTGATGGATAACGGTACGCTGACCGATAACAACGGACGCAAAGCAGACTTCCGTAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCG  >  NC_000913/924612‑925104

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