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,0120CT98.5% 1329.4 / 6.0 402D167N (GAT→AAT) yagMCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (205/191);  total (209/193)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 402 total aligned reads displayed.

GACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAA  >  CP009273/289867‑290153
                                                                                                                                                 |                                                                                                                                             
  cAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTaa                                                                                                                                           >  1:3437502/1‑148 (MQ=255)
    agagTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAAtt                                                                                                                                         <  1:327095/148‑1 (MQ=255)
    agagTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAAtt                                                                                                                                         <  2:2996441/148‑1 (MQ=255)
      agTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAAtttt                                                                                                                                       >  2:1163048/1‑148 (MQ=255)
        tAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTcc                                                                                                                                     <  2:152811/148‑1 (MQ=255)
          gagCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCt                                                                                                                                   >  1:88827/1‑148 (MQ=255)
             ccTGCACGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAg                                                                                                                                >  2:669361/1‑148 (MQ=255)
               tGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCg                                                                                                                              >  1:3576543/1‑148 (MQ=255)
                 cccGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAg                                                                                                                            <  2:3610090/148‑1 (MQ=255)
                  ccGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAgg                                                                                                                           <  1:3015391/148‑1 (MQ=255)
                     gCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATa                                                                                                                        >  2:3170565/1‑148 (MQ=255)
                      cATCAAAGTAATAACCCTCACAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAg                                                                                                                       <  1:3403942/148‑1 (MQ=255)
                        tCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGcc                                                                                                                     <  2:2667309/148‑1 (MQ=255)
                         cAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCt                                                                                                                    <  2:3528463/148‑1 (MQ=255)
                           aaGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCttt                                                                                                                  >  2:2189675/1‑148 (MQ=255)
                            aGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTc                                                                                                                 >  2:3227870/1‑148 (MQ=255)
                             gTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCa                                                                                                                >  1:2632860/1‑148 (MQ=255)
                             gTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTGTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCa                                                                                                                >  1:715948/1‑148 (MQ=255)
                               aataaCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACg                                                                                                              <  1:2167444/148‑1 (MQ=255)
                                ataaCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGt                                                                                                             >  2:805956/1‑148 (MQ=255)
                                   aCCCTGAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTGTAAATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGTTAGCCTTTCACGtata                                                                                                          <  2:1626092/148‑1 (MQ=255)
                                   aCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGtata                                                                                                          >  1:402355/1‑148 (MQ=255)
                                         aaaGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGtt                                                                                                    <  2:3559025/148‑1 (MQ=255)
                                          aaGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGttt                                                                                                   >  1:1048877/1‑148 (MQ=255)
                                           aGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGtttt                                                                                                  >  2:3132880/1‑148 (MQ=255)
                                                 ttCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATc                                                                                            >  1:1203926/1‑148 (MQ=255)
                                                 tactaacaatcagactaaTGTTTCAATAGCCAGAGACACTCCGATTCTTGACATCAAAATGTCTTTGATTGTAAATGTACATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATc                                                                                            <  2:2747233/141‑1 (MQ=255)
                                                     aTGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTc                                                                                        >  2:3805754/1‑148 (MQ=255)
                                                     aTGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTc                                                                                        <  2:1859353/148‑1 (MQ=255)
                                                           aagaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATg                                                                                  >  1:3394791/1‑148 (MQ=255)
                                                             gaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGaa                                                                                <  2:1941979/148‑1 (MQ=255)
                                                                taatGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTAATCTTTACCCATGAATtc                                                                             >  1:2024902/1‑148 (MQ=255)
                                                                  atGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATtcat                                                                           >  1:2305995/1‑148 (MQ=255)
                                                                    gTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATtcatca                                                                         <  2:3759211/148‑1 (MQ=255)
                                                                     tttGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTcatcat                                                                        >  2:96381/1‑148 (MQ=255)
                                                                     tttGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTcatcat                                                                        >  1:2843323/1‑148 (MQ=255)
                                                                      ttGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCAtt                                                                       <  1:523372/148‑1 (MQ=255)
                                                                      ttGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCAtt                                                                       >  1:2472144/1‑148 (MQ=255)
                                                                       tGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTa                                                                      <  2:1367520/148‑1 (MQ=255)
                                                                        gAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAg                                                                     >  1:1099953/1‑148 (MQ=255)
                                                                             tCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTc                                                                <  2:3633603/148‑1 (MQ=255)
                                                                              ccTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCa                                                               <  1:2100588/148‑1 (MQ=255)
                                                                               cTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAg                                                              >  1:3726226/1‑148 (MQ=255)
                                                                                tGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGt                                                             <  1:2785312/148‑1 (MQ=255)
                                                                                  agaCACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTg                                                           >  1:318963/1‑148 (MQ=255)
                                                                                   gaCACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGa                                                          >  2:819613/1‑148 (MQ=255)
                                                                                   gaCACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGa                                                          <  1:3235385/148‑1 (MQ=255)
                                                                                   gaCACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGa                                                          >  2:2378739/1‑148 (MQ=255)
                                                                                        tttGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTa                                                     <  2:3395155/148‑1 (MQ=255)
                                                                                          tGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAAt                                                   <  2:2535955/148‑1 (MQ=255)
                                                                                          tGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAAt                                                   <  2:1190972/148‑1 (MQ=255)
                                                                                           gATGCTGTACATCAACATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATc                                                  <  2:3497571/148‑1 (MQ=255)
                                                                                               cTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCa                                              <  2:2992365/148‑1 (MQ=255)
                                                                                               cTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCa                                              >  1:818878/1‑148 (MQ=255)
                                                                                               cTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCa                                              <  1:1013697/148‑1 (MQ=255)
                                                                                                tttACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCaa                                             <  2:3075349/148‑1 (MQ=255)
                                                                                                 ttACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAAt                                            >  1:1402199/1‑148 (MQ=255)
                                                                                                 ttACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAAt                                            >  1:1734673/1‑148 (MQ=255)
                                                                                                 ttACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAAt                                            >  1:25979/1‑148 (MQ=255)
                                                                                                  tACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATg                                           <  2:579212/148‑1 (MQ=255)
                                                                                                   aCATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATgg                                          <  2:2123314/148‑1 (MQ=255)
                                                                                                     aTCAACATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATgggg                                        <  2:1436205/148‑1 (MQ=255)
                                                                                                     aTCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATgggg                                        <  2:2258909/148‑1 (MQ=255)
                                                                                                       cAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGAt                                      >  1:1322895/1‑148 (MQ=255)
                                                                                                       cAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGAt                                      <  1:96381/148‑1 (MQ=255)
                                                                                                       cAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGAt                                      >  1:3130954/1‑148 (MQ=255)
                                                                                                         aaaTGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGAttt                                    >  2:909476/1‑148 (MQ=255)
                                                                                                            tGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTa                                 <  1:3633291/148‑1 (MQ=255)
                                                                                                            tGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTa                                 >  1:164989/1‑148 (MQ=255)
                                                                                                            tGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTa                                 >  1:572241/1‑148 (MQ=255)
                                                                                                             gTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTaa                                >  1:2693361/1‑148 (MQ=255)
                                                                                                             gTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGNAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTaa                                >  1:3663655/1‑148 (MQ=255)
                                                                                                              tCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAAt                               >  1:3108603/1‑148 (MQ=255)
                                                                                                                tttGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATgg                             <  2:282246/148‑1 (MQ=255)
                                                                                                                 ttGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGa                            >  1:3450701/1‑148 (MQ=255)
                                                                                                                  tGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGAc                           <  1:3292418/148‑1 (MQ=255)
                                                                                                                     ttGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCa                        <  1:86613/148‑1 (MQ=255)
                                                                                                                     ttGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCa                        >  2:1753296/1‑148 (MQ=255)
                                                                                                                      tgtaaatgtaaaTGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCat                       <  1:430302/148‑1 (MQ=255)
                                                                                                                      tgtaaatgtaaaTGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCat                       <  2:2073916/148‑1 (MQ=255)
                                                                                                                      tgtaaatgtaaaTGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCat                       <  1:256390/148‑1 (MQ=255)
                                                                                                                        taaatgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatga                     <  2:1051170/148‑1 (MQ=255)
                                                                                                                         aaatgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatgat                    >  1:1848391/1‑148 (MQ=255)
                                                                                                                          aatgtaactgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatgatg                   <  2:874674/148‑1 (MQ=255)
                                                                                                                           atgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTTTTTTTTCATCTTTCCCCCTGAATTCATCATTCGGGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGt                  <  1:350831/148‑1 (MQ=255)
                                                                                                                           atgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGt                  <  2:1266374/148‑1 (MQ=255)
                                                                                                                           atgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGt                  <  2:735473/148‑1 (MQ=255)
                                                                                                                           atgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGt                  <  1:3685314/148‑1 (MQ=255)
                                                                                                                           atgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTGTTAATGGACTCATGATGt                  >  2:2454603/1‑148 (MQ=255)
                                                                                                                            tgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGtt                 >  1:1504784/1‑148 (MQ=255)
                                                                                                                            tgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGtt                 >  2:2332613/1‑148 (MQ=255)
                                                                                                                            tgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGtt                 <  2:3725678/148‑1 (MQ=255)
                                                                                                                            tgtaaatgGCTTTTTACATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCAAGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTTCAATGGGGATTTTTAATGGACTCATGATGtt                 <  2:2917376/148‑1 (MQ=255)
                                                                                                                             gtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTc                >  1:3349002/1‑148 (MQ=255)
                                                                                                                               aaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGt              >  2:1728454/1‑148 (MQ=255)
                                                                                                                                    gCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTcc         <  2:414701/148‑1 (MQ=255)
                                                                                                                                     cTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCATTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCa        >  2:3322542/1‑148 (MQ=255)
                                                                                                                                      tttttCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAg       <  1:1910105/148‑1 (MQ=255)
                                                                                                                                       ttttCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGa      >  2:2895362/1‑148 (MQ=255)
                                                                                                                                        tttCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGAt     >  2:1782598/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAA  >  CP009273/289867‑290153

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: