Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 453,656 T→A 10.2% S259T (TCC→ACC)  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
*CP009273453,6560TA10.2% 869.1 / 72.2 274S259T (TCC→ACC) clpXATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base T (132/114);  new base A (6/22);  total (138/136)
Fisher's exact test for biased strand distribution p-value = 1.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 277 total aligned reads displayed.

GGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTG  >  CP009273/453508‑453805
                                                                                                                                                    |                                                                                                                                                     
     tGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCAcc                                                                                                                                                 <  2:1835977/150‑1 (MQ=255)
       ttGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCgt                                                                                                                                               >  1:633501/1‑150 (MQ=255)
         gAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCgtgt                                                                                                                                             >  2:2914669/1‑150 (MQ=255)
         gAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCATGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCgtgt                                                                                                                                             >  2:2141934/1‑150 (MQ=255)
             cTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGaa                                                                                                                                         >  2:249162/1‑150 (MQ=255)
             cTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGaa                                                                                                                                         <  2:2102359/150‑1 (MQ=255)
              tGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGaaa                                                                                                                                        >  1:2965104/1‑150 (MQ=255)
               gATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAAc                                                                                                                                       >  2:3502087/1‑150 (MQ=255)
                aTCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAAcc                                                                                                                                      <  1:3373754/150‑1 (MQ=255)
                aTCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAAcc                                                                                                                                      >  1:309768/1‑150 (MQ=255)
                 tCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCg                                                                                                                                     >  2:1487990/1‑150 (MQ=255)
                 tCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCg                                                                                                                                     >  2:3756478/1‑150 (MQ=255)
                 tCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCg                                                                                                                                     <  1:4117884/150‑1 (MQ=255)
                 tCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTGTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCg                                                                                                                                     >  2:3756429/1‑150 (MQ=255)
                   gAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGAGACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGc                                                                                                                                   <  2:594480/150‑1 (MQ=255)
                     aGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTc                                                                                                                                 <  2:4113863/150‑1 (MQ=255)
                     aGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTc                                                                                                                                 >  2:388/1‑150 (MQ=255)
                      ggTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTcc                                                                                                                                <  2:606186/150‑1 (MQ=255)
                      ggTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTcc                                                                                                                                <  2:4218611/150‑1 (MQ=255)
                      ggTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAAATCTGGCAGGTTGATACCTCTAAGATCCTGTTTAAATGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTcc                                                                                                                                <  2:4218847/150‑1 (MQ=255)
                       gTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCg                                                                                                                               <  1:2581888/150‑1 (MQ=255)
                       gTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCg                                                                                                                               <  1:1729201/150‑1 (MQ=255)
                              aGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTgg                                                                                                                        <  2:378852/150‑1 (MQ=255)
                                ctgctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATCAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGtt                                                                                                                      >  2:3791655/1‑150 (MQ=255)
                                ctgctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGtt                                                                                                                      >  1:2074648/1‑150 (MQ=255)
                                     gTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGgc                                                                                                                 >  1:2719429/1‑150 (MQ=255)
                                     gTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGgc                                                                                                                 <  2:1900410/150‑1 (MQ=255)
                                        ccACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGcgcg                                                                                                              <  1:2385572/150‑1 (MQ=255)
                                        ccACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGcgcg                                                                                                              >  2:3253010/1‑150 (MQ=255)
                                         cACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGa                                                                                                             >  1:1173229/1‑150 (MQ=255)
                                             gCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGt                                                                                                         <  1:1017706/150‑1 (MQ=255)
                                              cAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTa                                                                                                        <  2:1302894/150‑1 (MQ=255)
                                                  gTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAg                                                                                                    <  2:4214075/150‑1 (MQ=255)
                                                   tGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGc                                                                                                   <  2:2548571/150‑1 (MQ=255)
                                                        gTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGt                                                                                              >  2:1795394/1‑150 (MQ=255)
                                                         tAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTc                                                                                             >  2:419700/1‑150 (MQ=255)
                                                          aaaCATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTcc                                                                                            <  1:1291105/150‑1 (MQ=255)
                                                           aaCATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCg                                                                                           >  1:1968505/1‑150 (MQ=255)
                                                            aCATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGa                                                                                          <  1:1017531/150‑1 (MQ=255)
                                                               tCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACaa                                                                                       <  2:4046968/150‑1 (MQ=255)
                                                                        ggAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGa                                                                              <  1:1305186/150‑1 (MQ=255)
                                                                          aaTTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAg                                                                            >  1:605014/1‑150 (MQ=255)
                                                                            ttCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGc                                                                          >  2:1949020/1‑150 (MQ=255)
                                                                             tCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCg                                                                         >  2:4352644/1‑150 (MQ=255)
                                                                              cTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGa                                                                        <  1:183323/150‑1 (MQ=255)
                                                                              cTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGa                                                                        <  1:183056/150‑1 (MQ=255)
                                                                               ttGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAg                                                                       >  2:3078839/1‑150 (MQ=255)
                                                                                tGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgc                                                                      <  1:2274037/150‑1 (MQ=255)
                                                                                  cAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgctg                                                                    <  2:118577/150‑1 (MQ=255)
                                                                                   aGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgctgc                                                                   >  2:1529463/1‑150 (MQ=255)
                                                                                    ggTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgctgct                                                                  >  1:4228564/1‑150 (MQ=255)
                                                                                    ggTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgctgct                                                                  <  2:4528011/150‑1 (MQ=255)
                                                                                    ggTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgctgct                                                                  <  2:2971017/150‑1 (MQ=255)
                                                                                    ggTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgctgct                                                                  >  1:3534909/1‑150 (MQ=255)
                                                                                    ggTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgctgct                                                                  >  1:4228530/1‑150 (MQ=255)
                                                                                       tGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgc                                                               >  2:1050956/1‑150 (MQ=255)
                                                                                       tGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgc                                                               >  1:695835/1‑150 (MQ=255)
                                                                                        gATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgcg                                                              >  2:3499700/1‑150 (MQ=255)
                                                                                        gATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgcg                                                              >  2:2685954/1‑150 (MQ=255)
                                                                                           aCCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgcgctg                                                           >  2:2435840/1‑148 (MQ=255)
                                                                                           aCCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAg                                                           <  2:2527867/150‑1 (MQ=255)
                                                                                               ctAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTg                                                       <  1:2903002/150‑1 (MQ=255)
                                                                                               ctAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTg                                                       <  1:617343/150‑1 (MQ=255)
                                                                                                tAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCCTTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGa                                                      >  2:1310775/1‑150 (MQ=255)
                                                                                                tAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGa                                                      >  1:547451/1‑150 (MQ=255)
                                                                                                  aGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAAc                                                    >  1:852543/1‑150 (MQ=255)
                                                                                                   gATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAAcc                                                   >  1:3500030/1‑150 (MQ=255)
                                                                                                        tGTTTATTTGTGGCGGTGCGTTTGCCGGTCTTGATCAACTGATTTCCCACCGTGTAGAAACAGGCTCCGGAATTTGTTTTTGCGCGACTGTAAAAGCGAAGTCCGACAAAGCAAGCGTAGTCGATCTTCTGGCGCAGGTTGAACCGGAAg                                              <  1:2874832/150‑1 (MQ=255)
                                                                                                        tGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAg                                              <  2:644331/150‑1 (MQ=255)
                                                                                                        tGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAg                                              >  2:563644/1‑150 (MQ=255)
                                                                                                        tGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAg                                              >  1:79174/1‑150 (MQ=255)
                                                                                                          tttATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGAt                                            <  2:473909/150‑1 (MQ=255)
                                                                                                            tATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCt                                          >  2:2654095/1‑150 (MQ=255)
                                                                                                             aTTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTg                                         >  1:3964386/1‑150 (MQ=255)
                                                                                                              tttGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGa                                        >  2:657669/1‑150 (MQ=255)
                                                                                                                  tgGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCaa                                    <  1:1897040/150‑1 (MQ=255)
                                                                                                                    gCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGt                                  <  1:2816516/150‑1 (MQ=255)
                                                                                                                     cGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGtt                                 <  1:765734/150‑1 (MQ=255)
                                                                                                                      ggTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCATGttt                                >  2:1254532/1‑150 (MQ=255)
                                                                                                                      ggTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGttt                                >  1:4375545/1‑150 (MQ=255)
                                                                                                                       gTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTg                               <  1:3000624/150‑1 (MQ=255)
                                                                                                                        tGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTgg                              >  1:1042416/1‑150 (MQ=255)
                                                                                                                        tGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTgg                              <  1:498103/150‑1 (MQ=255)
                                                                                                                        tGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTgg                              <  1:498316/150‑1 (MQ=255)
                                                                                                                         gCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGATGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGt                             <  1:260959/150‑1 (MQ=255)
                                                                                                                           gTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCt                           <  2:1777912/150‑1 (MQ=255)
                                                                                                                             ttGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGACGATCTGATCAAGTTTGGTCTTa                         >  2:3326085/1‑150 (MQ=255)
                                                                                                                             ttGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTa                         <  2:1414396/150‑1 (MQ=255)
                                                                                                                             ttGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTa                         >  2:743421/1‑150 (MQ=255)
                                                                                                                             ttGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTa                         <  1:1465991/150‑1 (MQ=255)
                                                                                                                                 cGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATccc                     >  1:4170740/1‑150 (MQ=255)
                                                                                                                                  ggTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCt                    <  1:3338855/150‑1 (MQ=255)
                                                                                                                                  ggTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCt                    <  1:4442903/150‑1 (MQ=255)
                                                                                                                                      tGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGt                >  2:2275326/1‑150 (MQ=255)
                                                                                                                                       ggATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGTGtt               >  2:1645331/1‑150 (MQ=255)
                                                                                                                                        gATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGttt              >  1:4186733/1‑150 (MQ=255)
                                                                                                                                        gATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGttt              <  1:4352483/150‑1 (MQ=255)
                                                                                                                                          tAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTAt            <  2:2880278/150‑1 (MQ=255)
                                                                                                                                            aaGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTg          <  1:3909279/150‑1 (MQ=255)
                                                                                                                                                gATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGgtcg      >  1:2143460/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTG  >  CP009273/453508‑453805

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