Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 232,860 G→T 100% 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,8600GT97.9% 1118.5 / 12.5 327D103Y (GAT→TAT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base T (152/168);  total (154/173)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 332 total aligned reads displayed.

CCCGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGC  >  CP009273/232713‑233006
                                                                                                                                                   |                                                                                                                                                  
cccGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTat                                                                                                                                                  <  2:1171405/150‑1 (MQ=255)
 ccGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCtata                                                                                                                                                 <  2:263887/150‑1 (MQ=255)
  cGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTatatat                                                                                                                                                >  2:1118384/1‑150 (MQ=255)
  cGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTatat                                                                                                                                                <  2:2066750/150‑1 (MQ=255)
    aTCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCg                                                                                                                                              >  2:94700/1‑150 (MQ=255)
    aTCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCg                                                                                                                                              >  2:2364331/1‑150 (MQ=255)
         cTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCttt                                                                                                                                         <  2:457153/150‑1 (MQ=255)
         cTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCttt                                                                                                                                         <  2:3382749/150‑1 (MQ=255)
             tggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATgg                                                                                                                                     <  2:3287571/150‑1 (MQ=255)
               gATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCTGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGAc                                                                                                                                   <  2:3874832/150‑1 (MQ=255)
                 tCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTa                                                                                                                                 <  2:2029480/150‑1 (MQ=255)
                 tCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTa                                                                                                                                 <  2:3709582/150‑1 (MQ=255)
                  ccGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTAc                                                                                                                                >  1:4012446/1‑150 (MQ=255)
                     gAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAg                                                                                                                             <  1:842231/150‑1 (MQ=255)
                     gAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAg                                                                                                                             >  1:2420233/1‑150 (MQ=255)
                     gAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAg                                                                                                                             <  1:1587009/150‑1 (MQ=255)
                       aGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGtt                                                                                                                           <  2:244540/150‑1 (MQ=255)
                       aGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGtt                                                                                                                           >  1:4219360/1‑150 (MQ=255)
                            ttGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGt                                                                                                                      <  2:1195404/150‑1 (MQ=255)
                             tGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGtt                                                                                                                     <  1:2063646/150‑1 (MQ=255)
                             tGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGtt                                                                                                                     >  2:1308231/1‑150 (MQ=255)
                                cGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCt                                                                                                                  <  2:677695/150‑1 (MQ=255)
                                cGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCt                                                                                                                  <  1:476899/150‑1 (MQ=255)
                                 gTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCtt                                                                                                                 >  2:151596/1‑150 (MQ=255)
                                         tATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGa                                                                                                         >  1:3894557/1‑150 (MQ=255)
                                         tATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGa                                                                                                         >  1:4127781/1‑150 (MQ=255)
                                               cGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTcc                                                                                                   <  2:2346688/150‑1 (MQ=255)
                                               cGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTcc                                                                                                   <  1:3133051/150‑1 (MQ=255)
                                               cGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTcc                                                                                                   <  1:3108003/150‑1 (MQ=255)
                                                  tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGaa                                                                                                >  2:4247400/1‑150 (MQ=255)
                                                  tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGaa                                                                                                <  1:4099092/150‑1 (MQ=255)
                                                   gAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAg                                                                                               <  1:2480480/150‑1 (MQ=255)
                                                     aTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGAc                                                                                             >  1:3517635/1‑150 (MQ=255)
                                                          tGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATa                                                                                        <  1:4351887/150‑1 (MQ=255)
                                                          tGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATa                                                                                        <  1:1004726/150‑1 (MQ=255)
                                                            cTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACt                                                                                      <  1:2913251/150‑1 (MQ=255)
                                                               gATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTc                                                                                   <  1:297083/150‑1 (MQ=255)
                                                                aTAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCt                                                                                  >  1:2197771/1‑150 (MQ=255)
                                                                     cccACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAg                                                                             >  2:700560/1‑150 (MQ=255)
                                                                      ccACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAgg                                                                            >  1:1077160/1‑150 (MQ=255)
                                                                      ccACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAgg                                                                            >  1:1129852/1‑150 (MQ=255)
                                                                        aCGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTc                                                                          >  1:3747937/1‑150 (MQ=255)
                                                                        aCGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTc                                                                          <  2:1567139/150‑1 (MQ=255)
                                                                           tttGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCAcc                                                                       <  2:3753583/150‑1 (MQ=255)
                                                                            ttGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCg                                                                      <  1:3173465/150‑1 (MQ=255)
                                                                             tGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGa                                                                     >  1:1820809/1‑150 (MQ=255)
                                                                              gCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGAt                                                                    <  1:2097490/150‑1 (MQ=255)
                                                                                 gAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGc                                                                 <  2:4488410/150‑1 (MQ=255)
                                                                                  aaGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGcc                                                                <  1:3276764/150‑1 (MQ=255)
                                                                                   aGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCt                                                               <  2:296117/150‑1 (MQ=255)
                                                                                   aGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCt                                                               >  2:2747861/1‑150 (MQ=255)
                                                                                    gTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCtt                                                              >  1:1628523/1‑150 (MQ=255)
                                                                                    gTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCtt                                                              <  2:1659892/150‑1 (MQ=255)
                                                                                      aGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACGTTCTGTAAGGTCACCGATAGCCTTGc                                                            >  1:2601894/1‑150 (MQ=255)
                                                                                       gCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCg                                                           >  2:4086849/1‑150 (MQ=255)
                                                                                         cGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGt                                                         <  1:3203824/150‑1 (MQ=255)
                                                                                          gATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTg                                                        >  1:3600684/1‑150 (MQ=255)
                                                                                          gATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTg                                                        >  1:3022568/1‑150 (MQ=255)
                                                                                            tGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGc                                                      <  1:1778438/150‑1 (MQ=255)
                                                                                             gAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCg                                                     >  1:1379392/1‑150 (MQ=255)
                                                                                              aGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGa                                                    <  1:4169134/150‑1 (MQ=255)
                                                                                              aGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGa                                                    <  2:396124/150‑1 (MQ=255)
                                                                                              aGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGa                                                    <  1:883033/150‑1 (MQ=255)
                                                                                                    tGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGaaaa                                              >  1:1164780/1‑150 (MQ=255)
                                                                                                     ggACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAAt                                             >  1:1984560/1‑150 (MQ=255)
                                                                                                     ggACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAAt                                             <  1:3341116/150‑1 (MQ=255)
                                                                                                         tataTTCGCGGCGCGGAGTTGGTTATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGttt                                         >  1:1649794/1‑150 (MQ=255)
                                                                                                         tataTTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGttt                                         <  2:2101603/150‑1 (MQ=255)
                                                                                                          atatTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTc                                        <  1:897195/150‑1 (MQ=255)
                                                                                                            atTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTccc                                      >  2:863164/1‑150 (MQ=255)
                                                                                                            atTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTccc                                      >  1:528090/1‑150 (MQ=255)
                                                                                                            atTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTccc                                      <  1:4416945/150‑1 (MQ=255)
                                                                                                             ttCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCg                                     <  2:1776776/150‑1 (MQ=255)
                                                                                                             ttCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCg                                     <  2:515748/150‑1 (MQ=255)
                                                                                                               cgcggcgcggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGt                                   >  1:1282568/1‑150 (MQ=255)
                                                                                                               cgcggcgcggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGt                                   <  1:3463770/150‑1 (MQ=255)
                                                                                                                    cgcggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAgcg                              >  2:2270499/1‑150 (MQ=255)
                                                                                                                      cggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCa                            <  1:2011258/150‑1 (MQ=255)
                                                                                                                       ggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGTAAAATGTTTCCCGGTAAGAGCaa                           >  2:256092/1‑150 (MQ=255)
                                                                                                                       ggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCaa                           >  1:1579669/1‑150 (MQ=255)
                                                                                                                        gAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAAc                          <  2:1379991/150‑1 (MQ=255)
                                                                                                                         aGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACa                         >  1:779069/1‑150 (MQ=255)
                                                                                                                          gTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAg                        >  2:2808782/1‑150 (MQ=255)
                                                                                                                            tGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGcc                      >  2:2844596/1‑150 (MQ=255)
                                                                                                                              gTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTc                    >  2:3652578/1‑150 (MQ=255)
                                                                                                                                ggtCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGa                  >  1:3780195/3‑150 (MQ=255)
                                                                                                                                gATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGa                  <  2:3189809/150‑1 (MQ=255)
                                                                                                                                   ccATAACGCAGCGTTCTTTATCGGCTTTTTGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCATTACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGc               <  2:4424568/150‑1 (MQ=255)
                                                                                                                                    cATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCg              <  1:3396930/150‑1 (MQ=255)
                                                                                                                                        aCGCAGCGTTCTATATAGGCTTTATGGAATACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCTATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAt          <  2:2155485/150‑1 (MQ=255)
                                                                                                                                          gCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtgt        <  2:847891/150‑1 (MQ=255)
                                                                                                                                           cAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtgtg       <  1:958422/150‑1 (MQ=255)
                                                                                                                                            aGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGc      >  2:238197/1‑150 (MQ=255)
                                                                                                                                            aGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGc      >  1:1521431/1‑150 (MQ=255)
                                                                                                                                              cGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTc    <  1:3964527/150‑1 (MQ=255)
                                                                                                                                              cGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTc    <  2:2729852/150‑1 (MQ=255)
                                                                                                                                               gTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCg   <  1:383442/150‑1 (MQ=255)
                                                                                                                                               gTTCTATATAGGCGTTATGGACTACGAGTTTTCGTTTCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGGTGCGAGCAAAATGTTTCCCGGTAATCGCAACAGCCTCGATGCGTTATTGGCTCg   >  2:3533310/1‑150 (MQ=255)
                                                                                                                                                ttCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGc  >  1:2270292/1‑150 (MQ=255)
                                                                                                                                                ttCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCAATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGc  >  1:4587861/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCCGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGC  >  CP009273/232713‑233006

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: