Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 284,851 T→A 18.5% E8V (GAA→GTA)  yagI ← CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273284,8510TA18.5% 332.0 / 67.6 135E8V (GAA→GTA) yagICP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (50/60);  new base A (11/14);  total (61/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 135 total aligned reads displayed.

TTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAA  >  CP009273/284707‑284995
                                                                                                                                                |                                                                                                                                                
 tCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAAc                                                                                                                                              <  2:1808189/148‑1 (MQ=255)
 tCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAAc                                                                                                                                              <  2:448733/148‑1 (MQ=255)
    gATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCTGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGa                                                                                                                                           <  2:1182310/148‑1 (MQ=255)
    gATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGa                                                                                                                                           <  2:1894512/148‑1 (MQ=255)
    gATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGa                                                                                                                                           >  1:1370275/1‑148 (MQ=255)
    gATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGa                                                                                                                                           <  2:1894498/148‑1 (MQ=255)
       cGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTg                                                                                                                                        <  1:2215539/148‑1 (MQ=255)
        gATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGa                                                                                                                                       <  1:3204220/148‑1 (MQ=255)
        gATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGa                                                                                                                                       <  1:3123596/148‑1 (MQ=255)
         atatAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGaa                                                                                                                                      <  1:2314598/148‑1 (MQ=255)
         atatAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGaa                                                                                                                                      >  1:1760819/1‑148 (MQ=255)
         atatAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGAACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGaa                                                                                                                                      >  1:3243279/1‑148 (MQ=255)
           ataGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGaata                                                                                                                                    <  2:2514422/148‑1 (MQ=255)
              gCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATc                                                                                                                                 <  1:1049621/148‑1 (MQ=255)
              gCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATc                                                                                                                                 <  2:3044685/148‑1 (MQ=255)
                cGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCgg                                                                                                                               <  1:2250682/148‑1 (MQ=255)
                    aaGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATa                                                                                                                           <  1:529152/148‑1 (MQ=255)
                    aaGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATa                                                                                                                           <  2:1248847/148‑1 (MQ=255)
                    aaGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATa                                                                                                                           <  2:1856455/148‑1 (MQ=255)
                     aGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATaa                                                                                                                          >  2:2273308/1‑148 (MQ=255)
                     aGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATaa                                                                                                                          <  1:2880172/148‑1 (MQ=255)
                      gCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAAt                                                                                                                         >  1:277780/1‑148 (MQ=255)
                      gCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAAt                                                                                                                         >  1:2875427/1‑148 (MQ=255)
                         gCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCt                                                                                                                      <  2:2048682/148‑1 (MQ=255)
                         gCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCt                                                                                                                      <  2:2050125/148‑1 (MQ=255)
                             ggTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCACATCAGTTTGCTGATGTCGGTGATCTTAAGATCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCa                                                                                                                  >  1:2510330/1‑148 (MQ=255)
                              gTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCAt                                                                                                                 <  1:1042121/148‑1 (MQ=255)
                               tttttAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCagg                                                                                                                >  2:2441702/1‑146 (MQ=255)
                                  ttAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTcc                                                                                                             <  2:35593/148‑1 (MQ=255)
                                     gcagcGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACg                                                                                                          <  1:459559/148‑1 (MQ=255)
                                     gcagcGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACg                                                                                                          >  1:2846242/1‑148 (MQ=255)
                                      cagcGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGc                                                                                                         >  1:759245/1‑148 (MQ=255)
                                      cagcGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGc                                                                                                         <  1:560061/148‑1 (MQ=255)
                                        gcGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTg                                                                                                       <  1:2452469/148‑1 (MQ=255)
                                         cGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTgg                                                                                                      >  1:1767110/1‑148 (MQ=255)
                                          gAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGa                                                                                                     <  2:3129507/148‑1 (MQ=255)
                                            gTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATTATGCTGGCATGTCCACGCTGGAAt                                                                                                   <  1:2734851/148‑1 (MQ=255)
                                            gTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAAt                                                                                                   <  1:1078694/148‑1 (MQ=255)
                                             tGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGATGATGTCGGTGATCTTAAGCTCGGTGGACTGCTCGTTGAACAGGTCGAGGATCTGCAACTCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAAtt                                                                                                  >  2:140797/1‑148 (MQ=255)
                                                aGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAg                                                                                               >  1:235552/1‑148 (MQ=255)
                                                aGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAg                                                                                               >  1:235551/1‑148 (MQ=255)
                                                  ggTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGc                                                                                             >  1:3261263/1‑148 (MQ=255)
                                                  ggTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGc                                                                                             >  1:3261459/1‑148 (MQ=255)
                                                  ggTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATAAAGGc                                                                                             >  2:79981/1‑148 (MQ=255)
                                                     aCTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGaa                                                                                          >  2:512601/1‑148 (MQ=255)
                                                     aCTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGaa                                                                                          >  1:311467/1‑148 (MQ=255)
                                                         ttGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACg                                                                                      <  2:1722126/148‑1 (MQ=255)
                                                           gCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTa                                                                                    <  1:2157130/148‑1 (MQ=255)
                                                            cTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTTGAATTAAGGCGAAAACGTAc                                                                                   <  2:3104805/148‑1 (MQ=255)
                                                            cTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTAc                                                                                   <  2:2387887/148‑1 (MQ=255)
                                                             tCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTAcc                                                                                  <  2:1693310/148‑1 (MQ=255)
                                                                  cccATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCttttt                                                                             <  1:2348402/148‑1 (MQ=255)
                                                                  cccATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCttttt                                                                             >  2:2615350/1‑148 (MQ=255)
                                                                   ccATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGCCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTc                                                                            <  2:450781/148‑1 (MQ=255)
                                                                    cATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCg                                                                           <  1:3107072/148‑1 (MQ=255)
                                                                       cAGTTTGCTGATGGCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGa                                                                        >  2:624115/1‑148 (MQ=255)
                                                                        aGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGAt                                                                       <  2:338705/148‑1 (MQ=255)
                                                                          tttGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGa                                                                     >  2:299428/1‑148 (MQ=255)
                                                                          tttGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACCTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGa                                                                     >  1:782113/1‑148 (MQ=255)
                                                                            tGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAACACGTACCTTTTTCGGGATGaaa                                                                   >  2:802811/1‑148 (MQ=255)
                                                                                gATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATc                                                               >  1:2140022/1‑148 (MQ=255)
                                                                                  tGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCaa                                                             <  1:2592151/148‑1 (MQ=255)
                                                                                   gTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAAt                                                            <  2:1616411/148‑1 (MQ=255)
                                                                                   gTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAAt                                                            >  1:153317/1‑148 (MQ=255)
                                                                                    tCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATg                                                           <  2:1236697/148‑1 (MQ=255)
                                                                                       gTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGaaa                                                        <  1:546304/148‑1 (MQ=255)
                                                                                        tGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAAt                                                       >  1:2612598/1‑148 (MQ=255)
                                                                                         gATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATg                                                      <  1:2372801/148‑1 (MQ=255)
                                                                                           tCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGa                                                    <  1:1912801/148‑1 (MQ=255)
                                                                                              tAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGcc                                                 >  1:994987/1‑148 (MQ=255)
                                                                                              tAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGcc                                                 >  1:2890875/1‑148 (MQ=255)
                                                                                              tAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGcc                                                 <  1:2698436/148‑1 (MQ=255)
                                                                                                    cGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAAGGAAATGGAGCCGGTGtt                                           >  1:656116/1‑148 (MQ=255)
                                                                                                      gTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtct                                         <  2:2344143/148‑1 (MQ=255)
                                                                                                      gTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtct                                         <  2:277780/148‑1 (MQ=255)
                                                                                                      gTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtct                                         <  2:361770/148‑1 (MQ=255)
                                                                                                        ggCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtctcc                                       <  1:3298779/148‑1 (MQ=255)
                                                                                                          ccTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtctccct                                     >  2:1168745/1‑148 (MQ=255)
                                                                                                          ccTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtctccct                                     >  2:2988790/1‑148 (MQ=255)
                                                                                                              cTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtctccctctcc                                 >  2:1367245/1‑148 (MQ=255)
                                                                                                               tCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtctccctctccc                                <  2:315932/148‑1 (MQ=255)
                                                                                                                cGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAAGCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAATACGTACCTGTTTCGGGATGAAATCTCAAGGACGTGGAGCCGGTGTtctccctctccgt                               >  2:1846005/1‑146 (MQ=255)
                                                                                                                  ttGAACAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtgt                             >  1:1010371/1‑148 (MQ=255)
                                                                                                                     aaCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTggg                          >  2:1397400/1‑148 (MQ=255)
                                                                                                                      aCAGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGga                         <  1:787196/148‑1 (MQ=255)
                                                                                                                       cAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCACTGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgag                        >  1:2502480/1‑148 (MQ=255)
                                                                                                                       cAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgag                        <  2:2905110/148‑1 (MQ=255)
                                                                                                                       cAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgag                        <  1:1016552/148‑1 (MQ=255)
                                                                                                                       cAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgag                        <  1:1016590/148‑1 (MQ=255)
                                                                                                                          gTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                     <  2:2140544/148‑1 (MQ=255)
                                                                                                                          gTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                     <  1:1397400/148‑1 (MQ=255)
                                                                                                                          gTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGATGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                     <  1:2273308/148‑1 (MQ=255)
                                                                                                                          gTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                     <  1:387148/148‑1 (MQ=255)
                                                                                                                          gTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                     >  1:2333654/1‑148 (MQ=255)
                                                                                                                               ggATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                >  2:2141964/1‑148 (MQ=255)
                                                                                                                                gATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCggg               >  2:687315/1‑148 (MQ=255)
                                                                                                                                   cTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATTAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg            <  2:655675/148‑1 (MQ=255)
                                                                                                                                     gCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAACGCGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAg          >  2:1991401/1‑148 (MQ=255)
                                                                                                                                      cAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg         >  2:2588050/1‑148 (MQ=255)
                                                                                                                                             cGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCaa  >  2:1877443/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAA  >  CP009273/284707‑284995

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: