Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 288,249 T→A 100% T137T (ACA→ACT)  yagK ← CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273288,2490TA95.6% 719.9 / 16.8 225T137T (ACA→ACT) yagKCP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base T (5/5);  new base A (83/132);  total (88/137)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 228 total aligned reads displayed.

GGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGG  >  CP009273/288104‑288393
                                                                                                                                                 |                                                                                                                                                
ggTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCagt                                                                                                                                                <  2:2218913/148‑1 (MQ=255)
    gTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGatc                                                                                                                                            <  2:2494674/148‑1 (MQ=255)
     tCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGatca                                                                                                                                           >  2:27165/1‑148 (MQ=255)
      ccAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGatcat                                                                                                                                          <  1:2504317/148‑1 (MQ=255)
      ccAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGatcat                                                                                                                                          >  1:2393081/1‑148 (MQ=255)
        aGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATcc                                                                                                                                        <  1:229789/148‑1 (MQ=255)
          tCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCct                                                                                                                                      <  1:2772725/148‑1 (MQ=255)
              gTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTaa                                                                                                                                  <  1:1540303/148‑1 (MQ=255)
              gTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTaa                                                                                                                                  <  1:2084428/148‑1 (MQ=255)
               ttAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAg                                                                                                                                 <  2:2915099/148‑1 (MQ=255)
               ttAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAg                                                                                                                                 >  1:2829412/1‑148 (MQ=255)
                tAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAgg                                                                                                                                <  1:1319338/148‑1 (MQ=255)
                  aCAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGtt                                                                                                                              <  2:2884284/148‑1 (MQ=255)
                  aCAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGtt                                                                                                                              >  1:2407272/1‑148 (MQ=255)
                      cgcCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCg                                                                                                                          <  1:104126/148‑1 (MQ=255)
                        ccTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTc                                                                                                                        <  1:1927394/148‑1 (MQ=255)
                            gTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGa                                                                                                                    <  1:2541266/148‑1 (MQ=255)
                             tAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGAt                                                                                                                   >  2:2022465/1‑148 (MQ=255)
                              aGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATc                                                                                                                  <  2:1153423/148‑1 (MQ=255)
                               gTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCg                                                                                                                 <  1:1932120/148‑1 (MQ=255)
                                 tCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTa                                                                                                               >  1:520910/1‑148 (MQ=255)
                                  ctgctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTaa                                                                                                              <  1:1245712/148‑1 (MQ=255)
                                  ctgctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTaa                                                                                                              <  2:2113802/148‑1 (MQ=255)
                                   tgctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAAt                                                                                                             <  2:528465/148‑1 (MQ=255)
                                   tgctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAAt                                                                                                             <  2:181366/148‑1 (MQ=255)
                                     ctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATcc                                                                                                           >  1:411897/1‑148 (MQ=255)
                                       ggAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATcccc                                                                                                         >  2:2940581/1‑148 (MQ=255)
                                        gAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCg                                                                                                        <  2:22020/148‑1 (MQ=255)
                                            tCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggt                                                                                                    >  2:2634365/1‑148 (MQ=255)
                                                tACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgta                                                                                                >  2:3295078/1‑148 (MQ=255)
                                                tACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgta                                                                                                >  2:369185/1‑148 (MQ=255)
                                                 aCTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgtag                                                                                               >  1:1372392/1‑148 (MQ=255)
                                                 aCTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgtag                                                                                               >  1:870015/1‑148 (MQ=255)
                                                 aCTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgtag                                                                                               >  2:951327/1‑148 (MQ=255)
                                                  cTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgtagt                                                                                              >  1:635170/1‑148 (MQ=255)
                                                   tATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgtagta                                                                                             >  1:3003077/1‑148 (MQ=255)
                                                    aTTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTaa                                                                                            >  2:598426/1‑148 (MQ=255)
                                                     ttGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAg                                                                                           >  2:675705/1‑148 (MQ=255)
                                                     ttGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAg                                                                                           >  2:1206376/1‑148 (MQ=255)
                                                         tgtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTc                                                                                       <  2:1975666/148‑1 (MQ=255)
                                                          gtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCt                                                                                      >  2:1634449/1‑148 (MQ=255)
                                                          gtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCt                                                                                      >  2:1355148/1‑148 (MQ=255)
                                                              aGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttgt                                                                                  <  1:3200877/148‑1 (MQ=255)
                                                              aGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttgt                                                                                  >  1:2802312/1‑148 (MQ=255)
                                                                cacaTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttgttg                                                                                <  2:2995865/148‑1 (MQ=255)
                                                                 acaTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGa                                                                               <  2:155314/148‑1 (MQ=255)
                                                                     aCTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaacaa                                                                           <  1:1147365/148‑1 (MQ=255)
                                                                     aCTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaacaa                                                                           >  1:1918138/1‑148 (MQ=255)
                                                                      cTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaacaac                                                                          >  2:1440119/1‑148 (MQ=255)
                                                                       ttACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAacaaca                                                                         >  2:3090261/1‑148 (MQ=255)
                                                                         aCAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAgg                                                                       <  1:917611/148‑1 (MQ=255)
                                                                            gTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCaa                                                                    <  2:682140/148‑1 (MQ=255)
                                                                            gTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCaa                                                                    <  2:2342710/148‑1 (MQ=255)
                                                                                tCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAtgt                                                                >  2:572326/1‑148 (MQ=255)
                                                                                tCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAtgt                                                                <  2:413885/148‑1 (MQ=255)
                                                                                tCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAtgt                                                                <  2:2077046/148‑1 (MQ=255)
                                                                                 cAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAtgtg                                                               <  1:949762/148‑1 (MQ=255)
                                                                                 cAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAtgtg                                                               <  1:1614835/148‑1 (MQ=255)
                                                                                  aGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGa                                                              <  2:120163/148‑1 (MQ=255)
                                                                                    gAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATa                                                            >  1:246728/1‑148 (MQ=255)
                                                                                         tGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGc                                                       >  1:2865915/1‑148 (MQ=255)
                                                                                         tGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGc                                                       >  1:2643119/1‑148 (MQ=255)
                                                                                          ggACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCa                                                      <  2:3206876/148‑1 (MQ=255)
                                                                                           gaccagaccaGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCAt                                                     >  2:128255/1‑148 (MQ=255)
                                                                                               agaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAc                                                 <  1:625181/148‑1 (MQ=255)
                                                                                               agaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAc                                                 <  2:337194/148‑1 (MQ=255)
                                                                                                gaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAcc                                                >  1:2301126/1‑148 (MQ=255)
                                                                                                 accagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCt                                               >  2:3233086/1‑148 (MQ=255)
                                                                                                     gggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAtt                                           <  2:174248/148‑1 (MQ=255)
                                                                                                      ggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTc                                          >  1:2034017/1‑148 (MQ=255)
                                                                                                         gATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGa                                       >  2:582157/1‑148 (MQ=255)
                                                                                                          aTCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGaa                                      >  1:28195/1‑148 (MQ=255)
                                                                                                          aTCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGaa                                      >  2:1607351/1‑148 (MQ=255)
                                                                                                          aTCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGaa                                      >  2:3060168/1‑148 (MQ=255)
                                                                                                           tCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAAt                                     >  2:180116/1‑148 (MQ=255)
                                                                                                               ctctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATActc                                 <  1:1867877/148‑1 (MQ=255)
                                                                                                               ctctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATActc                                 >  2:919504/1‑148 (MQ=255)
                                                                                                                 ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCtt                               >  2:1850940/1‑148 (MQ=255)
                                                                                                                 ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCtt                               <  2:868462/148‑1 (MQ=255)
                                                                                                                 ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCtt                               <  1:305455/148‑1 (MQ=255)
                                                                                                                   ctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTc                             >  1:702670/1‑148 (MQ=255)
                                                                                                                      cAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGcc                          >  2:2784439/1‑148 (MQ=255)
                                                                                                                        gTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCa                        <  1:50940/148‑1 (MQ=255)
                                                                                                                         tctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCat                       <  2:882504/148‑1 (MQ=255)
                                                                                                                          ctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCata                      <  1:1206376/148‑1 (MQ=255)
                                                                                                                          ctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCata                      >  2:1367425/1‑148 (MQ=255)
                                                                                                                              gggCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataat                  >  1:3092012/1‑148 (MQ=255)
                                                                                                                              gggCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataat                  >  2:18928/1‑148 (MQ=255)
                                                                                                                               ggCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataata                 <  1:1595733/148‑1 (MQ=255)
                                                                                                                                 cgcTGTACCAGGCTCCAGTGATCATACCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATACAGAATACTCTTTCGCCCATATaataaa               >  2:716515/1‑148 (MQ=255)
                                                                                                                                    tGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAataaataa            <  1:2118143/148‑1 (MQ=255)
                                                                                                                                     gTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAAtaaataaa           <  2:2393081/148‑1 (MQ=255)
                                                                                                                                      tACCAGGCTCCAGTGATCATCCCTCTAATGTTATCGTGCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAg          <  1:1747069/148‑1 (MQ=255)
                                                                                                                                      tACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCAGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACATGATTCAGAATACTCTTTAGCCCATATAAtaaataaat          >  1:1492947/1‑147 (MQ=255)
                                                                                                                                        ccAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAggg        >  2:1879927/1‑148 (MQ=255)
                                                                                                                                        ccAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAggg        <  2:448678/148‑1 (MQ=255)
                                                                                                                                         cAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGc       >  1:2704071/1‑148 (MQ=255)
                                                                                                                                          aGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCa      <  2:1691093/148‑1 (MQ=255)
                                                                                                                                            gCTCCAGTGATCATCCCTCTAAGCATATCGTCCTGATCGTAATCCCCGAGGTGGTAGTACGCGTCTTTGTTGAACCACAGGCAAATGTGATAGTGGAATAAACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAAc    >  2:1299871/1‑148 (MQ=255)
                                                                                                                                              tCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACgg  <  1:24427/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGG  >  CP009273/288104‑288393

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: ‑