Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 232,860 G→T 90.4% D103Y (GAT→TAT)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273232,8600GT90.4% 970.6 / 90.3 311D103Y (GAT→TAT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base G (16/14);  new base T (151/130);  total (167/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 312 total aligned reads displayed.

CCGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCG  >  CP009273/232714‑233005
                                                                                                                                                  |                                                                                                                                                 
ccGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTGTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTat                                                                                                                                                 <  2:3984589/149‑1 (MQ=255)
 cGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGACGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCtata                                                                                                                                                >  2:339348/1‑149 (MQ=255)
 cGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGata                                                                                                                                                >  2:1333176/1‑149 (MQ=255)
   aTCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATc                                                                                                                                              <  1:3086206/149‑1 (MQ=255)
   aTCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATc                                                                                                                                              >  2:4289307/1‑149 (MQ=255)
   aTCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATc                                                                                                                                              <  2:557875/149‑1 (MQ=255)
    tCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCg                                                                                                                                             >  2:2026851/1‑149 (MQ=255)
    tCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCg                                                                                                                                             >  2:3241212/1‑149 (MQ=255)
     cGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCgg                                                                                                                                            >  2:2228343/1‑149 (MQ=255)
      ggCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGc                                                                                                                                           >  1:4062653/1‑149 (MQ=255)
         tggtggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCttt                                                                                                                                        <  1:2783478/149‑1 (MQ=255)
         tggtggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCttt                                                                                                                                        >  1:2393383/1‑149 (MQ=255)
          ggtggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTa                                                                                                                                       <  1:3511324/149‑1 (MQ=255)
            tggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATg                                                                                                                                     >  1:1763704/1‑149 (MQ=255)
             ggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATgg                                                                                                                                    >  1:1863120/1‑149 (MQ=255)
              gATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGa                                                                                                                                   <  2:3045724/149‑1 (MQ=255)
                tCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACt                                                                                                                                 <  1:1015845/149‑1 (MQ=255)
                 ccGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTa                                                                                                                                <  1:4209909/149‑1 (MQ=255)
                  cGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTAc                                                                                                                               <  1:516646/149‑1 (MQ=255)
                    gAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGa                                                                                                                             >  1:3548762/1‑149 (MQ=255)
                     aaGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAg                                                                                                                            <  1:1725117/149‑1 (MQ=255)
                      aGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGt                                                                                                                           >  2:4746932/1‑149 (MQ=255)
                      aGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGt                                                                                                                           >  2:1164595/1‑149 (MQ=255)
                             ggCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGtt                                                                                                                    <  1:1175829/149‑1 (MQ=255)
                                gTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCt                                                                                                                 >  1:4499185/1‑149 (MQ=255)
                                 tACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCtt                                                                                                                <  2:1626046/149‑1 (MQ=255)
                                  aCATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTa                                                                                                               >  1:3248853/1‑149 (MQ=255)
                                  aCATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTa                                                                                                               >  2:2759776/1‑149 (MQ=255)
                                   cATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTaa                                                                                                              >  2:4988608/1‑149 (MQ=255)
                                    aTGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAg                                                                                                             <  1:4929113/149‑1 (MQ=255)
                                      ggTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgcg                                                                                                           >  1:4343558/1‑149 (MQ=255)
                                      ggTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgcg                                                                                                           >  2:4416453/1‑149 (MQ=255)
                                       gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGATTGCCGAAGTAGCCGATGAGTTCATGGACTATAGTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgcgc                                                                                                          <  2:2732765/149‑1 (MQ=255)
                                        tATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGcgcg                                                                                                         >  2:3775041/1‑149 (MQ=255)
                                         aTTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGa                                                                                                        <  2:2070677/149‑1 (MQ=255)
                                           tGCCGATGAATTTTTGCTCGATAAGCCCACGTGTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGata                                                                                                      <  1:3614081/149‑1 (MQ=255)
                                            gCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGatat                                                                                                     <  1:3951692/149‑1 (MQ=255)
                                              cGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTc                                                                                                   <  1:3707446/149‑1 (MQ=255)
                                                     tttttGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGAc                                                                                            <  2:699274/149‑1 (MQ=255)
                                                      ttttGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGAcc                                                                                           <  2:1058550/149‑1 (MQ=255)
                                                      ttttGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGAcc                                                                                           <  2:3677636/149‑1 (MQ=255)
                                                        ttGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCaa                                                                                         <  1:3289822/149‑1 (MQ=255)
                                                         tGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAAt                                                                                        <  1:679577/149‑1 (MQ=255)
                                                          gCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATa                                                                                       <  2:4863607/149‑1 (MQ=255)
                                                             cGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACtt                                                                                    <  1:1923774/149‑1 (MQ=255)
                                                              gATAAGCCCACGTGTGCCGAAGTAGCCGATGAGTTCATGGACTAGATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACttt                                                                                   <  1:1909807/149‑1 (MQ=255)
                                                               aTAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTc                                                                                  <  1:4993329/149‑1 (MQ=255)
                                                               aTAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTc                                                                                  <  2:4151825/149‑1 (MQ=255)
                                                                tAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCt                                                                                 <  1:733743/149‑1 (MQ=255)
                                                                  aGCCCATGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGt                                                                               >  2:3084760/1‑149 (MQ=255)
                                                                  aGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGt                                                                               >  2:3527346/1‑149 (MQ=255)
                                                                   gCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTa                                                                              <  2:5141253/149‑1 (MQ=255)
                                                                     ccACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAg                                                                            >  1:4276639/1‑149 (MQ=255)
                                                                       aCGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGt                                                                          <  2:4499050/149‑1 (MQ=255)
                                                                       aCGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGt                                                                          >  2:952821/1‑149 (MQ=255)
                                                                          tttGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCAc                                                                       <  1:4461905/149‑1 (MQ=255)
                                                                            tGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCg                                                                     >  2:2245264/1‑149 (MQ=255)
                                                                             gCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGa                                                                    <  1:3431290/149‑1 (MQ=255)
                                                                              ccGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGAt                                                                   <  1:1978233/149‑1 (MQ=255)
                                                                               cGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATACCGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATa                                                                  <  2:1896136/149‑1 (MQ=255)
                                                                                gAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAg                                                                 >  2:4594443/1‑149 (MQ=255)
                                                                                 aaGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGc                                                                >  2:3900173/1‑149 (MQ=255)
                                                                                   gTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCt                                                              >  2:1582277/1‑149 (MQ=255)
                                                                                         gATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATATCCTTGCGGt                                                        >  1:1318903/1‑149 (MQ=255)
                                                                                         gATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGt                                                        >  1:2322138/1‑149 (MQ=255)
                                                                                         gATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGt                                                        >  2:1237136/1‑149 (MQ=255)
                                                                                          aTGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTg                                                       >  2:4838751/1‑149 (MQ=255)
                                                                                           tGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTgg                                                      >  1:2161343/1‑149 (MQ=255)
                                                                                            gAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGc                                                     >  2:4382330/1‑149 (MQ=255)
                                                                                                 cATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGa                                                <  2:355845/149‑1 (MQ=255)
                                                                                                  aTGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGaa                                               <  2:1104271/149‑1 (MQ=255)
                                                                                                   tGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGaaa                                              <  2:1863120/149‑1 (MQ=255)
                                                                                                    ggACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGaaaa                                             <  1:1375375/149‑1 (MQ=255)
                                                                                                      aCTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATg                                           >  1:4817703/1‑149 (MQ=255)
                                                                                                          tatTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTc                                       >  1:2669856/1‑149 (MQ=255)
                                                                                                            ttCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTccc                                     >  2:1290453/1‑149 (MQ=255)
                                                                                                             tCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCg                                    <  1:3788150/149‑1 (MQ=255)
                                                                                                                  gcgtggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAACATGTTTCCCGGTAAg                               >  1:608425/1‑149 (MQ=255)
                                                                                                                       gAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCaa                          <  2:2561730/149‑1 (MQ=255)
                                                                                                                       gAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCaa                          >  2:3809875/1‑149 (MQ=255)
                                                                                                                        aGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAAc                         <  1:4450730/149‑1 (MQ=255)
                                                                                                                            ggTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGcc                     <  2:1000818/149‑1 (MQ=255)
                                                                                                                            ggTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGcc                     <  1:3366746/149‑1 (MQ=255)
                                                                                                                            ggTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGcc                     <  1:2119352/149‑1 (MQ=255)
                                                                                                                             gTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCt                    <  1:4675623/149‑1 (MQ=255)
                                                                                                                              tGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTc                   >  1:3339029/1‑149 (MQ=255)
                                                                                                                              tGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTc                   >  1:2807944/1‑149 (MQ=255)
                                                                                                                                aTCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGa                 <  2:4005536/149‑1 (MQ=255)
                                                                                                                                aTCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGa                 <  2:4972933/149‑1 (MQ=255)
                                                                                                                                 tCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGAt                <  1:3151827/149‑1 (MQ=255)
                                                                                                                                     tAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGt            <  2:141929/149‑1 (MQ=255)
                                                                                                                                       aCGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTa          <  2:613555/149‑1 (MQ=255)
                                                                                                                                       aCGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTa          <  2:4343558/149‑1 (MQ=255)
                                                                                                                                       aCGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTa          <  1:2853142/149‑1 (MQ=255)
                                                                                                                                        cGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAt         <  2:2577912/149‑1 (MQ=255)
                                                                                                                                         gCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTTCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtg        <  1:1142887/149‑1 (MQ=255)
                                                                                                                                         gCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtg        <  1:3043709/149‑1 (MQ=255)
                                                                                                                                         gCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtg        <  1:4694558/149‑1 (MQ=255)
                                                                                                                                         gCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtg        <  1:1052715/149‑1 (MQ=255)
                                                                                                                                           aGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtgtg      >  2:1776359/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCG  >  CP009273/232714‑233005

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: