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.3% 694.0 / ‑4.5 206T137T (ACA→ACT) yagKCP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (83/121);  total (85/121)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 206 total aligned reads displayed.

GGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTA  >  CP009273/288104‑288395
                                                                                                                                                 |                                                                                                                                                  
ggTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGa                                                                                                                                                >  2:4476329/1‑150 (MQ=255)
 gTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGat                                                                                                                                               <  2:2076778/150‑1 (MQ=255)
   aGTCCAGTAGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGatca                                                                                                                                             <  2:1994221/150‑1 (MQ=255)
    gTCCAGTCGGGTTAACAGCACCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGatcat                                                                                                                                            <  2:2582280/150‑1 (MQ=255)
      ccAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATcc                                                                                                                                          >  2:223885/1‑150 (MQ=255)
      ccAGTCGGGTTAACAGCGCCTGGTAGTTCGGCTGGAAGTCAGTACTAGTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATcc                                                                                                                                          >  2:2553159/1‑150 (MQ=255)
          tCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCctct                                                                                                                                      <  2:4802436/150‑1 (MQ=255)
          tCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCctct                                                                                                                                      >  2:4780009/1‑150 (MQ=255)
            gggTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTaa                                                                                                                                    <  2:590970/150‑1 (MQ=255)
               ttAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTCAGGt                                                                                                                                 >  2:4338306/1‑150 (MQ=255)
                 aaCAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTTTTGGTGTCTAGCACATACTTACTGTTTTCAGGAAAATGGAGCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGATGTACCAGTCTCCAGTGATCATCCCTCTAATGTTa                                                                                                                               <  2:2742243/150‑1 (MQ=255)
                  aCAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTAt                                                                                                                              <  1:2611431/150‑1 (MQ=255)
                     gcgcCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGt                                                                                                                           >  1:486051/1‑150 (MQ=255)
                      cgcCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTc                                                                                                                          <  2:116422/150‑1 (MQ=255)
                      cgcCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTc                                                                                                                          <  2:4536457/150‑1 (MQ=255)
                           ggTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGAt                                                                                                                     <  2:304796/150‑1 (MQ=255)
                           ggTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGAt                                                                                                                     <  2:2499846/150‑1 (MQ=255)
                           ggTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGAt                                                                                                                     <  2:4966428/150‑1 (MQ=255)
                           ggTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGAt                                                                                                                     >  2:1380104/1‑150 (MQ=255)
                               gTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTa                                                                                                                 >  2:1858927/1‑150 (MQ=255)
                                    gctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATccc                                                                                                            <  2:638467/150‑1 (MQ=255)
                                      tgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCg                                                                                                          <  2:1454044/150‑1 (MQ=255)
                                      tgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCg                                                                                                          <  2:4059372/150‑1 (MQ=255)
                                      tgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCg                                                                                                          <  1:3253000/150‑1 (MQ=255)
                                       ggAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGa                                                                                                         >  1:4343479/1‑150 (MQ=255)
                                        gAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATAGTAATCCCCGAg                                                                                                        <  1:2327593/150‑1 (MQ=255)
                                          aGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggt                                                                                                      >  2:3702309/1‑150 (MQ=255)
                                           gTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggtg                                                                                                     >  1:1310206/1‑150 (MQ=255)
                                            tCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggtgg                                                                                                    >  2:4979477/1‑150 (MQ=255)
                                             cAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggtggt                                                                                                   <  2:17140/150‑1 (MQ=255)
                                              aGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgta                                                                                                  >  2:269897/1‑150 (MQ=255)
                                               gTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCCGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgtag                                                                                                 <  2:802931/150‑1 (MQ=255)
                                                  cTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTaa                                                                                              <  1:3543321/150‑1 (MQ=255)
                                                  cTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTaa                                                                                              <  1:963109/150‑1 (MQ=255)
                                                  cTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTaa                                                                                              <  1:3456504/150‑1 (MQ=255)
                                                     ttGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCg                                                                                           >  1:2841092/1‑150 (MQ=255)
                                                     ttGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCg                                                                                           <  1:125321/150‑1 (MQ=255)
                                                     ttGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCg                                                                                           <  1:4416619/150‑1 (MQ=255)
                                                      tGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGt                                                                                          >  1:2159695/1‑150 (MQ=255)
                                                      tGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGt                                                                                          >  1:3692941/1‑150 (MQ=255)
                                                          gtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCttt                                                                                      <  1:396117/150‑1 (MQ=255)
                                                            cTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttgt                                                                                    <  1:1052135/150‑1 (MQ=255)
                                                                cacaTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaa                                                                                >  1:3152342/1‑150 (MQ=255)
                                                                      cTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAg                                                                          <  1:4866659/150‑1 (MQ=255)
                                                                        tACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGc                                                                        <  1:2566759/150‑1 (MQ=255)
                                                                        tACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGc                                                                        <  1:2516241/150‑1 (MQ=255)
                                                                          cAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCaa                                                                      <  2:4590683/150‑1 (MQ=255)
                                                                           aGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCaaa                                                                     >  1:2664388/1‑150 (MQ=255)
                                                                              tttCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAtgt                                                                  <  2:4794927/150‑1 (MQ=255)
                                                                               ttCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAtgtg                                                                 <  1:4682221/150‑1 (MQ=255)
                                                                                   ggAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAg                                                             <  2:1756577/150‑1 (MQ=255)
                                                                                   ggAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAg                                                             >  2:1277146/1‑150 (MQ=255)
                                                                                   ggAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAg                                                             <  2:3956511/150‑1 (MQ=255)
                                                                                    gAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGt                                                            >  1:4120951/1‑150 (MQ=255)
                                                                                    gAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGt                                                            <  2:1945539/150‑1 (MQ=255)
                                                                                     aaaaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTg                                                           <  2:3792732/150‑1 (MQ=255)
                                                                                      aaaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTgg                                                          <  1:2539637/150‑1 (MQ=255)
                                                                                       aaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGc                                                         >  2:3552671/1‑150 (MQ=255)
                                                                                          ggACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCAtt                                                      >  1:311606/1‑150 (MQ=255)
                                                                                           gaccagaccaGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCAttt                                                     <  2:486119/150‑1 (MQ=255)
                                                                                            accagaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTa                                                    <  2:1830652/150‑1 (MQ=255)
                                                                                              cagaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAcc                                                  >  1:3855263/1‑150 (MQ=255)
                                                                                              cagaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAcc                                                  <  2:2874534/150‑1 (MQ=255)
                                                                                               agaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCt                                                 >  1:555093/1‑150 (MQ=255)
                                                                                               agaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCt                                                 >  1:4290336/1‑150 (MQ=255)
                                                                                                gaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTg                                                <  2:3257214/150‑1 (MQ=255)
                                                                                                gaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTg                                                <  2:1646995/150‑1 (MQ=255)
                                                                                                  ccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAt                                              <  1:634531/150‑1 (MQ=255)
                                                                                                   cagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAtt                                             <  1:1730728/150‑1 (MQ=255)
                                                                                                   cagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAtt                                             <  2:2643626/150‑1 (MQ=255)
                                                                                                    agGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTc                                            >  1:992323/1‑150 (MQ=255)
                                                                                                     gggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCa                                           <  1:1124743/150‑1 (MQ=255)
                                                                                                     gggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCa                                           <  2:2116531/150‑1 (MQ=255)
                                                                                                       gTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGa                                         <  2:4905527/150‑1 (MQ=255)
                                                                                                       gTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGa                                         >  2:2188676/1‑150 (MQ=255)
                                                                                                            cGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATAct                                    <  2:3128179/150‑1 (MQ=255)
                                                                                                               ctctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCtt                                 <  1:5013353/150‑1 (MQ=255)
                                                                                                                 ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTc                               <  1:2243821/150‑1 (MQ=255)
                                                                                                                 ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTc                               <  2:2793371/150‑1 (MQ=255)
                                                                                                                  tctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCg                              <  1:1547408/150‑1 (MQ=255)
                                                                                                                  tctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCg                              <  2:4482056/150‑1 (MQ=255)
                                                                                                                     ccAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGccc                           <  1:1944864/150‑1 (MQ=255)
                                                                                                                      cAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCa                          >  2:4021094/1‑150 (MQ=255)
                                                                                                                       aGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCat                         <  2:1285690/150‑1 (MQ=255)
                                                                                                                        gTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCata                        >  1:2930744/1‑150 (MQ=255)
                                                                                                                        gTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCata                        <  1:1643156/150‑1 (MQ=255)
                                                                                                                        gTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCata                        <  1:4626972/150‑1 (MQ=255)
                                                                                                                         tctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCatat                       <  1:2781350/150‑1 (MQ=255)
                                                                                                                         tctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCatat                       <  1:4476329/150‑1 (MQ=255)
                                                                                                                          ctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCAtata                      <  2:4137316/150‑1 (MQ=255)
                                                                                                                          ctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCAtata                      >  1:3091153/1‑150 (MQ=255)
                                                                                                                          ctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCAtata                      <  1:2427154/150‑1 (MQ=255)
                                                                                                                          ctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCAtata                      >  1:4694355/1‑150 (MQ=255)
                                                                                                                            cAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataat                    <  1:4780009/150‑1 (MQ=255)
                                                                                                                             aGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataata                   <  1:4400767/150‑1 (MQ=255)
                                                                                                                              gggCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATAtaataa                  >  1:2047761/1‑150 (MQ=255)
                                                                                                                                       aCCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGc         <  2:3252159/150‑1 (MQ=255)
                                                                                                                                          aGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAAc      <  2:2129731/150‑1 (MQ=255)
                                                                                                                                            gCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACgg    >  1:3325353/1‑150 (MQ=255)
                                                                                                                                              tCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTa  <  1:3010613/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTA  >  CP009273/288104‑288395

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