Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,151,750 G→T 91.6% 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,7500GT91.6% 876.6 / 59.1 311R178L (CGC→CTC) holBDNA polymerase III, delta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (19/9);  new base T (144/137);  total (165/146)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 317 total aligned reads displayed.

TGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCGG  >  CP009273/1151603‑1151896
                                                                                                                                                   |                                                                                                                                                  
tGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCActcg                                                                                                                                                  <  1:1927333/150‑1 (MQ=255)
tGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCActcg                                                                                                                                                  >  1:4409002/1‑150 (MQ=255)
tGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCAcgcg                                                                                                                                                  <  1:965580/150‑1 (MQ=255)
 gAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGa                                                                                                                                                 >  2:2251307/1‑150 (MQ=255)
    aaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGt                                                                                                                                              <  1:4188249/150‑1 (MQ=255)
       gCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGAc                                                                                                                                           >  2:1306081/1‑150 (MQ=255)
        cTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACa                                                                                                                                          <  1:115445/150‑1 (MQ=255)
        cTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACa                                                                                                                                          <  1:1433990/150‑1 (MQ=255)
          tGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAAt                                                                                                                                        <  1:2570635/150‑1 (MQ=255)
            aaGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGt                                                                                                                                      >  1:2280109/1‑150 (MQ=255)
                  caccaGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAgg                                                                                                                                >  1:1064221/1‑150 (MQ=255)
                     cagcagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATg                                                                                                                             >  2:3714522/1‑150 (MQ=255)
                      agcagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGc                                                                                                                            >  2:2197569/1‑150 (MQ=255)
                      agcagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGc                                                                                                                            >  1:362128/1‑150 (MQ=255)
                        cagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAt                                                                                                                          <  1:1679656/150‑1 (MQ=255)
                        cagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAt                                                                                                                          <  1:4241558/150‑1 (MQ=255)
                        cagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAt                                                                                                                          <  1:3585653/150‑1 (MQ=255)
                        cagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAAAGTCACACGATGCAt                                                                                                                          >  1:4645623/1‑150 (MQ=255)
                           aaaCTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAttcc                                                                                                                       <  2:473621/150‑3 (MQ=255)
                             aCTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACtt                                                                                                                     <  1:3784831/150‑1 (MQ=255)
                                  gTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGc                                                                                                                >  1:522210/1‑150 (MQ=255)
                                    ttttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCAt                                                                                                              >  2:3239551/1‑150 (MQ=255)
                                    ttttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCAt                                                                                                              <  2:3975822/150‑1 (MQ=255)
                                     tttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGGACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGAAGCCTTACATTCCGCaat                                                                                                             >  1:4443243/1‑148 (MQ=255)
                                     tttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCAtt                                                                                                             <  2:3739293/150‑1 (MQ=255)
                                      ttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTg                                                                                                            <  1:3229197/150‑1 (MQ=255)
                                        ccTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcg                                                                                                          <  1:3500150/150‑1 (MQ=255)
                                         cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGTGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                                                         >  2:4343797/1‑150 (MQ=255)
                                         cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                                                         >  1:720166/1‑150 (MQ=255)
                                          tGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCt                                                                                                        >  1:3854608/1‑150 (MQ=255)
                                               aCCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCCCAGGATGCATTACTTGCCGCATTGCGCTTAAgc                                                                                                   >  2:4571858/1‑150 (MQ=255)
                                                cccGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAgcg                                                                                                  <  1:1855385/150‑1 (MQ=255)
                                                     gAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGt                                                                                             >  1:299522/1‑150 (MQ=255)
                                                      aGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGtt                                                                                            >  2:1967017/1‑150 (MQ=255)
                                                       gCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTc                                                                                           >  1:2310622/1‑150 (MQ=255)
                                                       gCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTc                                                                                           <  2:4434300/150‑1 (MQ=255)
                                                       gCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCATTCGCCGCCGCCGGAACAGTAACCCGTGAACAGGCTTTAACTCGAAGTGACAAAGTCACACGATGCATTAATTTCCCCATTTCGCTTAATCGCCGGTTc                                                                                           >  1:4064708/1‑150 (MQ=255)
                                                         cTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTTACAAAGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGc                                                                                         >  1:4144209/1‑150 (MQ=255)
                                                         cTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGc                                                                                         >  1:4382954/1‑150 (MQ=255)
                                                         cTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCAGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGc                                                                                         >  1:1981331/1‑150 (MQ=255)
                                                           gAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCt                                                                                       <  2:4503901/150‑1 (MQ=255)
                                                            aaCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTg                                                                                      >  1:2627228/1‑150 (MQ=255)
                                                            aaCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTAACATGCGCCGCCCCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAAAGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTg                                                                                      >  1:1286348/1‑150 (MQ=255)
                                                             aCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTgg                                                                                     >  2:345335/1‑150 (MQ=255)
                                                                tttACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGgcgc                                                                                  >  2:4370444/1‑150 (MQ=255)
                                                                    cTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCa                                                                              <  2:1826199/150‑1 (MQ=255)
                                                                     tGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCAc                                                                             >  2:4196493/1‑150 (MQ=255)
                                                                     tGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCAc                                                                             >  2:2359327/1‑150 (MQ=255)
                                                                      ggCAACATTCCGTAGTCGTGGTCGGGTACAGTCCCGTGCGCCGCCGCCGGCACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTCCGTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACt                                                                            <  2:564207/150‑1 (MQ=255)
                                                                      ggCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACt                                                                            >  2:1661729/1‑150 (MQ=255)
                                                                       gCAACATTACGTAGTCGTTGTCGGTTACATTCCCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTg                                                                           <  1:1628171/150‑1 (MQ=255)
                                                                       gCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTg                                                                           <  1:393697/150‑1 (MQ=255)
                                                                       gCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTg                                                                           <  1:1629348/150‑1 (MQ=255)
                                                                          aCATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCg                                                                        >  2:430791/1‑150 (MQ=255)
                                                                              tactttgTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTGGAAGTGAAAATGTCACAGGATGCATTACTTGCCGAATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgttgt                                                                    <  1:4688622/144‑1 (MQ=255)
                                                                                cGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGttt                                                                  <  1:1019577/150‑1 (MQ=255)
                                                                                cGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGttt                                                                  <  1:939286/150‑1 (MQ=255)
                                                                                 gTAGTCGTTGTCGGTTACAATACCTTGCGCCGCCGCCCGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAAAGTCACAGGATGCATTAAATTCCGCATTGCGCATAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTc                                                                 >  2:1609214/1‑150 (MQ=255)
                                                                                  tAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCa                                                                >  1:1840405/1‑150 (MQ=255)
                                                                                   aGTCGTTGTCGGTTACATTAACTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAATCGCCGGTTCGCCTGGCGCGGCACAGGCGTTGTTTAAg                                                               >  1:4506576/1‑150 (MQ=255)
                                                                                        ttGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGgaga                                                          >  1:3308315/1‑150 (MQ=255)
                                                                                        ttGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGgaga                                                          >  2:3079650/1‑150 (MQ=255)
                                                                                          gTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATa                                                        >  1:769065/1‑150 (MQ=255)
                                                                                           tCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATaa                                                       >  2:1089042/1‑150 (MQ=255)
                                                                                            cGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAAc                                                      >  2:1155340/1‑150 (MQ=255)
                                                                                             ggTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACt                                                     >  1:2769758/1‑150 (MQ=255)
                                                                                               ttACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTgg                                                   <  1:1409462/150‑1 (MQ=255)
                                                                                                tACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGc                                                  >  2:3523669/1‑150 (MQ=255)
                                                                                                 aCATTACCTTGCGCCGCCGCCGGAGCAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCa                                                 <  2:1849176/150‑1 (MQ=255)
                                                                                                  cATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAg                                                <  2:2519360/150‑1 (MQ=255)
                                                                                                  cATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAg                                                >  2:3284521/1‑150 (MQ=255)
                                                                                                     tACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCt                                             <  1:1468675/150‑1 (MQ=255)
                                                                                                      aCCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTc                                            >  2:3995844/1‑150 (MQ=255)
                                                                                                        cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGt                                          <  2:2648337/150‑1 (MQ=255)
                                                                                                        cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGt                                          <  2:2632190/150‑1 (MQ=255)
                                                                                                         ttGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTg                                         <  2:3775099/150‑1 (MQ=255)
                                                                                                          tGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGa                                        >  1:1128773/1‑150 (MQ=255)
                                                                                                             gccgccgccGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAAc                                     <  2:2220892/150‑1 (MQ=255)
                                                                                                              ccgccgccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACa                                    <  1:3714800/150‑1 (MQ=255)
                                                                                                               cgccgccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACAt                                   >  1:4594394/1‑150 (MQ=255)
                                                                                                                   gccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATtgtg                               >  1:2381951/1‑150 (MQ=255)
                                                                                                                    ccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTgtgt                              >  1:88891/1‑150 (MQ=255)
                                                                                                                      ggAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATCACTGGCAGGCGCGTGAAACAGTGTGTCa                            >  1:2073386/1‑150 (MQ=255)
                                                                                                                      ggAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCa                            >  2:3879896/1‑150 (MQ=255)
                                                                                                                      ggAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCa                            >  2:1022907/1‑150 (MQ=255)
                                                                                                                      ggAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCa                            >  1:4751387/1‑150 (MQ=255)
                                                                                                                          cAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCg                        <  2:4645712/150‑1 (MQ=255)
                                                                                                                            gTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGtt                      <  1:2743108/150‑1 (MQ=255)
                                                                                                                            gTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGtt                      <  1:3108113/150‑1 (MQ=255)
                                                                                                                               cGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGc                   >  1:430550/1‑150 (MQ=255)
                                                                                                                                gCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCa                  >  2:2994846/1‑150 (MQ=255)
                                                                                                                                  cGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCata                <  2:2619079/150‑1 (MQ=255)
                                                                                                                                   gTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCatat               <  1:1258481/150‑1 (MQ=255)
                                                                                                                                    tGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCAtata              <  2:4117075/150‑1 (MQ=255)
                                                                                                                                      aCCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGc            >  2:703037/1‑150 (MQ=255)
                                                                                                                                         tGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTg         >  2:3733206/1‑150 (MQ=255)
                                                                                                                                           gCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGcc       >  2:914873/1‑150 (MQ=255)
                                                                                                                                            cTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCa      >  1:1769522/1‑150 (MQ=255)
                                                                                                                                             tttCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCAt     <  2:2773415/150‑1 (MQ=255)
                                                                                                                                             tttCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCATGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCAt     >  1:2497534/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCGG  >  CP009273/1151603‑1151896

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: