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,8600GT99.4% 1011.3 / ‑2.4 308D103Y (GAT→TAT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (127/179);  total (129/179)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 312 total aligned reads displayed.

CCCGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCG  >  CP009273/232713‑233005
                                                                                                                                                   |                                                                                                                                                 
cccGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTat                                                                                                                                                 <  2:1647634/150‑1 (MQ=255)
cccGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTat                                                                                                                                                 <  1:179924/150‑1 (MQ=255)
cccGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTat                                                                                                                                                 <  2:1912041/150‑1 (MQ=255)
 ccGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCtata                                                                                                                                                >  2:3414974/1‑150 (MQ=255)
  cGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTatat                                                                                                                                               <  2:1540642/150‑1 (MQ=255)
  cGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTatat                                                                                                                                               >  2:1811498/1‑150 (MQ=255)
   gATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATc                                                                                                                                              >  1:4401814/1‑150 (MQ=255)
   gATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATc                                                                                                                                              <  1:3395806/150‑1 (MQ=255)
     tCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCgg                                                                                                                                            <  1:2648815/150‑1 (MQ=255)
       ggCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCt                                                                                                                                          >  1:223452/1‑150 (MQ=255)
            gtggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATg                                                                                                                                     <  1:4562087/150‑1 (MQ=255)
             tggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATgg                                                                                                                                    <  2:1048805/150‑1 (MQ=255)
                aTCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACt                                                                                                                                 <  1:184732/150‑1 (MQ=255)
                   cGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACg                                                                                                                              >  1:2013046/1‑150 (MQ=255)
                   cGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACg                                                                                                                              <  1:3121699/150‑1 (MQ=255)
                    ggAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGa                                                                                                                             <  1:4285116/150‑1 (MQ=255)
                    ggAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGa                                                                                                                             >  1:3070684/1‑150 (MQ=255)
                        gCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGttt                                                                                                                         <  1:2432303/150‑1 (MQ=255)
                        gCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGttt                                                                                                                         >  1:1210190/1‑150 (MQ=255)
                        gCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGttt                                                                                                                         >  2:4004874/1‑150 (MQ=255)
                          cTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTc                                                                                                                       >  2:740621/1‑150 (MQ=255)
                          cTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTc                                                                                                                       >  2:68863/1‑150 (MQ=255)
                             tGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGATTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGtt                                                                                                                    <  1:2972410/150‑1 (MQ=255)
                              ggCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTg                                                                                                                   >  2:2989545/1‑150 (MQ=255)
                              ggCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGAGGAGGTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTg                                                                                                                   >  2:2339390/1‑150 (MQ=255)
                                cGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCt                                                                                                                 >  1:1538639/1‑150 (MQ=255)
                                 gTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCtt                                                                                                                <  1:3220843/150‑1 (MQ=255)
                                  tACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTa                                                                                                               >  2:1253411/1‑150 (MQ=255)
                                   aCATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTaa                                                                                                              <  2:1501481/150‑1 (MQ=255)
                                   aCATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTaa                                                                                                              >  1:3901283/1‑150 (MQ=255)
                                    cATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAg                                                                                                             <  1:2999790/150‑1 (MQ=255)
                                    cATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAg                                                                                                             <  2:5183859/150‑1 (MQ=255)
                                     aTGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgc                                                                                                            <  1:2864199/150‑1 (MQ=255)
                                     aTGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgc                                                                                                            <  2:3904070/150‑1 (MQ=255)
                                        gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGcgcg                                                                                                         >  1:3781520/1‑150 (MQ=255)
                                         tATTGCCGATGAATTTTTGCTCGATAAGCCCCCGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGa                                                                                                        <  2:2771898/150‑1 (MQ=255)
                                            tGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGatat                                                                                                     <  2:4749000/150‑1 (MQ=255)
                                              ccGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTc                                                                                                   <  1:3881264/150‑1 (MQ=255)
                                               cGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTcc                                                                                                  <  2:481825/150‑1 (MQ=255)
                                                gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCg                                                                                                 <  2:2465228/150‑1 (MQ=255)
                                                gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCg                                                                                                 <  2:4094426/150‑1 (MQ=255)
                                                   gAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAg                                                                                              <  1:2149003/150‑1 (MQ=255)
                                                         ttGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCTATATTCCGAAGACCAAt                                                                                        <  2:2850065/150‑1 (MQ=255)
                                                         ttGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAAt                                                                                        <  2:311116/150‑1 (MQ=255)
                                                           gCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGTCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATAc                                                                                      >  1:1162671/1‑150 (MQ=255)
                                                                ttAAGCCGTCGTTTGCCGAAGTAGCCGATGTGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTTCGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCt                                                                                 <  1:953970/149‑1 (MQ=255)
                                                                aTAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCt                                                                                 <  1:3384261/150‑1 (MQ=255)
                                                                aTAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCt                                                                                 <  2:4902693/150‑1 (MQ=255)
                                                                  aaGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGt                                                                               <  1:5027480/150‑1 (MQ=255)
                                                                  aaGCCCACGTTTGCCGAAGGAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGt                                                                               <  2:1389621/150‑1 (MQ=255)
                                                                    gCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTaa                                                                             <  1:740621/150‑1 (MQ=255)
                                                                    gCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTaa                                                                             <  1:4445373/150‑1 (MQ=255)
                                                                     cccACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAg                                                                            >  1:31562/1‑150 (MQ=255)
                                                                      ccACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAgg                                                                           <  2:1145938/150‑1 (MQ=255)
                                                                       cACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGt                                                                          >  1:2097130/1‑150 (MQ=255)
                                                                        aCGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTc                                                                         <  2:3851411/150‑1 (MQ=255)
                                                                          gTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCAc                                                                       >  2:1783003/1‑150 (MQ=255)
                                                                          gTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCAc                                                                       >  2:4439797/1‑150 (MQ=255)
                                                                          gTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCAc                                                                       <  1:269125/150‑1 (MQ=255)
                                                                            ttGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCg                                                                     <  2:2226696/150‑1 (MQ=255)
                                                                             tGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGa                                                                    <  2:3166090/150‑1 (MQ=255)
                                                                                cGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAg                                                                 >  1:157549/1‑150 (MQ=255)
                                                                                 gAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGc                                                                <  1:5138054/150‑1 (MQ=255)
                                                                                   aGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCt                                                              <  1:3488280/150‑1 (MQ=255)
                                                                                         cGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGt                                                        <  2:5325294/150‑1 (MQ=255)
                                                                                          gATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTg                                                       <  1:816245/150‑1 (MQ=255)
                                                                                             gAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCg                                                    <  1:1015261/150‑1 (MQ=255)
                                                                                             gAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCg                                                    <  1:4229104/150‑1 (MQ=255)
                                                                                             gAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCg                                                    <  2:859730/150‑1 (MQ=255)
                                                                                              aGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGa                                                   >  2:908245/1‑150 (MQ=255)
                                                                                                  cATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGaa                                               >  2:3600654/1‑150 (MQ=255)
                                                                                                     ggACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAAt                                            <  2:413080/150‑1 (MQ=255)
                                                                                                     ggACTATAGTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAAt                                            <  2:4221048/150‑1 (MQ=255)
                                                                                                       aCTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGt                                          >  2:805982/1‑150 (MQ=255)
                                                                                                        cTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGtt                                         <  2:2449277/150‑1 (MQ=255)
                                                                                                         tataTTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGttt                                        <  2:3894049/150‑1 (MQ=255)
                                                                                                          atatTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTc                                       >  1:3987293/1‑150 (MQ=255)
                                                                                                          atatTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTc                                       >  2:3744309/1‑150 (MQ=255)
                                                                                                            atTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTccc                                     >  1:484965/1‑150 (MQ=255)
                                                                                                             ttCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCg                                    <  2:1683406/150‑1 (MQ=255)
                                                                                                             ttCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCg                                    >  1:362439/1‑150 (MQ=255)
                                                                                                             ttCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCg                                    <  1:2326278/150‑1 (MQ=255)
                                                                                                                  ggcgcggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAg                               >  1:3252765/1‑150 (MQ=255)
                                                                                                                    cgcggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAgcg                             <  2:659958/150‑1 (MQ=255)
                                                                                                                        gAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAAc                         <  2:3958246/150‑1 (MQ=255)
                                                                                                                          gTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAg                       >  2:4323928/1‑150 (MQ=255)
                                                                                                                          gTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAg                       >  2:229279/1‑150 (MQ=255)
                                                                                                                             ggTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAATCGCAACAGCCt                    <  2:4655181/150‑1 (MQ=255)
                                                                                                                             ggTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCt                    >  2:3344186/1‑150 (MQ=255)
                                                                                                                             ggTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCt                    >  1:1729783/1‑150 (MQ=255)
                                                                                                                             ggTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCt                    >  1:1725944/1‑150 (MQ=255)
                                                                                                                              gTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTc                   <  2:2697828/150‑1 (MQ=255)
                                                                                                                              gTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCACAACAGCCTc                   <  1:4517810/150‑1 (MQ=255)
                                                                                                                               tGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCg                  <  1:526053/150‑1 (MQ=255)
                                                                                                                                gATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGa                 >  2:4389966/1‑150 (MQ=255)
                                                                                                                                  tCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATg               >  2:4671120/1‑150 (MQ=255)
                                                                                                                                       aaCGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTa          >  2:2924601/1‑150 (MQ=255)
                                                                                                                                         cGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtg        <  1:2580131/150‑1 (MQ=255)
                                                                                                                                          gCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtgt       <  1:2347841/150‑1 (MQ=255)
                                                                                                                                             gCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCt    <  2:2146530/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCCGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCG  >  CP009273/232713‑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: ‑