Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 458,742 T→A 5.6% intergenic (+42/‑146) clpX → / → lon ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease/DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913458,7420TA 5.6% 223.4 / 2.1 81intergenic (+42/‑146)clpX/lonATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease/DNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (33/43);  new base A (2/3);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTAT  >  NC_000913/458594‑458889
                                                                                                                                                    |                                                                                                                                                   
tGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGCTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGAtt                                                                                                                                                     >  1:1023704/1‑149 (MQ=255)
   aaGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTa                                                                                                                                                  >  1:948219/1‑149 (MQ=255)
   aaGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTa                                                                                                                                                  >  1:352792/1‑149 (MQ=255)
    aGACGTCGAACAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAt                                                                                                                                                 >  2:24941/1‑149 (MQ=255)
    aGACGTCGAAAAAGTGGTTATCTACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAt                                                                                                                                                 <  2:544927/149‑1 (MQ=255)
    aGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAt                                                                                                                                                 <  2:540674/149‑1 (MQ=255)
        ttctaaaaaTTGGTTATCGACGAGTAGGTAATTGATTTTCAAAGCAAACCTATGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAACAAGGGGGGATTTTATCTcc                                                                                                                                             <  1:856437/145‑1 (MQ=255)
        gTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATTATTAAGCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATATcc                                                                                                                                             >  2:608375/1‑149 (MQ=255)
        gTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTcc                                                                                                                                             >  2:675785/1‑149 (MQ=255)
          cGAAAAAGTGGGTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGAGTTATGGCAAGCCGGAAGCGCAACAGGCAACTGGTGAATAATTACCCATTACCATACAATTAGTTAACCAAAAAGGGGGTATTTTATCTcccc                                                                                                                                           >  2:1077524/1‑149 (MQ=255)
               aaGGGGGTATCGACGAGTCTGTAATTGATTGTAAAATCAAACCGTTGCTTATTAATGTCAAGCCGGAAGCTCAACAGGCATCTTGTGAATACTTAACCATTCCCATACATTTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCttttt                                                                                                                                      <  1:663082/149‑3 (MQ=255)
                aGTGGTTATCGACGAGTCGGTTATTGATGGTCAAATCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttg                                                                                                                                     <  2:1001270/149‑4 (MQ=255)
                aGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttg                                                                                                                                     <  2:983839/149‑4 (MQ=255)
                   ggTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAtttt                                                                                                                                  <  1:483829/149‑1 (MQ=255)
                   ggTTATCGACGAGGCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCTGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCACTTTAAtttt                                                                                                                                  <  1:949332/149‑1 (MQ=255)
                    gTTATCGACGAGTCGGTAATTGATGGTCTAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGCGGGGATTTTATCTCCCCTTTAAttttt                                                                                                                                 >  1:991553/1‑149 (MQ=255)
                    gTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGACGCACAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTGAAttttt                                                                                                                                 >  1:492271/1‑149 (MQ=255)
                     ttATCGCCGAGTCGGTACTTGATGGTCAACGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTc                                                                                                                                <  1:1029802/149‑1 (MQ=255)
                     ttATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAATCGCAACAGGCATCTGGTGAATAAATAACCATTCCCATACAATTAGATAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTc                                                                                                                                <  1:894216/149‑1 (MQ=255)
                     ttATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTc                                                                                                                                <  1:167132/149‑1 (MQ=255)
                     ttATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCAGCTGGGGAATAATTCACCATTCCCATACAATTAGTTAACCAAAACGGGGGGATTTTATCTCCCCTTTAATTTTTc                                                                                                                                <  1:857419/149‑1 (MQ=255)
                       aTCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCct                                                                                                                              <  1:147117/149‑1 (MQ=255)
                        tCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCctc                                                                                                                             >  2:316898/1‑149 (MQ=255)
                         cgacgaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCATTTAATTTTTCctct                                                                                                                            >  2:460994/1‑149 (MQ=255)
                            cgaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTAtt                                                                                                                         <  2:153072/149‑1 (MQ=255)
                             gaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCtttta                                                                                                                        >  2:580330/1‑148 (MQ=255)
                               gtaggTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtctc                                                                                                                      <  2:51058/146‑1 (MQ=255)
                                tctgTAATTGATTGTCAAAGCAAACCATTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCg                                                                                                                     <  2:216723/146‑1 (MQ=255)
                                       ttGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGa                                                                                                              <  2:948219/149‑1 (MQ=255)
                                         taTGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCTCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAAt                                                                                                            <  1:663877/148‑1 (MQ=255)
                                          aTGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAAtg                                                                                                           >  2:963071/1‑149 (MQ=255)
                                             gTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTgg                                                                                                        <  1:189313/149‑1 (MQ=255)
                                             gTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTgg                                                                                                        <  1:337793/149‑1 (MQ=255)
                                             gTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTgg                                                                                                        <  2:1058296/149‑1 (MQ=255)
                                              tCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACAATTCCCATACAATTAGTTAACCAAAAAGGTGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCTGCGTTGAATGTggg                                                                                                       >  2:73069/1‑149 (MQ=255)
                                                    cAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAAc                                                                                                 <  1:500507/149‑1 (MQ=255)
                                                               tgtattATGGCAAGCCTGAAGCTCAACAGGCATCTGGTTTATAATTAACCATTGCCATACAATTAGTTAAGCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTAATCTCGGCGTTTAATGTGGGGGAAACATCCCCAtata                                                                                      <  1:369914/145‑1 (MQ=255)
                                                               tGATTTATGGCAAGCCGGAAGCGCAACCGGCATCTGGTGAATAATTAACCATTCCCATACAAGTAGTTAACCAAAAAGGGGGGAGTTTATGTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAAGGTGGGGGAAAAATCCCCAtata                                                                                      >  2:1145309/1‑149 (MQ=255)
                                                                   ttATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGa                                                                                  <  2:57327/149‑1 (MQ=255)
                                                                    tATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGAc                                                                                 <  2:55055/149‑1 (MQ=255)
                                                                     aTGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACg                                                                                >  2:92015/1‑149 (MQ=255)
                                                                       gtcAAGTCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCTTTCCCATTATATTAGTTAACCAAATAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTa                                                                              <  1:1124226/147‑1 (MQ=255)
                                                                       ggCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTa                                                                              >  1:662668/1‑149 (MQ=255)
                                                                        gCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATGCAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGAGTTGAATGTGGGGGAAACATCCCCATATACTGACGTAc                                                                             >  2:623651/1‑149 (MQ=255)
                                                                        gCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTAc                                                                             >  2:299627/1‑149 (MQ=255)
                                                                             ccGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGGGTAGAATGTGGGGGAAACATCCCCATATACTGACGTACATCtt                                                                        >  1:341051/1‑149 (MQ=255)
                                                                                     gcAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATg                                                                <  2:302799/149‑1 (MQ=255)
                                                                                       aaCAGGCATCTGCTGAATAATTAACCATTCCCATACAATTACTAAACCAAAAAGGGGGAATTTTATCACCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTTGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGc                                                              <  2:169368/149‑1 (MQ=255)
                                                                                         cAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTGAATTCTCTGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGt                                                            >  2:170381/1‑149 (MQ=255)
                                                                                             cATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAg                                                        <  1:506902/149‑1 (MQ=255)
                                                                                             cATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAg                                                        <  2:949836/149‑1 (MQ=255)
                                                                                             cATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTCGACTGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGGGAAg                                                        <  1:692444/149‑1 (MQ=255)
                                                                                               tatGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAATGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGca                                                      <  2:145069/147‑1 (MQ=255)
                                                                                               tCTGGTGAATAATAAACCATTCCCATACCATTAGTAAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGca                                                      <  2:140941/149‑1 (MQ=255)
                                                                                                 tGGTGAATAATTAACCATTCCCATACAATTAGTTAACAAAAAAGGGGGTATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGcaca                                                    <  1:126949/149‑1 (MQ=255)
                                                                                                   gTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTAATCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGAACATGTTAATAGATGGCGTGAAGcaccga                                                  >  2:563358/1‑146 (MQ=255)
                                                                                                      aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgt                                               <  1:316898/149‑1 (MQ=255)
                                                                                                      aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGCGGGATTTTATCTCCCCTTTAATTTTTCCACTATACTCGGCGTAGAATGTGGTGGAAACATCCCCATATACTGACGTTCATGTTAATAGAAGTAGTGATGAACAGtatt                                               >  2:665352/1‑146 (MQ=255)
                                                                                                      aataaataACCATTCCCATACAATTAGTAAACCAAAAAGGGGGGATTTTATCTCCCCTTAAATTTTTCCACTATTCACGGCTTTGAATGAGGGGGAAACAACCCCATATACTGACGTACATGTTAAAAGATGGCGTGAAGCACAGTCgt                                               >  1:48018/1‑149 (MQ=255)
                                                                                                       ataatTAAACATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTAAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATAAACTGACGTACATGTTAATAGATGGCGTGAAGCACAGGCgtg                                              <  1:24941/149‑1 (MQ=255)
                                                                                                         aatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTc                                            >  2:393536/1‑149 (MQ=255)
                                                                                                           ttAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCAt                                          >  1:705930/1‑149 (MQ=255)
                                                                                                            tAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATc                                         >  2:516119/1‑149 (MQ=255)
                                                                                                                cATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGAt                                     <  1:963071/149‑1 (MQ=255)
                                                                                                                   tCCCATACAATTAGTTAACCAAAAAGGGTGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAc                                  <  1:747951/149‑1 (MQ=255)
                                                                                                                     ccATACAATTAGTTACCCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCt                                <  1:1136501/149‑1 (MQ=255)
                                                                                                                     ccATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGAGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCt                                >  1:43619/1‑149 (MQ=255)
                                                                                                                      cATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGTAACATCCCCATATACTGACGTACATGTTTATAGATGGAGTGATGCACAGTCGTGTCTTCTGATTACCTg                               >  2:27946/1‑149 (MQ=255)
                                                                                                                         aCAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCg                            <  2:242090/149‑1 (MQ=255)
                                                                                                                             ttAGTTAACCAAACAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGAGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGaaa                        <  2:732070/149‑1 (MQ=255)
                                                                                                                             ttAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGaaa                        <  1:92015/149‑1 (MQ=255)
                                                                                                                                gTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTa                     >  1:646667/1‑149 (MQ=255)
                                                                                                                                    aCCAAAAAGCGGGGATATTATCTCCCCTTTAAATTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCAAATATACTGACGAACATGTTAATAGATCGCGTGAAGCACAGTCGTGTCATATAATTACCTGCAGCAAATTAAAca                 >  2:690356/1‑148 (MQ=255)
                                                                                                                                     ccAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGCCGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGGCATCTGATTACCTGGCGGAAATTAAACTa                >  2:855722/1‑149 (MQ=255)
                                                                                                                                      cAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCTTATAATGACGTACTTTTATATAGATGGCGTGAAGCACAGGCGTTTCATCTGATTACAATGCGGAAATTaaccgaa               >  2:644802/1‑144 (MQ=255)
                                                                                                                                         aaaGGGGAGATAAAATCCCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAagag            >  1:967195/1‑149 (MQ=255)
                                                                                                                                           aGGGGGGATTTTATCACCCCTTTAATTTTTACTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAagagag          <  2:991553/149‑1 (MQ=255)
                                                                                                                                                 gATTTTATCTCCGCTTTAATTTTTCCTCTATTCTGGGGGTTGAATGTGGGGGAAACATCCCCATATACTGACGTAGATTTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGTTTACCTGGCGTAAATTAAACTAAGAGAGAGctct    <  2:1023704/149‑1 (MQ=255)
                                                                                                                                                 gATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGctct    >  1:212131/1‑149 (MQ=255)
                                                                                                                                                 gATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGctct    <  2:510892/149‑1 (MQ=255)
                                                                                                                                                  aTTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTa   <  2:607411/149‑1 (MQ=255)
                                                                                                                                                  aTTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTa   >  2:261407/1‑149 (MQ=255)
                                                                                                                                                   ttttATCTCCCCTTTAATTTTTCGTCTTTTCTGGGCGTTGAATGTGGGGGCAACCTCTCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTAt  <  1:73069/149‑1 (MQ=255)
                                                                                                                                                   ttttATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTAt  <  2:728720/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTAT  >  NC_000913/458594‑458889

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