Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 458,195 T→G 13.4% V257G (GTG→GGG)  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,1950TG13.4% 166.0 / 2.0 71V257G (GTG→GGG) clpXATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base T (12/46);  new base G (6/5);  total (20/51)
Fisher's exact test for biased strand distribution p-value = 2.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAG  >  NC_000913/457966‑458430
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
gtctcaaTCGATGAGATCGACATAGTTTCGCGAAAGTCAGACAACCCGTCCATTTCCCTGGGCGTTTCCGGTGAAGGCGTCCCGCAGGCACTGTTGAAACTGAGCGAAGGTACGGTCGCTGCTGTTCCACCGCAAGGGGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCa                                                                                                                                                                                                                                       <  1:1289911/233‑1 (MQ=255)
gTCTAAATCGATGAAAGCGACAAGGTTTCGCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGGTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCa                                                                                                                                                                                                                                       <  2:786145/239‑1 (MQ=255)
  cTACATCGATGAAAGCGACAAGATTCTCGGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCAc                                                                                                                                                                                                                                      <  2:1244916/238‑1 (MQ=255)
  cTAAATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCAcc                                                                                                                                                                                                                                     <  2:1123031/239‑1 (MQ=255)
   tACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTTGCGGTGCGGTTTCCGGTCTGGTTAAAAGGATTTTTCAccc                                                                                                                                                                                                                                    >  2:1235085/1‑238 (MQ=255)
      aTCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGCTTTTCTGTCTGTCTAAAGTGAGTTTCCAcccgct                                                                                                                                                                                                                                 >  2:970279/1‑235 (MQ=255)
           tGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTCTGGGCTTGGTAAAATGATTTcccctgtctctttt                                                                                                                                                                                                                             >  2:1015161/1‑227 (MQ=255)
           tGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGTAAAGAGATTTCCAACCctgtgggaa                                                                                                                                                                                                                            >  2:1139855/1‑234 (MQ=255)
                                   aGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGCCGGTCTGGGTAAAAGGATTTTCCCCCGTGGAGAAAACGGCCCCGGCCTTTGTTTTggg                                                                                                                                                                                                   >  2:537281/1‑238 (MQ=255)
                                      caacacaCCCGTCCATTACCCGAGGCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGgcgc                                                                                                                                                                                                 <  2:1002558/232‑1 (MQ=255)
                                        gACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGGCGGGTTAAATGGATTTCTCCCCGGGTGGCAAACGCACCCGGCCTTGGgttttcacagcg                                                                                                                                                                                              >  2:80187/1‑231 (MQ=255)
                                          cAACCCGTCCATTGCCCGATCCGGTTCCGGTGAAGGCGTACCGCAGGCACGGTTGAACCTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTATGGATAAAGTGATTTCCCTCCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACgg                                                                                                                                                                                            <  2:714740/239‑1 (MQ=255)
                                                    aTTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGaa                                                                                                                                                                                  <  2:949766/239‑1 (MQ=255)
                                                           gagaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGGGTTTCCGGGCTGGATAAAGTGAGTTTCACCCAGGTAGAAAACGGGCCCCGCCCTTGGTTTTTGCGGACGGTcaaacgcaatgctcgc                                                                                                                                                                           >  2:1070363/1‑227 (MQ=255)
                                                            agaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGGTTTTCCGGGCTTGGTAAAGTGGTTTTCCCCCCTGGAGGAACCGGCCCCCGCATTGGGTTTGGTGCGACGGTTAAACGCAAGTcccgca                                                                                                                                                                          >  2:781949/1‑235 (MQ=255)
                                                            agaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTTCGGTCTGGGTAAAAGGATTTCCCCCCGTGTAGAAACCCGCTCCGGCCTTGGTTTTGGCGCGACGGTCAAAGCGAAGTCCGACa                                                                                                                                                                          >  2:203872/1‑239 (MQ=255)
                                                                        tGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTTTGGCGGTTGGCTTGTCGGGCTCGGTCAAGCGGTTTTCCCCCCTGTTGGAAACGGCCCCTGCACTTGCTTTGGTCGCGCCGGAAAAAGCaaaagcgacaaacaaacgcgaa                                                                                                                                                               >  1:899939/1‑220 (MQ=255)
                                                                                           cTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTCTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGt                                                                                                                                           <  2:996759/239‑1 (MQ=255)
                                                                                             gTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTTCGTTTGCCGGACTGGGTAAAGTGCATTTCCACCGTGTAGAAACCGGCTCCGGCATTTGGTTTTGGGGGCCGGGAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGGCGCTGGCGCCGGTTg                                                                                                                                         >  2:201529/1‑239 (MQ=255)
                                                                                             gTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGGTTTGCCGGTCTTGATAAAGGGATTTTCCACCGTGTAGAAAACGGCTCCGGCATTGGTTTTTGCGCGGCGGGAAAAGAGAAGGCCGGCAAAACAAACGAAGGCGAGCTGCTGGGCCAGGGTg                                                                                                                                         >  2:686838/1‑239 (MQ=255)
                                                                                             gTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTTGCGGGCTGGATAAAATGATTTTCCCCCGCGTAGAAAACAGGCCCGGCATTTGGTTTGGGGCGGCGGTGAAAGCGAAGAGCGCCAAAAAAAGCGAAGGGGGGCTGCTggggggggggg                                                                                                                                         >  1:1196546/1‑230 (MQ=255)
                                                                                             gTTGAAACTGAACGAAGGTAGGTGAGTGGCTGTTCCACCGTAAGGTGGGCGTAAACAACCCGAGAAGGAATTCTTGGAGGTCGATACCTGTAAGATCCTGTTTATTTGGGGCGGTGCGGTTGCCGGTCTGGATAAAGTGAGTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCTCAGGTTg                                                                                                                                         <  1:201529/239‑1 (MQ=255)
                                                                                              tttgatgggctggggggTTCGGTGGCGTCTGTTCCCCCGAAAGGTGGGCCTAAAACTCCGCAGCAGGAATTCTTGCAGGTTGATTCCTCTTAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTTGATAAAGTGATTTCCCCCCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGa                                                                                                                                        <  1:760013/224‑1 (MQ=255)
                                                                                              ttGAAACTGATCGAGGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAGCATCCGCAGCAGGAATGCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGa                                                                                                                                        <  2:197806/239‑1 (MQ=255)
                                                                                              ttGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGGTTTGCCGGGCTGGGTAAAATGATTTTCCACCGTGTCGAAAACCGCCCCGGCATTGGGTTTGGTGCGCCGGGAAAAGAGAAGACCGGCAAAGCAAGCGAAGGCGGCCTGCTGTGGCCGGTTGa                                                                                                                                        >  2:1074885/1‑239 (MQ=255)
                                                                                              ttGAAACTGACTGAAGGTGCGTGGTCTGCTGTGCCACCGCACGGGTGGCGTTAACATCCGCATCAGGAATTCGTGAGGGTTGATACCGTTAAGCTCCTGTTGATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGGGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGa                                                                                                                                        <  1:781949/239‑1 (MQ=255)
                                                                                                   ccTGATCGAAGGGACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCgg                                                                                                                                   <  2:170455/238‑1 (MQ=255)
                                                                                                        tCGAAGGTACGGTAGCTGCGGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGAt                                                                                                                              <  2:537719/239‑1 (MQ=255)
                                                                                                         cGAAGGTTCGGTAGGTGCTGTCCCCCCGCAAGGTTGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGAt                                                                                                                              <  2:511365/238‑1 (MQ=255)
                                                                                                             tgttcggtggttgtttttCCCCCGCATGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTTTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGAt                                                                                                                         <  1:537281/223‑1 (MQ=255)
                                                                                                               tACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgcg                                                                                                                                                <  2:922143/214‑1 (MQ=255)
                                                                                                                 cGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGACCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGATAAATGATTTCTCCCCGTGTAGAAAACGTCTCCGGCGTTGGGTTTGGAGCGCGGGCAAAATAGAAAGCCGGCAAAGAAAGCGCAAGCGGGCGGGTGGGGCAGAGTGAAAAGGaaaaaatgaacaaa                                                                                                                     >  1:27905/1‑227 (MQ=255)
                                                                                                                   gtcgctcttgTTCCCCCGCAATTGTGGCGTAAAATTCCGTAGCAGGAATTCTTGACGGTTGATACCTCTAAGGTCCTGTTGTTTTGTGGCGGTGCGTTTGCCGGTCGGGATAAAGTGATTTCCCACCGGTTAGAAACCGGCTCCGGCATTGGTTTTGGCGGGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGtt                                                                                                                   <  1:644447/231‑1 (MQ=255)
                                                                                                                    tAGCTGCTGTTCCACCGGAAGGTGGGCGTAAACCTCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGtt                                                                                                                   <  2:1094244/238‑1 (MQ=255)
                                                                                                                    tAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACCTCCGCAGCAGGAATTCTTGCAGGTTTATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGtt                                                                                                                   <  2:1205243/238‑1 (MQ=255)
                                                                                                                     ggctgcgGTTCCCCCTGAAGGGGGGGGTAAACATCCGAAGCAGGAATTCTGTCAGGTTGATACCTCGAAGATCCTGTTTAGTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGttt                                                                                                                  <  1:834255/237‑1 (MQ=255)
                                                                                                                         gctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGTAAAGTg                                                                                                                                                                                                                                               >  2:32436/1‑111 (MQ=255)
                                                                                                                         gctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTg                                                                                                                                                                                                                                               <  1:32436/111‑1 (MQ=255)
                                                                                                                         gctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCt                                                                                                             <  2:1196546/239‑1 (MQ=255)
                                                                                                                               ccACCGGAATGTGGGCTGAACCATCCGCAGCAGGAAATCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATccc                                                                                                       <  2:870587/239‑1 (MQ=255)
                                                                                                                               ccACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATccc                                                                                                       <  2:568557/239‑1 (MQ=255)
                                                                                                                                cccctgaaGGTGGGTATAAACCTCCGGAGCAGTAATTCTTGCAGTTTCGTCCCTCTAAGAGCCTGTTTGTTTGTGGCGGCGCGTTTGGCGGGCCGGATAAAGTGAGTGCCCACCGTGTAGAAACCGGCTTCTGCCTTTGTTTGGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGACGAGCTGCTGGCGCAGGTTGAACCGGAATATCTGATCAAGTTTGGTCTTATCCCt                                                                                                      <  1:313096/233‑1 (MQ=255)
                                                                                                                                ccccGTAAGGTGGGCTTAACCATCCGGAGCAGGAATTCTTGCAGGTTGATTCCTCTAAGATCCTGTCTATTTGTGGCGGTGCGGTGTCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCt                                                                                                      <  1:260982/237‑1 (MQ=255)
                                                                                                                                 gccGCAAGCGGGGCGTAAACATCCGGAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCAGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTg                                                                                                     <  2:177259/238‑1 (MQ=255)
                                                                                                                                  ccGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGtt                                                                                                                   <  2:539900/224‑1 (MQ=255)
                                                                                                                                                              gggAATTTTTGCAGGTTGATACCTCTAAGATCATGTTTAATTGGTGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCa                                                                        <  2:958429/238‑1 (MQ=255)
                                                                                                                                                              agGAATTCTTGCAGGTTGATTCCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCg                                                                          <  2:382926/237‑1 (MQ=255)
                                                                                                                                                                 aaTTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATATAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAAc                                                                      <  2:983329/238‑1 (MQ=255)
                                                                                                                                                                 aaTTCTTGAAGGTTGATTCCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCCCCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGATAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCa                                                                        <  2:673463/236‑1 (MQ=255)
                                                                                                                                                                      ttGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgcg                                                                                                                                                <  2:1208849/159‑1 (MQ=255)
                                                                                                                                                                           ggTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGCCGGTCTGGATAAAAGGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGGTTTTGCGCGACGGTAAAAGCGAAGTCCGCCAAAACAAGCGAAGGCGAGCGGCTGGCGCAGGGTGCACCGGGAGCACCGAGACAGCTTGGTACTCTCCCTGAGTTGTTTTGGCTTTTTCCGGTTGtgcggacagtgtttgaaa                                                           >  1:1121084/1‑222 (MQ=255)
                                                                                                                                                                             ttGATACCTCTAAGGTCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTtgaatgaa                                                            <  2:823613/236‑1 (MQ=255)
                                                                                                                                                                                    ctctAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTc                                                                                                              <  2:168042/179‑1 (MQ=255)
                                                                                                                                                                                    ctctAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCgaag                                                 <  2:857284/239‑1 (MQ=255)
                                                                                                                                                                                     tctAAGGTCCTGTTTATTTGGGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCgaaga                                                <  2:974074/239‑1 (MQ=255)
                                                                                                                                                                                           aTCATGTTTATTGGGTGCGGTGCGTTTGCCGGTCTGGATAAAGGTATTACCGACCGTGTAGAAACCGGCTACGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGctct                                          <  2:435958/239‑1 (MQ=255)
                                                                                                                                                                                              cTGTTTATTTGTGGCGGTGCGTTTTCCGGTCTGGGTAAAGTGATTTTCCACCGTGTAGAAAACGGCTCCGGCATTGGGTTTGGCGCGACGGGAAAAACGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGGCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTGTTGGCAACGTTGAATTAACTGAGCGAGGAAGCTCTGAt                                       >  2:326541/1‑239 (MQ=255)
                                                                                                                                                                                               tGTTTATTTGGGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGAtt                                      <  2:1105375/239‑1 (MQ=255)
                                                                                                                                                                                                 tttATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCa                                    <  2:910734/239‑1 (MQ=255)
                                                                                                                                                                                                 tttATTTGGGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCa                                    <  2:1156287/239‑1 (MQ=255)
                                                                                                                                                                                                         tgGCGGTGCGTTTGCCGGTCTGGATAAAGTGCTTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGGCGGACGGGAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATATGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCa                            <  1:203872/239‑1 (MQ=255)
                                                                                                                                                                                                         ggggggTGCGTTTGGCGGTCGGGAGAAAGGGAGTTCCCACGGTGTAGAAACCGGCTCAGGCATTGGGTTTGGCCCGGGCGTAAAAGCGAAGTCCGAAAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTTTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCa                            <  1:1074885/235‑1 (MQ=255)
                                                                                                                                                                                                           ggggTGCTTTTTCTGGGCTGAGTAAAGGTATTTCACACCGTGTAGAATCCGGCTTGGGCATTGGTGTCGGCGCGACGGTAAAAGCGAAGTCCGTCAAAGCAAGCGAAGGTGAGATGCTGTCGCAGGTTGCACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGGTGAATGAACTGAGCGAAGAAGCGCTGATTCGGACCATCaaa                          <  1:210137/237‑1 (MQ=255)
                                                                                                                                                                                                             ggggtgttttgcggtttttgATAAAGGGATTGCCCACCGTGTAGAAACCGGCTCCGGCATGGGTTTTGGCGCGCATGTAAAAGCGAAGGCCGACAAAGCAAGCGAAGGCGAGCTTCTGGCTCAGGTTGAACCTGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCaaa                          <  1:686838/219‑1 (MQ=255)
                                                                                                                                                                                                                         ggTCTGGGTAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGc            <  2:260512/239‑1 (MQ=255)
                                                                                                                                                                                                                           tCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAACGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGAGCAAGTTTGGTCTTATCCCTTAGTTTATTGGTCTACTTCCGGTTGTTGCCAAGTTGCAGCAACTGAGCGAAGAAGCTCTGATTCACATCCTAAAAGAGCCGCAAACCGcca          >  1:364237/1‑238 (MQ=255)
                                                                                                                                                                                                                            cTGGAAAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGccc          <  2:1093492/238‑1 (MQ=255)
                                                                                                                                                                                                                             tGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGCAGGCGAGCTGCTGGCGCAGGGTGAACCGGAAGATCTGATCAAGTTTGGGCTAATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAGTGAAATAGCGGAATAAGCCCGGGTTCAGAACCCAAAACAGACAAAAAACGcccaa        >  1:685890/1‑237 (MQ=255)
                                                                                                                                                                                                                               ggATAAAGTGGTTTCCCTCCGGGTAGAAAGCGGGTCCGGCATTGGTTTTGGCGCGTCGGTAAAAGCGAAGTGCGAAAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGACGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACTCCCTGa       <  2:428283/239‑1 (MQ=255)
                                                                                                                                                                                                                                    aaGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAg  <  2:933306/239‑1 (MQ=255)
                                                                                                                                                                                                                                    aaGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTGTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAg  <  2:1066597/239‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
GTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAG  >  NC_000913/457966‑458430

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