Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 458,189 A→G 15.9% D255G (GAT→GGT)  clpX → ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913458,1890AG15.9% 173.3 / 4.6 68D255G (GAT→GGT) clpXATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base A (5/49);  new base G (9/2);  total (17/51)
Fisher's exact test for biased strand distribution p-value = 2.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGA  >  NC_000913/457952‑458425
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
cacaGGCTGGTATTGTCTTCATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGAt                                                                                                                                                                                                                                               <  1:994742/239‑1 (MQ=255)
    gggtttttttGTCTACATCGATGAAACCGACCAGATTTTTGGCAAGCCAGAACACCCGTCCATGACCCGAGGCGTTGCCGGTGAAGGCGTACGGCCGGCCATGGTGGAGCGGATCGAAGGTACGGTAGCTGCTGTTCCACAGCAAGGTGGGCGTAAACATCCGGAGCAGGAATTCTTGCAGGTTCATAACTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAg                                                                                                                                                                                                                                           <  2:995581/231‑1 (MQ=255)
    agagacagttGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAg                                                                                                                                                                                                                                           <  1:606801/231‑1 (MQ=255)
      gTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAg                                                                                                                                                                                                                                           <  1:189820/237‑1 (MQ=255)
              gtctcaaTCGATGAGATCGACATAGTTTCGCGAAAGTCAGACAACCCGTCCATTTCCCTGGGCGTTTCCGGTGAAGGCGTCCCGCAGGCACTGTTGAAACTGAGCGAAGGTACGGTCGCTGCTGTTCCACCGCAAGGGGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCa                                                                                                                                                                                                                                 <  2:1289911/233‑1 (MQ=255)
              gTCTAAATCGATGAAAGCGACAAGGTTTCGCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGGTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCa                                                                                                                                                                                                                                 <  1:786145/239‑1 (MQ=255)
                cTACATCGATGAAAGCGACAAGATTCTCGGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCAc                                                                                                                                                                                                                                <  1:1244916/238‑1 (MQ=255)
                cTAAATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCAcc                                                                                                                                                                                                                               <  1:1123031/239‑1 (MQ=255)
                 tACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTTGCGGTGCGGTTTCCGGTCTGGTTAAAAGGATTTTTCAccc                                                                                                                                                                                                                              >  1:1235085/1‑238 (MQ=255)
                    aTCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGCTTTTCTGTCTGTCTAAAGTGAGTTTCCAcccgct                                                                                                                                                                                                                           >  1:970279/1‑235 (MQ=255)
                         tGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTCTGGGCTTGGTAAAATGATTTcccctgtctctttt                                                                                                                                                                                                                       >  1:1015161/1‑227 (MQ=255)
                         tGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGTAAAGAGATTTCCAACCctgtgggaa                                                                                                                                                                                                                      >  1:1139855/1‑234 (MQ=255)
                                                 aGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGCCGGTCTGGGTAAAAGGATTTTCCCCCGTGGAGAAAACGGCCCCGGCCTTTGTTTTggg                                                                                                                                                                                              >  1:537281/1‑238 (MQ=255)
                                                    caacacaCCCGTCCATTACCCGAGGCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGgcgc                                                                                                                                                                                            <  1:1002558/232‑1 (MQ=255)
                                                      gACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGGCGGGTTAAATGGATTTCTCCCCGGGTGGCAAACGCACCCGGCCTTGGgttttcacagcg                                                                                                                                                                                         >  1:80187/1‑231 (MQ=255)
                                                        cAACCCGTCCATTGCCCGATCCGGTTCCGGTGAAGGCGTACCGCAGGCACGGTTGAACCTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTATGGATAAAGTGATTTCCCTCCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACgg                                                                                                                                                                                       <  1:714740/239‑1 (MQ=255)
                                                                 cATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATaaaa                                                                                                                                                                                                                                           >  1:517912/1‑178 (MQ=255)
                                                                 cATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAg                                                                                                                                                                                                                                           <  2:517912/179‑1 (MQ=255)
                                                                  aTTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGaa                                                                                                                                                                             <  1:949766/239‑1 (MQ=255)
                                                                          agaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGGTTTTCCGGGCTTGGTAAAGTGGTTTTCCCCCCTGGAGGAACCGGCCCCCGCATTGGGTTTGGTGCGACGGTTAAACGCAAGTcccgca                                                                                                                                                                     >  1:781949/1‑235 (MQ=255)
                                                                          agaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTTCGGTCTGGGTAAAAGGATTTCCCCCCGTGTAGAAACCCGCTCCGGCCTTGGTTTTGGCGCGACGGTCAAAGCGAAGTCCGACa                                                                                                                                                                     >  1:203872/1‑239 (MQ=255)
                                                                                                         cTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTCTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGt                                                                                                                                      <  1:996759/239‑1 (MQ=255)
                                                                                                           gTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTTCGTTTGCCGGACTGGGTAAAGTGCATTTCCACCGTGTAGAAACCGGCTCCGGCATTTGGTTTTGGGGGCCGGGAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGGCGCTGGCGCCGGTTg                                                                                                                                    >  1:201529/1‑239 (MQ=255)
                                                                                                           gTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGGTTTGCCGGTCTTGATAAAGGGATTTTCCACCGTGTAGAAAACGGCTCCGGCATTGGTTTTTGCGCGGCGGGAAAAGAGAAGGCCGGCAAAACAAACGAAGGCGAGCTGCTGGGCCAGGGTg                                                                                                                                    >  1:686838/1‑239 (MQ=255)
                                                                                                           gTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTTGCGGGCTGGATAAAATGATTTTCCCCCGCGTAGAAAACAGGCCCGGCATTTGGTTTGGGGCGGCGGTGAAAGCGAAGAGCGCCAAAAAAAGCGAAGGGGGGCTGCTggggggggggg                                                                                                                                    >  2:1196546/1‑230 (MQ=255)
                                                                                                           gTTGAAACTGAACGAAGGTAGGTGAGTGGCTGTTCCACCGTAAGGTGGGCGTAAACAACCCGAGAAGGAATTCTTGGAGGTCGATACCTGTAAGATCCTGTTTATTTGGGGCGGTGCGGTTGCCGGTCTGGATAAAGTGAGTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCTCAGGTTg                                                                                                                                    <  2:201529/239‑1 (MQ=255)
                                                                                                            ttGAAACTGATCGAGGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAGCATCCGCAGCAGGAATGCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGa                                                                                                                                   <  1:197806/239‑1 (MQ=255)
                                                                                                            ttGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGGTTTGCCGGGCTGGGTAAAATGATTTTCCACCGTGTCGAAAACCGCCCCGGCATTGGGTTTGGTGCGCCGGGAAAAGAGAAGACCGGCAAAGCAAGCGAAGGCGGCCTGCTGTGGCCGGTTGa                                                                                                                                   >  1:1074885/1‑239 (MQ=255)
                                                                                                            ttGAAACTGACTGAAGGTGCGTGGTCTGCTGTGCCACCGCACGGGTGGCGTTAACATCCGCATCAGGAATTCGTGAGGGTTGATACCGTTAAGCTCCTGTTGATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGGGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGa                                                                                                                                   <  2:781949/239‑1 (MQ=255)
                                                                                                                 ccTGATCGAAGGGACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCgg                                                                                                                              <  1:170455/238‑1 (MQ=255)
                                                                                                                      tCGAAGGTACGGTAGCTGCGGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGAt                                                                                                                         <  1:537719/239‑1 (MQ=255)
                                                                                                                       cGAAGGTTCGGTAGGTGCTGTCCCCCCGCAAGGTTGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGAt                                                                                                                         <  1:511365/238‑1 (MQ=255)
                                                                                                                             tACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgcg                                                                                                                                           <  1:922143/214‑1 (MQ=255)
                                                                                                                                 gtcgctcttgTTCCCCCGCAATTGTGGCGTAAAATTCCGTAGCAGGAATTCTTGACGGTTGATACCTCTAAGGTCCTGTTGTTTTGTGGCGGTGCGTTTGCCGGTCGGGATAAAGTGATTTCCCACCGGTTAGAAACCGGCTCCGGCATTGGTTTTGGCGGGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGtt                                                                                                              <  2:644447/231‑1 (MQ=255)
                                                                                                                                  tAGCTGCTGTTCCACCGGAAGGTGGGCGTAAACCTCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGtt                                                                                                              <  1:1094244/238‑1 (MQ=255)
                                                                                                                                  tAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACCTCCGCAGCAGGAATTCTTGCAGGTTTATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGtt                                                                                                              <  1:1205243/238‑1 (MQ=255)
                                                                                                                                   ggctgcgGTTCCCCCTGAAGGGGGGGGTAAACATCCGAAGCAGGAATTCTGTCAGGTTGATACCTCGAAGATCCTGTTTAGTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGttt                                                                                                             <  2:834255/237‑1 (MQ=255)
                                                                                                                                       gctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGTAAAGTg                                                                                                                                                                                                                                         >  1:32436/1‑111 (MQ=255)
                                                                                                                                       gctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTg                                                                                                                                                                                                                                         <  2:32436/111‑1 (MQ=255)
                                                                                                                                       gctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCt                                                                                                        <  1:1196546/239‑1 (MQ=255)
                                                                                                                                             ccACCGGAATGTGGGCTGAACCATCCGCAGCAGGAAATCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATccc                                                                                                  <  1:870587/239‑1 (MQ=255)
                                                                                                                                             ccACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATccc                                                                                                  <  1:568557/239‑1 (MQ=255)
                                                                                                                                              cccctgaaGGTGGGTATAAACCTCCGGAGCAGTAATTCTTGCAGTTTCGTCCCTCTAAGAGCCTGTTTGTTTGTGGCGGCGCGTTTGGCGGGCCGGATAAAGTGAGTGCCCACCGTGTAGAAACCGGCTTCTGCCTTTGTTTGGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGACGAGCTGCTGGCGCAGGTTGAACCGGAATATCTGATCAAGTTTGGTCTTATCCCt                                                                                                 <  2:313096/233‑1 (MQ=255)
                                                                                                                                              ccccGTAAGGTGGGCTTAACCATCCGGAGCAGGAATTCTTGCAGGTTGATTCCTCTAAGATCCTGTCTATTTGTGGCGGTGCGGTGTCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCt                                                                                                 <  2:260982/237‑1 (MQ=255)
                                                                                                                                               gccGCAAGCGGGGCGTAAACATCCGGAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCAGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTg                                                                                                <  1:177259/238‑1 (MQ=255)
                                                                                                                                                ccGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGtt                                                                                                              <  1:539900/224‑1 (MQ=255)
                                                                                                                                                                            gggAATTTTTGCAGGTTGATACCTCTAAGATCATGTTTAATTGGTGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCa                                                                   <  1:958429/238‑1 (MQ=255)
                                                                                                                                                                            agGAATTCTTGCAGGTTGATTCCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCg                                                                     <  1:382926/237‑1 (MQ=255)
                                                                                                                                                                               aaTTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATATAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAAc                                                                 <  1:983329/238‑1 (MQ=255)
                                                                                                                                                                               aaTTCTTGAAGGTTGATTCCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCCCCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGATAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCa                                                                   <  1:673463/236‑1 (MQ=255)
                                                                                                                                                                                    ttGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgcg                                                                                                                                           <  1:1208849/159‑1 (MQ=255)
                                                                                                                                                                                           ttGATACCTCTAAGGTCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTtgaatgaa                                                       <  1:823613/236‑1 (MQ=255)
                                                                                                                                                                                                  ctctAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTc                                                                                                         <  1:168042/179‑1 (MQ=255)
                                                                                                                                                                                                  ctctAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCgaag                                            <  1:857284/239‑1 (MQ=255)
                                                                                                                                                                                                   tctAAGGTCCTGTTTATTTGGGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCgaaga                                           <  1:974074/239‑1 (MQ=255)
                                                                                                                                                                                                         aTCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGa                                                                                                                                                                                                                                                <  1:942758/38‑1 (MQ=255)
                                                                                                                                                                                                         aTCATGTTTATTGGGTGCGGTGCGTTTGCCGGTCTGGATAAAGGTATTACCGACCGTGTAGAAACCGGCTACGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGctct                                     <  1:435958/239‑1 (MQ=255)
                                                                                                                                                                                                            cTGTTTATTTGTGGCGGTGCGTTTTCCGGTCTGGGTAAAGTGATTTTCCACCGTGTAGAAAACGGCTCCGGCATTGGGTTTGGCGCGACGGGAAAAACGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGGCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTGTTGGCAACGTTGAATTAACTGAGCGAGGAAGCTCTGAt                                  >  1:326541/1‑239 (MQ=255)
                                                                                                                                                                                                             tGTTTATTTGGGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGAtt                                 <  1:1105375/239‑1 (MQ=255)
                                                                                                                                                                                                               tttATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCa                               <  1:910734/239‑1 (MQ=255)
                                                                                                                                                                                                               tttATTTGGGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCa                               <  1:1156287/239‑1 (MQ=255)
                                                                                                                                                                                                                       tgGCGGTGCGTTTGCCGGTCTGGATAAAGTGCTTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGGCGGACGGGAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATATGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCa                       <  2:203872/239‑1 (MQ=255)
                                                                                                                                                                                                                       ggggggTGCGTTTGGCGGTCGGGAGAAAGGGAGTTCCCACGGTGTAGAAACCGGCTCAGGCATTGGGTTTGGCCCGGGCGTAAAAGCGAAGTCCGAAAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTTTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCa                       <  2:1074885/235‑1 (MQ=255)
                                                                                                                                                                                                                         ggggTGCTTTTTCTGGGCTGAGTAAAGGTATTTCACACCGTGTAGAATCCGGCTTGGGCATTGGTGTCGGCGCGACGGTAAAAGCGAAGTCCGTCAAAGCAAGCGAAGGTGAGATGCTGTCGCAGGTTGCACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGGTGAATGAACTGAGCGAAGAAGCGCTGATTCGGACCATCaaa                     <  2:210137/237‑1 (MQ=255)
                                                                                                                                                                                                                                       ggTCTGGGTAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGc       <  1:260512/239‑1 (MQ=255)
                                                                                                                                                                                                                                         tCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAACGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGAGCAAGTTTGGTCTTATCCCTTAGTTTATTGGTCTACTTCCGGTTGTTGCCAAGTTGCAGCAACTGAGCGAAGAAGCTCTGATTCACATCCTAAAAGAGCCGCAAACCGcca     >  2:364237/1‑238 (MQ=255)
                                                                                                                                                                                                                                          cTGGAAAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGccc     <  1:1093492/238‑1 (MQ=255)
                                                                                                                                                                                                                                           tGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGCAGGCGAGCTGCTGGCGCAGGGTGAACCGGAAGATCTGATCAAGTTTGGGCTAATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAGTGAAATAGCGGAATAAGCCCGGGTTCAGAACCCAAAACAGACAAAAAACGcccaa   >  2:685890/1‑237 (MQ=255)
                                                                                                                                                                                                                                            ggATAAAGTGGTTTCCCTCCGGGTAGAAAGCGGGTCCGGCATTGGTTTTGGCGCGTCGGTAAAAGCGAAGTGCGAAAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGACGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACTCCCTGa  <  1:428283/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
CACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGA  >  NC_000913/457952‑458425

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