Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,354,467 T→A 5.3% P131P (CCA→CCT)  dsbD ← fused thiol:disulfide interchange protein: activator of DsbC/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,354,4670TA 5.3% 645.4 / 20.2 209P131P (CCA→CCT) dsbDfused thiol:disulfide interchange protein: activator of DsbC/conserved protein
Reads supporting (aligned to +/- strand):  ref base T (79/119);  new base A (5/6);  total (84/125)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

Only 100 of 210 total aligned reads displayed.

CGTAAAGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAG  >  CP009273/4354320‑4354611
                                                                                                                                                   |                                                                                                                                                
cGTAAAGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCt                                                                                                                                                  >  2:2330460/1‑148 (MQ=255)
  tAAAGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTgg                                                                                                                                                >  1:11245/1‑148 (MQ=255)
  tAAAGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTgg                                                                                                                                                >  1:1900724/1‑148 (MQ=255)
  tAAAGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTgg                                                                                                                                                >  2:1984559/1‑148 (MQ=255)
   aaaGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGc                                                                                                                                               >  1:1214113/1‑148 (MQ=255)
   aaaGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGc                                                                                                                                               >  1:2965236/1‑148 (MQ=255)
    aaGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCg                                                                                                                                              >  2:3190833/1‑148 (MQ=255)
    aaGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCg                                                                                                                                              >  2:1244291/1‑148 (MQ=255)
     aGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAAGTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCgg                                                                                                                                             <  1:1850046/148‑1 (MQ=255)
      ggCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGa                                                                                                                                            <  2:184561/148‑1 (MQ=255)
        cGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATa                                                                                                                                          <  1:1797158/148‑1 (MQ=255)
        cGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATa                                                                                                                                          <  1:968441/148‑1 (MQ=255)
         gATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATaa                                                                                                                                         <  2:1353083/148‑1 (MQ=255)
           tACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACa                                                                                                                                       <  1:1155944/148‑1 (MQ=255)
            aCCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAg                                                                                                                                      >  2:2743052/1‑148 (MQ=255)
             ccAAGACCGATCAACAACGCCCAGAGCGCGGAAACGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGa                                                                                                                                     <  2:2261678/148‑1 (MQ=255)
               aaTACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGaaa                                                                                                                                   >  2:2446209/1‑148 (MQ=255)
               aaTACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGaaa                                                                                                                                   <  2:941434/148‑1 (MQ=255)
                aTACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAAc                                                                                                                                  <  2:3280754/148‑1 (MQ=255)
                   ccGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATGACAGGAACCgg                                                                                                                               >  2:1839152/1‑148 (MQ=255)
                    cGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGc                                                                                                                              <  1:2289124/148‑1 (MQ=255)
                     gATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCa                                                                                                                             >  2:1109296/1‑148 (MQ=255)
                      aTCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCAt                                                                                                                            >  1:640672/1‑148 (MQ=255)
                       tCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATc                                                                                                                           <  1:395454/148‑1 (MQ=255)
                       tCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATc                                                                                                                           >  1:3199815/1‑148 (MQ=255)
                           caacGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGca                                                                                                                       >  1:811982/1‑148 (MQ=255)
                            aacGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGcac                                                                                                                      >  1:1462967/1‑148 (MQ=255)
                              cGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAg                                                                                                                    <  2:1545148/148‑1 (MQ=255)
                              cGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAg                                                                                                                    <  2:2080191/148‑1 (MQ=255)
                               gCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGTACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGc                                                                                                                   <  2:813765/148‑1 (MQ=255)
                                 ccAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGccc                                                                                                                 >  2:218426/1‑148 (MQ=255)
                                   agagCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTg                                                                                                               <  1:763338/148‑1 (MQ=255)
                                       cgcgGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAg                                                                                                           >  1:2513256/1‑148 (MQ=255)
                                         cgGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGt                                                                                                         <  2:1382286/148‑1 (MQ=255)
                                          ggAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTg                                                                                                        <  1:839662/148‑1 (MQ=255)
                                           gAAATGGGCAAGTGCTCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGGGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACATAAACCGGCATGAGCACAGACCTGGTAGGTGa                                                                                                       <  2:1129470/148‑1 (MQ=255)
                                           gAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGa                                                                                                       >  2:3075792/1‑148 (MQ=255)
                                             aaaGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACa                                                                                                     >  1:703533/1‑148 (MQ=255)
                                             aaaGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCAGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACa                                                                                                     >  1:2886017/1‑148 (MQ=255)
                                               aGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTACCGGTACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGt                                                                                                   <  2:3265769/148‑1 (MQ=255)
                                                gggCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGtt                                                                                                  >  2:1499315/1‑148 (MQ=255)
                                                     agtGCGCGGTGGGCTGCTCGTGCTGCGGAACCGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGt                                                                                             <  1:2186649/146‑1 (MQ=255)
                                                     aTTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGt                                                                                             <  2:863536/148‑1 (MQ=255)
                                                       tGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCAGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCg                                                                                           <  2:1749519/148‑1 (MQ=255)
                                                           cgGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTcc                                                                                       <  1:3345614/148‑1 (MQ=255)
                                                               gggCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCa                                                                                   >  1:3251403/1‑148 (MQ=255)
                                                               gggCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCAGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCa                                                                                   <  1:3302709/148‑1 (MQ=255)
                                                                 gctgctCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACt                                                                                 >  1:2716009/1‑148 (MQ=255)
                                                                   tgctCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCg                                                                               <  1:104647/148‑1 (MQ=255)
                                                                     ctctTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGcc                                                                             <  1:3426658/148‑1 (MQ=255)
                                                                     ctctTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGcc                                                                             <  2:1216743/148‑1 (MQ=255)
                                                                          gctgcGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGtt                                                                        <  2:1484631/148‑1 (MQ=255)
                                                                            tgcGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGa                                                                      <  1:2102670/148‑1 (MQ=255)
                                                                               ggAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATgg                                                                   <  2:2046772/148‑1 (MQ=255)
                                                                               ggAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATgg                                                                   <  1:744915/148‑1 (MQ=255)
                                                                                  aCAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGa                                                                <  1:1150841/148‑1 (MQ=255)
                                                                                      acacAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACggg                                                            >  2:3112285/1‑148 (MQ=255)
                                                                                         cacgCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAg                                                         >  2:1713365/4‑148 (MQ=255)
                                                                                          aGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGc                                                        >  2:3451835/1‑148 (MQ=255)
                                                                                           ggCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCg                                                       >  2:2850083/1‑148 (MQ=255)
                                                                                           ggCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCg                                                       <  1:2181066/148‑1 (MQ=255)
                                                                                                tgGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGc                                                  >  2:3052508/1‑148 (MQ=255)
                                                                                                   cgcTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGa                                               <  2:2545220/148‑1 (MQ=255)
                                                                                                    gcTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGAt                                              <  1:469463/148‑1 (MQ=255)
                                                                                                    gcTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTGTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGAt                                              <  2:3351645/148‑1 (MQ=255)
                                                                                                     cTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCAGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATc                                             >  2:726293/1‑148 (MQ=255)
                                                                                                      tGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATcg                                            <  1:1195829/148‑1 (MQ=255)
                                                                                                      tGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATcg                                            <  1:1262491/148‑1 (MQ=255)
                                                                                                      tGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATcg                                            <  2:2219658/148‑1 (MQ=255)
                                                                                                        cGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATcgcg                                          <  1:2768297/148‑1 (MQ=255)
                                                                                                        cGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATcgcg                                          <  1:674901/148‑1 (MQ=255)
                                                                                                        cGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATcgcg                                          <  2:1078912/148‑1 (MQ=255)
                                                                                                           tgttgGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTa                                       >  2:491680/1‑148 (MQ=255)
                                                                                                             ttgGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTaaa                                     >  1:3160585/1‑148 (MQ=255)
                                                                                                              tgGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAAt                                    <  1:3155721/148‑1 (MQ=255)
                                                                                                              tgGAGACCACTTCGCTTAACGTAACGTTTTTGGTTTCTGGCGTATAACAGAAACCGGCATCAGCACAGCCCTGGTATGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTatat                                    <  2:1218772/148‑1 (MQ=255)
                                                                                                                  gACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCg                                <  1:2351067/148‑1 (MQ=255)
                                                                                                                    ccACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCt                              <  1:618197/148‑1 (MQ=255)
                                                                                                                    ccACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCt                              <  2:2161912/148‑1 (MQ=255)
                                                                                                                    ccACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCt                              >  1:1617546/1‑148 (MQ=255)
                                                                                                                     cACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCtt                             <  1:1392819/148‑1 (MQ=255)
                                                                                                                     cACTTCGATTAACGGAACGGTTTTGGATTCTGGCGGATAAAAGAAAACGGCATCAGCACAGCCCTGGTAGGTGACAGTAAACGTCGCTCCCGCACTAGCCGGGTTGATGGTGACTGGAAGCTTCAGACGATCGCGGTAAAtctatctt                             >  2:2421786/1‑143 (MQ=255)
                                                                                                                        ttCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTg                          >  1:1658509/1‑148 (MQ=255)
                                                                                                                        ttCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTg                          <  1:2151989/148‑1 (MQ=255)
                                                                                                                         tCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGc                         >  2:2878852/1‑148 (MQ=255)
                                                                                                                              tAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTaa                    >  2:3413970/1‑148 (MQ=255)
                                                                                                                                   gaacggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTc               <  2:1652703/148‑1 (MQ=255)
                                                                                                                                     acggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCAt             <  1:2895530/148‑1 (MQ=255)
                                                                                                                                     acggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCAt             <  1:949958/148‑1 (MQ=255)
                                                                                                                                     acggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCAt             <  1:2489805/148‑1 (MQ=255)
                                                                                                                                        gTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCtt          >  1:2440964/1‑148 (MQ=255)
                                                                                                                                        gTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCtt          <  2:640672/148‑1 (MQ=255)
                                                                                                                                         ttttGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTc         <  1:1284269/148‑1 (MQ=255)
                                                                                                                                         ttttGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTc         <  2:2621162/148‑1 (MQ=255)
                                                                                                                                         ttttGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTc         >  1:377855/1‑148 (MQ=255)
                                                                                                                                         ttttGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTc         <  1:748486/148‑1 (MQ=255)
                                                                                                                                          tttGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCa        <  2:3477560/148‑1 (MQ=255)
                                                                                                                                              gTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGcc    <  2:1462967/148‑1 (MQ=255)
                                                                                                                                              gTTTCAGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGcc    <  1:1700399/148‑1 (MQ=255)
                                                                                                                                                ttCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAg  <  1:3075792/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGTAAAGGCGATACCAATACCGATCAACAACGCCCAGAGCGCGGAAAAGGGCAATTGCGCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAG  >  CP009273/4354320‑4354611

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: ‑