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,2490TA99.1% 1091.2 / ‑1.2 329T137T (ACA→ACT) yagKCP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base A (152/174);  total (154/175)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 335 total aligned reads displayed.

CAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGT  >  CP009273/288102‑288394
                                                                                                                                                   |                                                                                                                                                 
  ggTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGa                                                                                                                                               >  2:5019413/1‑150 (MQ=255)
   gTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCTGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGat                                                                                                                                              >  2:421426/1‑150 (MQ=255)
    tAGTCCAGTCGGGTTAACAGCGCCTGGTAGGTCTGCTGGTAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGatc                                                                                                                                             <  1:5327924/150‑1 (MQ=255)
      gTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGatcat                                                                                                                                           >  2:4525450/1‑150 (MQ=255)
       tCCAGTCGTGTTAACAGCGACTGGTAGTTCTGCTGGAAGTCATTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGAACAGACCAGTTTGATAGTATCTCTCCAGTCTCATGGCGATGTACCAGGATACAGTGatcatc                                                                                                                                          <  2:1427199/150‑1 (MQ=255)
        ccAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATcc                                                                                                                                         >  1:1877330/1‑150 (MQ=255)
          aGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCct                                                                                                                                       >  2:1969286/1‑150 (MQ=255)
          aGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCct                                                                                                                                       >  2:4928091/1‑150 (MQ=255)
           gTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCctc                                                                                                                                      <  1:5265496/150‑1 (MQ=255)
            tCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCctct                                                                                                                                     <  2:3275852/150‑1 (MQ=255)
                gTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAgg                                                                                                                                 >  2:736857/1‑150 (MQ=255)
                    aCAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTAt                                                                                                                             <  2:5268279/150‑1 (MQ=255)
                     cAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATc                                                                                                                            <  2:2358093/150‑1 (MQ=255)
                        cgcCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTc                                                                                                                         >  1:630685/1‑150 (MQ=255)
                          ccTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCt                                                                                                                       <  2:624456/150‑1 (MQ=255)
                           cTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTg                                                                                                                      >  1:4552487/1‑150 (MQ=255)
                           cTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGAGCATCCCTCTAAGGTTATCGTCCTg                                                                                                                      >  2:4440747/1‑150 (MQ=255)
                            tGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGa                                                                                                                     >  1:3792278/1‑150 (MQ=255)
                             ggTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGAt                                                                                                                    >  1:4914078/1‑150 (MQ=255)
                             ggTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGAt                                                                                                                    >  1:3659105/1‑150 (MQ=255)
                              gTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATc                                                                                                                   <  2:2844957/150‑1 (MQ=255)
                              gTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGAACCAGGCTCCAGTGATCATCCATCTAAGGTTATCGTCATGATc                                                                                                                   >  2:3734286/1‑150 (MQ=255)
                                aGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGt                                                                                                                 >  2:3789882/1‑150 (MQ=255)
                                    ctgctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATc                                                                                                             <  2:4657103/150‑1 (MQ=255)
                                    ctgctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATc                                                                                                             <  1:724338/150‑1 (MQ=255)
                                      gctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATccc                                                                                                           <  1:1047869/150‑1 (MQ=255)
                                           aaGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAgg                                                                                                      >  2:2715547/1‑150 (MQ=255)
                                            aGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggt                                                                                                     >  2:2669675/1‑150 (MQ=255)
                                            aGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggt                                                                                                     <  2:2293726/150‑1 (MQ=255)
                                             gTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggtg                                                                                                    <  2:2039324/150‑1 (MQ=255)
                                                aGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgta                                                                                                 >  1:4971757/1‑150 (MQ=255)
                                                     tATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAg                                                                                            >  2:296101/1‑150 (MQ=255)
                                                     tATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAg                                                                                            <  1:3416460/150‑1 (MQ=255)
                                                      aTTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGc                                                                                           >  2:4322246/1‑150 (MQ=255)
                                                        tGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGGTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGt                                                                                         >  2:362155/1‑150 (MQ=255)
                                                            gtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCttt                                                                                     >  1:299112/1‑150 (MQ=255)
                                                               tAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttgtt                                                                                  >  2:2923945/1‑150 (MQ=255)
                                                               tAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttgtt                                                                                  <  2:3985685/150‑1 (MQ=255)
                                                                 gCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGa                                                                                <  1:1624630/150‑1 (MQ=255)
                                                                  cacaTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaa                                                                               >  1:440865/1‑150 (MQ=255)
                                                                  aacaTAATTACAGTTTTAAGGAAAATGGACCAGACCATGGTGATCGTATCTCTCCAGTCTCAGGGCGCTATACCAGGCTCCAGTGATCATCCATCTAAGGTTATAGTCCTGATCGTAAACAACGAGGTGGTAGTAAGCGTCTTGGTTGaa                                                                               >  2:1004818/2‑150 (MQ=255)
                                                                    caTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaaca                                                                             >  2:2400864/1‑150 (MQ=255)
                                                                         ttACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAgg                                                                        <  1:3180712/150‑1 (MQ=255)
                                                                            cAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCaa                                                                     >  1:1361154/1‑150 (MQ=255)
                                                                             aGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCaaa                                                                    >  1:1231718/1‑150 (MQ=255)
                                                                             aGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCaaa                                                                    >  1:5079642/1‑150 (MQ=255)
                                                                                 ttCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAtgtg                                                                <  2:3265858/150‑1 (MQ=255)
                                                                                      gAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGt                                                           <  1:4467001/150‑1 (MQ=255)
                                                                                       aaaaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTg                                                          >  2:1541847/1‑150 (MQ=255)
                                                                                       aaaaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTg                                                          >  2:2592112/1‑150 (MQ=255)
                                                                                         aaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGc                                                        <  1:921327/150‑1 (MQ=255)
                                                                                            ggACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCAtt                                                     <  1:2293427/150‑1 (MQ=255)
                                                                                             gaccagaccaGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCAttt                                                    <  2:19621/150‑1 (MQ=255)
                                                                                                cagaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACACCAGGCAAATGTGATAGTGGCATTTAcc                                                 >  1:5134454/1‑150 (MQ=255)
                                                                                                cagaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAcc                                                 <  2:3638871/150‑1 (MQ=255)
                                                                                                  gaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTg                                               <  1:4113480/150‑1 (MQ=255)
                                                                                                  gaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTg                                               <  2:2633063/150‑1 (MQ=255)
                                                                                                  gaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTg                                               <  1:4225739/150‑1 (MQ=255)
                                                                                                   accagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGa                                              <  2:5155770/150‑1 (MQ=255)
                                                                                                    ccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAt                                             <  2:3792287/150‑1 (MQ=255)
                                                                                                    ccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAt                                             <  1:2894165/150‑1 (MQ=255)
                                                                                                    ccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAt                                             <  1:2623421/150‑1 (MQ=255)
                                                                                                     cagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAtt                                            <  1:1958264/150‑1 (MQ=255)
                                                                                                     cagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGGACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAtt                                            <  2:1790721/150‑1 (MQ=255)
                                                                                                      agGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTc                                           >  1:4148031/1‑150 (MQ=255)
                                                                                                        ggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAg                                         <  1:2155078/150‑1 (MQ=255)
                                                                                                         gTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGa                                        >  2:1202518/1‑150 (MQ=255)
                                                                                                           gATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAAt                                      >  1:5170692/1‑150 (MQ=255)
                                                                                                           gATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAAt                                      <  1:2520191/150‑1 (MQ=255)
                                                                                                             tCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCAGCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATAc                                    <  2:797196/150‑1 (MQ=255)
                                                                                                              cGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATAct                                   >  1:3379387/1‑150 (MQ=255)
                                                                                                              cGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATAct                                   <  1:1840276/150‑1 (MQ=255)
                                                                                                               gTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATActc                                  <  2:4794717/150‑1 (MQ=255)
                                                                                                               gTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATActc                                  >  2:1663378/1‑150 (MQ=255)
                                                                                                                 ctctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCtt                                >  1:594616/1‑150 (MQ=255)
                                                                                                                 ctctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCtt                                >  1:742103/1‑150 (MQ=255)
                                                                                                                   ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTc                              >  2:5184433/1‑150 (MQ=255)
                                                                                                                   ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTc                              <  2:4914939/150‑1 (MQ=255)
                                                                                                                    tctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCg                             <  2:4459808/150‑1 (MQ=255)
                                                                                                                    tctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCg                             >  2:2543517/1‑150 (MQ=255)
                                                                                                                    tctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCg                             <  1:1207459/150‑1 (MQ=255)
                                                                                                                      tcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGcc                           <  2:1307211/150‑1 (MQ=255)
                                                                                                                        cAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCa                         >  1:4036718/1‑150 (MQ=255)
                                                                                                                         aGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCat                        <  1:2085158/150‑1 (MQ=255)
                                                                                                                         aGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCat                        >  1:4673148/1‑150 (MQ=255)
                                                                                                                            ctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGGGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCAtata                     >  1:3867211/1‑150 (MQ=255)
                                                                                                                            ctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTNTCGCCCAtata                     <  1:2766072/150‑1 (MQ=255)
                                                                                                                              cAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataat                   <  1:3521/150‑1 (MQ=255)
                                                                                                                               aGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataata                  >  2:3059190/1‑150 (MQ=255)
                                                                                                                               aGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataata                  <  2:4689206/150‑1 (MQ=255)
                                                                                                                                  gcgcTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaat               <  1:3739133/150‑1 (MQ=255)
                                                                                                                                      tGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAg           >  2:4113262/1‑150 (MQ=255)
                                                                                                                                      tGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAg           <  2:630743/150‑1 (MQ=255)
                                                                                                                                       gTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAgg          <  2:1279898/150‑1 (MQ=255)
                                                                                                                                         aCCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGc        <  1:756085/150‑1 (MQ=255)
                                                                                                                                         aCCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGc        <  2:91477/150‑1 (MQ=255)
                                                                                                                                          ccAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCa       <  1:4353025/150‑1 (MQ=255)
                                                                                                                                           cAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCaa      <  2:4268801/150‑1 (MQ=255)
                                                                                                                                            aGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAAc     <  1:1522948/150‑1 (MQ=255)
                                                                                                                                               cTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGt  <  2:3508927/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGT  >  CP009273/288102‑288394

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