Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 290,012 C→T 100% D167N (GAT→AAT)  yagM ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273290,0120CT96.2% 699.7 / 11.4 213D167N (GAT→AAT) yagMCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base T (87/117);  total (90/123)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 213 total aligned reads displayed.

TGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAACGAAC  >  CP009273/289866‑290158
                                                                                                                                                  |                                                                                                                                                  
  aCAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAAt                                                                                                                                               >  1:2810019/1‑150 (MQ=255)
   cAAGAGTAGAGACTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAAtt                                                                                                                                              <  1:1382310/150‑1 (MQ=255)
    aaGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAAttt                                                                                                                                             >  1:1419197/1‑150 (MQ=255)
     ggagTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAAtttt                                                                                                                                            >  2:1068248/2‑150 (MQ=255)
     agagTAGAGCCTGCCCGGCATCAAAGTAATACCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAAtttt                                                                                                                                            <  2:3401326/150‑1 (MQ=255)
       agTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTcc                                                                                                                                          <  2:473710/150‑1 (MQ=255)
         tAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCt                                                                                                                                        <  1:1274280/150‑1 (MQ=255)
         tAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCt                                                                                                                                        >  1:1879995/1‑150 (MQ=255)
          agagCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCtt                                                                                                                                       >  1:3145165/1‑150 (MQ=255)
           gagCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTa                                                                                                                                      >  1:3992254/1‑150 (MQ=255)
           gagCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTa                                                                                                                                      <  1:2198203/150‑1 (MQ=255)
            agCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAg                                                                                                                                     <  1:4144767/150‑1 (MQ=255)
                tGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAg                                                                                                                                 >  1:2916732/1‑150 (MQ=255)
                tGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAg                                                                                                                                 >  1:3467572/1‑150 (MQ=255)
                  cccGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGa                                                                                                                               >  2:1203516/1‑150 (MQ=255)
                       cATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGcc                                                                                                                          <  2:1041904/150‑1 (MQ=255)
                        aTCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGATTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCt                                                                                                                         <  2:2907614/150‑1 (MQ=255)
                         tCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCtt                                                                                                                        <  1:99641/150‑1 (MQ=255)
                           aaaGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTc                                                                                                                      >  2:2345028/1‑150 (MQ=255)
                           aaaGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTc                                                                                                                      <  2:167410/150‑1 (MQ=255)
                           aaaGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTc                                                                                                                      <  2:670094/150‑1 (MQ=255)
                           aaaGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTc                                                                                                                      <  2:670202/150‑1 (MQ=255)
                             aGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCAc                                                                                                                    <  2:1591442/150‑1 (MQ=255)
                             aGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCAc                                                                                                                    >  2:3465987/1‑150 (MQ=255)
                              gTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACg                                                                                                                   <  1:2234150/150‑1 (MQ=255)
                              gTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCANATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACg                                                                                                                   <  1:2681656/150‑1 (MQ=255)
                                    aCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAAc                                                                                                             >  2:4048187/1‑150 (MQ=255)
                                    aCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAAc                                                                                                             <  2:3640213/150‑1 (MQ=255)
                                      ccTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTg                                                                                                           >  1:1842993/1‑150 (MQ=255)
                                      ccTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCGTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTAGAACTg                                                                                                           >  2:2944605/1‑150 (MQ=255)
                                        tCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGtt                                                                                                         <  1:754167/150‑1 (MQ=255)
                                         cAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGttt                                                                                                        >  1:2827243/1‑150 (MQ=255)
                                         cAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGttt                                                                                                        <  2:1573827/150‑1 (MQ=255)
                                           aaGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGttttt                                                                                                      >  1:1516277/1‑150 (MQ=255)
                                            aGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGtttttt                                                                                                     <  2:4219997/150‑1 (MQ=255)
                                            aGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGtttttt                                                                                                     <  2:4220133/150‑1 (MQ=255)
                                             ggTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTc                                                                                                    <  1:5005867/150‑1 (MQ=255)
                                              cttgttCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCa                                                                                                   <  1:3894434/149‑1 (MQ=255)
                                               ttgttCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCAt                                                                                                  <  1:4062865/150‑1 (MQ=255)
                                                tgttCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATc                                                                                                 >  2:3711540/1‑150 (MQ=255)
                                                 gttCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCt                                                                                                <  2:4927562/150‑1 (MQ=255)
                                                     tATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTcc                                                                                            <  2:1879995/150‑1 (MQ=255)
                                                     tATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTcc                                                                                            <  1:1178809/150‑1 (MQ=255)
                                                       tGAATAAGAATAATGTTTTAATATCCTGAGACAATTTGATGCTTTACATCAAAATGTCTTTGATCGTAAGTTTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTAACGTATAACTGTTTTTTCATCTTTcccc                                                                                          <  1:887778/150‑1 (MQ=255)
                                                       tGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTcccc                                                                                          <  1:3541713/150‑1 (MQ=255)
                                                         aataagaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCAt                                                                                        <  2:280379/150‑1 (MQ=255)
                                                          ataagaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATg                                                                                       >  1:563309/1‑150 (MQ=255)
                                                             agaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAAt                                                                                    <  2:87087/150‑1 (MQ=255)
                                                               aataatGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATtc                                                                                  >  1:351761/1‑150 (MQ=255)
                                                                   atGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATtcatcg                                                                              <  2:3791713/150‑2 (MQ=255)
                                                                    tGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTcatcat                                                                             >  1:2080887/1‑150 (MQ=255)
                                                                     gTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCAtt                                                                            <  2:3672325/150‑1 (MQ=255)
                                                                        tGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGc                                                                         <  1:1223850/150‑1 (MQ=255)
                                                                        tGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGc                                                                         <  1:3727847/150‑1 (MQ=255)
                                                                          aaTATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGt                                                                       >  1:552952/1‑150 (MQ=255)
                                                                             atCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCa                                                                    >  1:2056844/1‑150 (MQ=255)
                                                                              tCCTGATACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAg                                                                   <  2:4844189/150‑1 (MQ=255)
                                                                               ccTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGt                                                                  >  2:111291/1‑150 (MQ=255)
                                                                               ccTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGt                                                                  >  2:811379/1‑150 (MQ=255)
                                                                                cTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGtt                                                                 >  1:4192127/1‑150 (MQ=255)
                                                                                 tGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTg                                                                >  1:2311932/1‑150 (MQ=255)
                                                                                 tGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTg                                                                >  1:2245647/1‑150 (MQ=255)
                                                                                 tGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTg                                                                <  1:4951125/150‑1 (MQ=255)
                                                                                  gagaCACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGa                                                               >  2:5012608/1‑150 (MQ=255)
                                                                                        cTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTaa                                                         >  1:4391721/1‑150 (MQ=255)
                                                                                        cTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTaa                                                         >  1:2936566/1‑150 (MQ=255)
                                                                                        cTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTaa                                                         >  2:2128135/1‑150 (MQ=255)
                                                                                        cTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTaa                                                         <  2:4580672/150‑1 (MQ=255)
                                                                                         tttGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAAt                                                        <  1:3061550/150‑1 (MQ=255)
                                                                                               gCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCaa                                                  <  2:332445/150‑1 (MQ=255)
                                                                                               gCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCaa                                                  <  2:598679/150‑1 (MQ=255)
                                                                                                cTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAAt                                                 <  2:3992254/150‑1 (MQ=255)
                                                                                                cTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAAt                                                 >  1:5142049/1‑150 (MQ=255)
                                                                                                  ttACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATgg                                               >  1:3234278/1‑150 (MQ=255)
                                                                                                    aCATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATgggg                                             <  2:3779200/150‑1 (MQ=255)
                                                                                                        cAAAATTTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGGTCAGTTGAAAGTAATCTGCAATGGGGAttt                                         <  1:5074871/150‑1 (MQ=255)
                                                                                                                cTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGa                                 >  1:706971/1‑150 (MQ=255)
                                                                                                                  ttGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACt                               <  1:3790145/150‑1 (MQ=255)
                                                                                                                    gATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCa                             >  2:1617370/1‑150 (MQ=255)
                                                                                                                     aTTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCat                            >  1:1069592/1‑150 (MQ=255)
                                                                                                                     aTTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCat                            >  1:4787282/1‑150 (MQ=255)
                                                                                                                      ttGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatg                           <  1:3808831/150‑1 (MQ=255)
                                                                                                                      ttGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatg                           <  2:1419197/150‑1 (MQ=255)
                                                                                                                        gtaaatgtaaatGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatgat                         <  2:410872/150‑1 (MQ=255)
                                                                                                                         taaatgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatgata                        <  2:5057544/150‑2 (MQ=255)
                                                                                                                          aaatgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGt                       <  2:900728/150‑1 (MQ=255)
                                                                                                                          aaatgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGt                       >  1:1924738/1‑150 (MQ=255)
                                                                                                                            atgtaaatgGCATTTACCATATTGTAATTTTCCCTTAGGGAGGTTAGTCTTTCACGTATAACTGTTTTTTCATCTTTCCCCCTGTATTCATCATTAGCGTTCATTTGAAAGTAATCTGCAATGGGGATTTTTACTGGACTCATTATGTTc                     <  2:1872457/150‑1 (MQ=255)
                                                                                                                             tgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCg                    <  2:1103835/150‑1 (MQ=255)
                                                                                                                              gtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGt                   >  1:4927393/1‑150 (MQ=255)
                                                                                                                               taaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGtt                  >  2:4148820/1‑150 (MQ=255)
                                                                                                                               taaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGtt                  <  2:226282/150‑1 (MQ=255)
                                                                                                                                   tgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTcc              >  2:1259022/1‑150 (MQ=255)
                                                                                                                                    ggCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCa             <  1:3873500/150‑1 (MQ=255)
                                                                                                                                     gCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAg            <  1:3465987/150‑1 (MQ=255)
                                                                                                                                      cTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGa           >  2:3187073/1‑150 (MQ=255)
                                                                                                                                         tttCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTa        <  1:4578869/150‑1 (MQ=255)
                                                                                                                                         tttCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTa        <  2:1826971/150‑1 (MQ=255)
                                                                                                                                          ttCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTaa       >  2:2400912/1‑150 (MQ=255)
                                                                                                                                           tCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAAc      <  1:5056576/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAACGAAC  >  CP009273/289866‑290158

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: