Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,151,750 G→T 86.0% R178L (CGC→CTC)  holB → DNA polymerase III, delta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,151,7500GT86.0% 663.1 / 106.5 273R178L (CGC→CTC) holBDNA polymerase III, delta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (25/14);  new base T (127/104);  total (155/118)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 279 total aligned reads displayed.

TGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCGGGC  >  CP009273/1151603‑1151898
                                                                                                                                                   |                                                                                                                                                    
tGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTTActcg                                                                                                                                                    >  1:2891322/1‑150 (MQ=255)
   aaaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAg                                                                                                                                                 >  2:266127/1‑150 (MQ=255)
    aaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCAGGCTTTCACTCGAAGt                                                                                                                                                >  1:1828300/1‑150 (MQ=255)
     aCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTg                                                                                                                                               <  2:3253835/150‑1 (MQ=255)
      cGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGa                                                                                                                                              >  2:702465/1‑150 (MQ=255)
             agagCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTc                                                                                                                                       <  2:3476365/150‑1 (MQ=255)
               agCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACAGTACCTTGCGCCGCCGCCGGCACAGTACGCCGTGACCAGGCTTTCACTCGAAGTGACAAAGTcac                                                                                                                                     >  1:2323313/1‑150 (MQ=255)
                gCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTcaca                                                                                                                                    <  1:371322/150‑1 (MQ=255)
                 ccaccaGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTGTCACTCGAAGTGACAATGTCACAg                                                                                                                                   >  2:3019154/1‑150 (MQ=255)
                  caccaGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAgg                                                                                                                                  >  2:3627987/1‑150 (MQ=255)
                  caccaGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAgg                                                                                                                                  <  2:2427798/150‑1 (MQ=255)
                     cagcagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAATTGACAATGTCACAGGATg                                                                                                                               >  1:665139/1‑150 (MQ=255)
                        cagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAt                                                                                                                            >  1:2795203/1‑150 (MQ=255)
                          gAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTa                                                                                                                          >  1:300215/1‑150 (MQ=255)
                            aaCTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCCCAGGATGCATTACt                                                                                                                        >  2:4216452/1‑150 (MQ=255)
                            aaCTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACt                                                                                                                        >  1:679880/1‑150 (MQ=255)
                               ttGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGc                                                                                                                     <  1:4102436/150‑1 (MQ=255)
                                tGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGcc                                                                                                                    >  2:1163181/1‑150 (MQ=255)
                                tGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGAAGCATTACTTgac                                                                                                                    >  2:3432420/1‑148 (MQ=255)
                                 ggaatatCCTGGCTACCCGCGAGCCTGAACGTGTACTGGCAACATTACGTATAAGATTTCGGTTACAATAACATGCGCCGCCGCCGCAACAGTACGACGAGACCAGTCATTCACTCGAAGTTACAAATTCACAGCATGCAAAACTTGACg                                                                                                                   >  2:1491414/7‑150 (MQ=255)
                                      ttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTg                                                                                                              <  1:3787979/150‑1 (MQ=255)
                                      ttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTg                                                                                                              <  2:2348338/150‑1 (MQ=255)
                                         cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                                                           >  2:4418196/1‑150 (MQ=255)
                                          tGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCt                                                                                                          >  2:2896817/1‑150 (MQ=255)
                                          tGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCt                                                                                                          >  1:1743655/1‑150 (MQ=255)
                                             cTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTaa                                                                                                       >  2:913553/1‑150 (MQ=255)
                                               aCCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAgc                                                                                                     <  2:1731132/150‑1 (MQ=255)
                                                cccGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAgcg                                                                                                    >  2:1986138/1‑150 (MQ=255)
                                                    cgAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCgg                                                                                                >  2:468490/1‑150 (MQ=255)
                                                    cgAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCgg                                                                                                >  2:4364280/1‑150 (MQ=255)
                                                        ccTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCg                                                                                            >  1:622837/1‑150 (MQ=255)
                                                        ccTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCg                                                                                            >  2:4111613/1‑150 (MQ=255)
                                                           gAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCt                                                                                         <  2:1866011/150‑1 (MQ=255)
                                                            aaCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTg                                                                                        <  1:1963512/150‑1 (MQ=255)
                                                            aaCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTGGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTg                                                                                        <  1:1963552/150‑1 (MQ=255)
                                                                 ttACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGcgcg                                                                                   <  2:3806972/150‑1 (MQ=255)
                                                                    cTGGCAACATTACATAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTAAGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAATACATGCCGCATTGCGCTTAAGCGCCGGTTCGCCAGGCGCGGCa                                                                                >  1:3001529/1‑150 (MQ=255)
                                                                        cAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTgg                                                                            <  2:3796565/150‑1 (MQ=255)
                                                                        cAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTgg                                                                            <  1:3926586/150‑1 (MQ=255)
                                                                        cAACATTACGTAGTCGTTGTCGGTTAAATTACCATGCGCCGCCGCCGGAACAGTACCCCGTGACCAGGCTTTAACACGAAGTGACAAAGTCCCAGGATGCAATACATGCCGCATTGCGCTTAAGCGCCCGTTCGCCAGGCGCGGCACAgg                                                                            >  1:2306200/1‑150 (MQ=255)
                                                                         aaCATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGc                                                                           <  2:539278/150‑1 (MQ=255)
                                                                         aaCATTACGTAGTCGTTGTCGGTTACATTACCATTCGCCGCCGCCGGAACAGTACGCCGTGACCTTGCTTTCACGCGAAGTGACAAATTCACAGGATTCATTACTTGCCGCATTGCGCTTAATCGCCGGTTAGCCAGGCGCGGCACTGGc                                                                           >  1:518230/1‑150 (MQ=255)
                                                                          aCATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCg                                                                          <  1:1484678/150‑1 (MQ=255)
                                                                           cATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgt                                                                         <  1:4401152/150‑1 (MQ=255)
                                                                            aTTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgtt                                                                        <  1:3840435/150‑1 (MQ=255)
                                                                            aTTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgtt                                                                        >  2:4364584/1‑150 (MQ=255)
                                                                            aTTACGTAGTCGTTGTCGGTTACATTACCTTGAGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgtt                                                                        <  1:2595263/150‑1 (MQ=255)
                                                                             ttACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgttg                                                                       <  1:3766353/150‑1 (MQ=255)
                                                                              tACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgttgt                                                                      >  1:3297699/1‑150 (MQ=255)
                                                                               aCGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgttgtt                                                                     <  2:1728752/150‑1 (MQ=255)
                                                                               aCGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgttgtt                                                                     <  1:3627987/150‑1 (MQ=255)
                                                                               aCGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgttgtt                                                                     <  1:122574/150‑1 (MQ=255)
                                                                                cGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGttt                                                                    >  2:2999476/1‑150 (MQ=255)
                                                                                  tAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCa                                                                  >  2:1896858/1‑150 (MQ=255)
                                                                                   aGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAg                                                                 >  2:3725780/1‑150 (MQ=255)
                                                                                    gTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAgg                                                                <  1:3111786/150‑1 (MQ=255)
                                                                                     tCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTCAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAggg                                                               >  2:1409810/1‑150 (MQ=255)
                                                                                     tCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAggg                                                               <  2:2284581/150‑1 (MQ=255)
                                                                                     tCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAggg                                                               <  1:1896123/150‑1 (MQ=255)
                                                                                     tCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAggg                                                               <  2:1807384/150‑1 (MQ=255)
                                                                                       gTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGgag                                                             <  1:4047313/150‑1 (MQ=255)
                                                                                         tGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGAt                                                           >  2:2560216/1‑150 (MQ=255)
                                                                                         tGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGAt                                                           <  1:4365548/150‑1 (MQ=255)
                                                                                            cGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAAc                                                        <  1:1811961/150‑1 (MQ=255)
                                                                                               ttACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTgg                                                     <  1:74505/150‑1 (MQ=255)
                                                                                               ttACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTgg                                                     >  2:2219328/1‑150 (MQ=255)
                                                                                                tACATTACCTTGAGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGc                                                    >  2:2095864/1‑150 (MQ=255)
                                                                                                 aCATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCa                                                   >  1:3976680/1‑150 (MQ=255)
                                                                                                    ttACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGc                                                >  2:1708037/1‑150 (MQ=255)
                                                                                                    ttACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGc                                                <  2:506530/150‑1 (MQ=255)
                                                                                                    ttACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGc                                                >  2:1933189/1‑150 (MQ=255)
                                                                                                      aCCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTc                                              >  1:2198397/1‑150 (MQ=255)
                                                                                                       ccTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCg                                             <  1:1068854/150‑1 (MQ=255)
                                                                                                        cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGt                                            >  1:1342149/1‑150 (MQ=255)
                                                                                                          tgagccgccgccGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGa                                          <  1:241998/147‑1 (MQ=255)
                                                                                                             gccgccgccGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAAc                                       <  2:4270284/150‑1 (MQ=255)
                                                                                                                  cgccgGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATtgt                                  >  1:1545007/1‑150 (MQ=255)
                                                                                                                   gccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATtgtg                                 <  2:1199235/150‑1 (MQ=255)
                                                                                                                    ccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTgtgt                                <  1:4256427/150‑1 (MQ=255)
                                                                                                                     cgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTc                               <  1:4304288/150‑1 (MQ=255)
                                                                                                                          cAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCg                          <  1:1624740/150‑1 (MQ=255)
                                                                                                                           aGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGt                         >  2:3929144/1‑150 (MQ=255)
                                                                                                                           aGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGt                         >  2:4035364/1‑150 (MQ=255)
                                                                                                                             tACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTg                       >  1:2653051/1‑150 (MQ=255)
                                                                                                                               cGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGc                     >  1:3425622/1‑150 (MQ=255)
                                                                                                                                   gTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCatat                 <  2:12807/150‑1 (MQ=255)
                                                                                                                                   gTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCCTTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACAGTGTGTCAGGCGTTGGCatat                 >  1:2140185/1‑150 (MQ=255)
                                                                                                                                    tGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCAtata                <  1:3035959/150‑1 (MQ=255)
                                                                                                                                      aCCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGc              >  1:1781969/1‑150 (MQ=255)
                                                                                                                                      aCCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGc              <  1:982663/150‑1 (MQ=255)
                                                                                                                                      aCCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGc              <  2:283598/150‑1 (MQ=255)
                                                                                                                                       ccTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCg             >  2:1628209/1‑150 (MQ=255)
                                                                                                                                       ccTGGCGTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCg             <  2:3061996/150‑1 (MQ=255)
                                                                                                                                             tttCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCAt       >  1:4529283/1‑150 (MQ=255)
                                                                                                                                              ttCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATc      >  1:2452824/1‑150 (MQ=255)
                                                                                                                                               tCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCg     <  1:3500248/150‑1 (MQ=255)
                                                                                                                                               tCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCg     <  1:3584189/150‑1 (MQ=255)
                                                                                                                                               tCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCg     >  2:1533038/1‑150 (MQ=255)
                                                                                                                                                cACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCgg    >  2:1136446/1‑150 (MQ=255)
                                                                                                                                                  cgcgAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCGGGc  <  1:1000200/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCGGGC  >  CP009273/1151603‑1151898

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: