Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 232,860 G→T 92.3% 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,8600GT92.3% 923.0 / 60.5 285D103Y (GAT→TAT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base G (13/9);  new base T (130/133);  total (143/142)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 285 total aligned reads displayed.

CCCGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGC  >  CP009273/232713‑233006
                                                                                                                                                   |                                                                                                                                                  
cccGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTat                                                                                                                                                  >  2:2550974/1‑150 (MQ=255)
cccGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTat                                                                                                                                                  <  2:1351816/150‑1 (MQ=255)
 ccGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGata                                                                                                                                                 <  1:1254006/150‑1 (MQ=255)
  cGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTatat                                                                                                                                                <  1:2540147/150‑1 (MQ=255)
  cGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTatat                                                                                                                                                <  1:1325011/150‑1 (MQ=255)
  cGATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTatat                                                                                                                                                >  2:2557269/1‑150 (MQ=255)
   gATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATc                                                                                                                                               >  2:2795434/1‑150 (MQ=255)
   gATCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATc                                                                                                                                               <  2:4310027/150‑1 (MQ=255)
    aTCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCg                                                                                                                                              >  1:793012/1‑150 (MQ=255)
     tCGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCgg                                                                                                                                             >  1:306416/1‑150 (MQ=255)
      cGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGc                                                                                                                                            >  1:3226227/1‑150 (MQ=255)
      cGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGc                                                                                                                                            <  2:4348095/150‑1 (MQ=255)
      cGGCTGGTGGATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGc                                                                                                                                            <  1:3036070/150‑1 (MQ=255)
            gtggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATg                                                                                                                                      >  1:633195/1‑150 (MQ=255)
            gtggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATg                                                                                                                                      >  2:1244275/1‑150 (MQ=255)
            gtggATCCGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATg                                                                                                                                      >  1:1542208/1‑150 (MQ=255)
                  ccGGAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTAc                                                                                                                                <  1:3178327/150‑1 (MQ=255)
                   cGGAAGCCTTTGGCGTACATGGCATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACg                                                                                                                               >  2:2765535/1‑150 (MQ=255)
                    ggAAGCCTTTGTCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGa                                                                                                                              <  1:1721410/150‑1 (MQ=255)
                     gAAGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAg                                                                                                                             <  2:2200700/150‑1 (MQ=255)
                       aGCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGtt                                                                                                                           >  2:1528518/1‑150 (MQ=255)
                        gCCTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGttt                                                                                                                          <  1:3820870/150‑1 (MQ=255)
                         ccTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGTCTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGtttt                                                                                                                         <  1:758504/150‑1 (MQ=255)
                         ccTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGtttt                                                                                                                         >  2:886033/1‑150 (MQ=255)
                          cTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTc                                                                                                                        >  1:397693/1‑150 (MQ=255)
                          cTTTGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTc                                                                                                                        >  2:308937/1‑150 (MQ=255)
                           tttGGCGTACATGGTATTGCCGATGAATTTGTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCg                                                                                                                       <  2:2390225/150‑1 (MQ=255)
                             tGGCGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGtt                                                                                                                     <  1:2687607/150‑1 (MQ=255)
                                cGTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCt                                                                                                                  <  2:2329348/150‑1 (MQ=255)
                                 gTACATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCtt                                                                                                                 >  1:911946/1‑150 (MQ=255)
                                   aCATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTaa                                                                                                               <  2:4276017/150‑1 (MQ=255)
                                    cATGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAg                                                                                                              >  2:3072218/1‑150 (MQ=255)
                                     aTGGTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgc                                                                                                             <  1:4375719/150‑1 (MQ=255)
                                       ggTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgcgc                                                                                                           >  2:3276776/1‑150 (MQ=255)
                                        gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGcgcg                                                                                                          <  2:536768/150‑1 (MQ=255)
                                             gCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATAtt                                                                                                     <  1:3049734/150‑1 (MQ=255)
                                                gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCg                                                                                                  <  2:3508585/150‑1 (MQ=255)
                                                gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCg                                                                                                  >  1:1842404/1‑150 (MQ=255)
                                                gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCg                                                                                                  <  2:1118057/150‑1 (MQ=255)
                                                   gAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAg                                                                                               >  2:3718684/1‑150 (MQ=255)
                                                    aaTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa                                                                                              >  1:4473212/1‑150 (MQ=255)
                                                    aaTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa                                                                                              <  1:4532750/150‑1 (MQ=255)
                                                    aaTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa                                                                                              <  1:3853362/150‑1 (MQ=255)
                                                    aaTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa                                                                                              >  1:4231858/1‑150 (MQ=255)
                                                     aTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATACCGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAgcc                                                                                             <  1:1958603/150‑3 (MQ=255)
                                                       ttttGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCa                                                                                           <  2:305240/150‑1 (MQ=255)
                                                       ttttGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCa                                                                                           <  1:2490971/150‑1 (MQ=255)
                                                          tGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATa                                                                                        <  2:747722/150‑1 (MQ=255)
                                                          tGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATa                                                                                        <  2:1542999/150‑1 (MQ=255)
                                                           gCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATAc                                                                                       <  2:3348899/150‑1 (MQ=255)
                                                            cTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACt                                                                                      <  1:116737/150‑1 (MQ=255)
                                                            cTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACt                                                                                      >  1:4066824/1‑150 (MQ=255)
                                                              cGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACttt                                                                                    <  2:2541408/150‑1 (MQ=255)
                                                                       cACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGt                                                                           >  2:95778/1‑150 (MQ=255)
                                                                       cACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCAGTAAGGt                                                                           >  1:4170650/1‑150 (MQ=255)
                                                                          gTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCAc                                                                        >  1:3532123/1‑150 (MQ=255)
                                                                          gTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCAc                                                                        >  1:219888/1‑150 (MQ=255)
                                                                             tGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGa                                                                     <  2:4035049/150‑1 (MQ=255)
                                                                             tGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGa                                                                     >  1:2059967/1‑150 (MQ=255)
                                                                                cGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAg                                                                  >  2:2597801/1‑150 (MQ=255)
                                                                                  aaGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGcc                                                                >  1:712426/1‑150 (MQ=255)
                                                                                    gTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCtt                                                              <  2:4646221/150‑1 (MQ=255)
                                                                                     tAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTg                                                             >  2:2870391/1‑150 (MQ=255)
                                                                                          gATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTg                                                        <  1:958608/150‑1 (MQ=255)
                                                                                            tGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGc                                                      <  2:3201664/150‑1 (MQ=255)
                                                                                              aGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGa                                                    >  1:4725589/1‑150 (MQ=255)
                                                                                                ttCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAgg                                                  <  1:2793844/150‑1 (MQ=255)
                                                                                                 tCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGa                                                 >  1:3311029/1‑150 (MQ=255)
                                                                                                   aTGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGaaa                                               >  2:4738264/1‑150 (MQ=255)
                                                                                                    tGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGaaca                                              >  2:2681921/1‑148 (MQ=255)
                                                                                                      gACTATAGTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATg                                            <  1:236461/150‑1 (MQ=255)
                                                                                                        cTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGtt                                          >  2:1177134/1‑150 (MQ=255)
                                                                                                        cTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGtt                                          >  1:3363/1‑150 (MQ=255)
                                                                                                          atcttCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTc                                        <  2:4242060/147‑1 (MQ=255)
                                                                                                          atatTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTc                                        >  1:1295233/1‑150 (MQ=255)
                                                                                                            atTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTccc                                      >  1:2154630/1‑150 (MQ=255)
                                                                                                            atTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTccc                                      >  2:1913423/1‑150 (MQ=255)
                                                                                                                gcggcgcggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTa                                  <  2:3299449/150‑1 (MQ=255)
                                                                                                                gcggcgcggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTa                                  <  2:2021089/150‑1 (MQ=255)
                                                                                                                    cgcggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAgcg                              >  2:451677/1‑150 (MQ=255)
                                                                                                                      cggAGTTGGTTATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCa                            >  1:1083129/1‑150 (MQ=255)
                                                                                                                      cggAGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCa                            >  1:2971187/1‑150 (MQ=255)
                                                                                                                         aGTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACa                         >  2:3918162/1‑150 (MQ=255)
                                                                                                                          gTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAATACCAATACTTTCTGTAAGGTCACAGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAg                        >  2:293498/1‑150 (MQ=255)
                                                                                                                          gTTGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAg                        <  2:2481729/150‑1 (MQ=255)
                                                                                                                            tGGTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGcc                      <  2:633195/150‑1 (MQ=255)
                                                                                                                              gTGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTc                    <  2:2090013/150‑1 (MQ=255)
                                                                                                                               tGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCg                   <  1:2577641/150‑1 (MQ=255)
                                                                                                                               tGATCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCg                   >  1:4580035/1‑150 (MQ=255)
                                                                                                                                 aTCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGAt                 <  2:793012/150‑1 (MQ=255)
                                                                                                                                 aTCCATAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGAt                 <  1:1032280/150‑1 (MQ=255)
                                                                                                                                      tAACGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGtt            >  1:4416299/1‑150 (MQ=255)
                                                                                                                                       aaCGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTa           <  1:3726116/150‑1 (MQ=255)
                                                                                                                                       aaCGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTa           <  2:4725589/150‑1 (MQ=255)
                                                                                                                                        aCGCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAt          <  1:1528518/150‑1 (MQ=255)
                                                                                                                                          gCAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtgt        <  2:3599240/150‑1 (MQ=255)
                                                                                                                                           cAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtgtg       <  1:2765535/150‑1 (MQ=255)
                                                                                                                                           cAGCGTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTAtgtg       >  1:3516715/1‑150 (MQ=255)
                                                                                                                                               gTTCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCg   >  1:1429249/1‑150 (MQ=255)
                                                                                                                                                ttCTATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTTTCTGTAAGGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGc  >  2:1618142/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: