Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,706,450 C→T 14.1% intergenic (+35/‑576) nth → / → dtpA DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,706,4500CT14.1% 791.4 / 93.3 270intergenic (+35/‑576)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base C (118/114);  new base T (16/22);  total (134/136)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 277 total aligned reads displayed.

CTGCCACCATTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATC  >  CP009273/1706301‑1706599
                                                                                                                                                     |                                                                                                                                                     
cTGCCACCATTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAAc                                                                                                                                                       <  2:4571817/150‑1 (MQ=255)
cTGCCACCATTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAAc                                                                                                                                                       >  1:1317098/1‑150 (MQ=255)
  gCCACCATTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACct                                                                                                                                                     >  2:4236201/1‑150 (MQ=255)
   ccaccaTTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGTCTCTTGTATTATTGTAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTTAAGAAAAGGGGTAACACCGATTACCCCATTGATAATctt                                                                                                                                                    <  1:3033178/150‑1 (MQ=255)
    caccaTTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACcttt                                                                                                                                                   <  1:2192420/150‑1 (MQ=255)
      ccaTTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACctttct                                                                                                                                                 >  1:2288632/1‑150 (MQ=255)
         ttGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTa                                                                                                                                              >  2:2286494/1‑150 (MQ=255)
         ttGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTa                                                                                                                                              >  1:1492964/1‑150 (MQ=255)
           ggTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATc                                                                                                                                            <  2:4311734/150‑1 (MQ=255)
            gTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATcc                                                                                                                                           <  1:1053446/150‑1 (MQ=255)
             ttGATCCTGCACGGGCGTTATACCTGCATTGCCCGCTAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGTGTAACACCGATTACCCCATTGATAATCTTTCTTTATCct                                                                                                                                          <  1:1815209/150‑1 (MQ=255)
               gATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCctct                                                                                                                                        <  2:1902599/150‑1 (MQ=255)
               gATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCctct                                                                                                                                        <  1:2031609/150‑1 (MQ=255)
                aTCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCtt                                                                                                                                       <  2:3709003/150‑1 (MQ=255)
                aTCCTGCACGGGCGTTATACCTGCATTGCACGCAAGCCCCTCTGTGTCACTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGATTACATCTGAAGAAAAGGGTTAACACCGATAACCCCATTGATAACCTTTCTTTATCCTCtt                                                                                                                                       >  1:2661092/1‑150 (MQ=255)
                 tCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCttt                                                                                                                                      <  1:1783526/150‑1 (MQ=255)
                    tGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTaa                                                                                                                                   >  2:2759877/1‑150 (MQ=255)
                     gCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTaaa                                                                                                                                  >  1:1547057/1‑150 (MQ=255)
                      cACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTaaaa                                                                                                                                 <  1:3633619/150‑1 (MQ=255)
                       aCGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAAc                                                                                                                                <  1:2361715/150‑1 (MQ=255)
                            cGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCtt                                                                                                                           <  2:2063472/150‑1 (MQ=255)
                                ataCCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTa                                                                                                                       >  2:1639691/1‑150 (MQ=255)
                                 taCCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTaa                                                                                                                      >  2:3602165/1‑150 (MQ=255)
                                  aCCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAAc                                                                                                                     >  1:3841114/1‑150 (MQ=255)
                                   ccTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACg                                                                                                                    <  2:676764/150‑1 (MQ=255)
                                        aTTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAat                                                                                                               <  2:953455/150‑1 (MQ=255)
                                         ttGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAatg                                                                                                              <  2:3333350/150‑1 (MQ=255)
                                         ttGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAatg                                                                                                              >  1:563250/1‑150 (MQ=255)
                                           gCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAatgat                                                                                                            <  1:1959840/150‑1 (MQ=255)
                                           gCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAatgat                                                                                                            <  1:1959780/150‑1 (MQ=255)
                                            cccGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAatgatg                                                                                                           >  1:591669/1‑150 (MQ=255)
                                            cccGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAatgatg                                                                                                           <  2:3872008/150‑1 (MQ=255)
                                               gCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCa                                                                                                        <  1:3291108/150‑1 (MQ=255)
                                                 aaGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAACAGTGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCAtt                                                                                                      <  2:887625/150‑1 (MQ=255)
                                                 aaGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCAtt                                                                                                      >  1:2822408/1‑150 (MQ=255)
                                                 aaGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCAtt                                                                                                      >  2:2224167/1‑150 (MQ=255)
                                                 aaGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCAtt                                                                                                      <  1:2094266/150‑1 (MQ=255)
                                                 aaGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCAtt                                                                                                      <  2:1590091/150‑1 (MQ=255)
                                                    ccccGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCt                                                                                                   >  2:1495515/1‑150 (MQ=255)
                                                     cccGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCtt                                                                                                  <  2:2450194/150‑1 (MQ=255)
                                                        gCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCATTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGc                                                                                               >  1:3406040/1‑150 (MQ=255)
                                                          tgtgGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGt                                                                                             <  2:1746295/150‑1 (MQ=255)
                                                           gtgGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGtt                                                                                            <  2:605120/150‑1 (MQ=255)
                                                              gCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTatc                                                                                         >  1:4428975/1‑150 (MQ=255)
                                                                  ttGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTAtcatca                                                                                     <  2:2933571/150‑1 (MQ=255)
                                                                     tattatTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTGAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGt                                                                                  >  2:3065444/1‑150 (MQ=255)
                                                                       ttattGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTaa                                                                                <  1:2444251/150‑1 (MQ=255)
                                                                         attGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATc                                                                              >  1:4587044/1‑150 (MQ=255)
                                                                          ttGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTAACCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATcc                                                                             <  1:3570657/150‑1 (MQ=255)
                                                                            gAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGt                                                                           <  2:3984729/150‑1 (MQ=255)
                                                                                  cTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCAGTGCTTAGCGTTATCATCAGGTAATCCGTTTGCgg                                                                     >  1:2850629/1‑150 (MQ=255)
                                                                                      gtgAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATaa                                                                 >  1:528004/1‑150 (MQ=255)
                                                                                         aaTACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCa                                                              >  1:1607695/1‑150 (MQ=255)
                                                                                         aaTACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCa                                                              >  1:1145820/1‑150 (MQ=255)
                                                                                          aTACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCaa                                                             <  1:156900/150‑1 (MQ=255)
                                                                                           tACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCaaa                                                            <  2:3243240/150‑1 (MQ=255)
                                                                                            aCAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTAGCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCAGTTTGATGATAACCaaaa                                                           >  1:494661/1‑150 (MQ=255)
                                                                                            aCAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCGTTCTTTATCCACTTTTAAAACATCTATTTAACGACAATGATGCCAATGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACAaaaa                                                           >  2:2955083/1‑150 (MQ=255)
                                                                                              aaaGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAAt                                                         <  1:708059/150‑1 (MQ=255)
                                                                                               aaGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATg                                                        <  1:4239304/150‑1 (MQ=255)
                                                                                                 gagaAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCa                                                      <  1:3740586/150‑1 (MQ=255)
                                                                                                    aaaGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGt                                                   <  2:2777281/150‑1 (MQ=255)
                                                                                                     aaGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGtt                                                  <  2:2152019/150‑1 (MQ=255)
                                                                                                     aaGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGtt                                                  >  2:3420418/1‑150 (MQ=255)
                                                                                                       ggtGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTaa                                                >  2:1023891/3‑150 (MQ=255)
                                                                                                        ttGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAAt                                               <  2:77639/150‑1 (MQ=255)
                                                                                                        ttGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAAt                                               <  1:472940/150‑1 (MQ=255)
                                                                                                            cATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtt                                           >  1:2611517/1‑150 (MQ=255)
                                                                                                            cATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtt                                           >  2:4145861/1‑150 (MQ=255)
                                                                                                            cATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtt                                           >  1:3431859/1‑150 (MQ=255)
                                                                                                               cTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttt                                        <  1:1239706/150‑1 (MQ=255)
                                                                                                                tGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttt                                       <  2:2432342/150‑1 (MQ=255)
                                                                                                                tGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttt                                       <  2:400315/150‑1 (MQ=255)
                                                                                                                tGAAGAAAAGGGGTAACACCGATTACCCCAATGATAACCTTTCTTTATCCTCTTTTAAAACATCTTATTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtatttt                                       <  2:2432574/150‑1 (MQ=255)
                                                                                                                  aagaaAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTa                                     <  1:678963/150‑1 (MQ=255)
                                                                                                                       aaGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGc                                >  2:4076781/1‑150 (MQ=255)
                                                                                                                         ggggTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGa                              >  2:2735532/1‑150 (MQ=255)
                                                                                                                             tAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAAtt                          <  2:496064/150‑1 (MQ=255)
                                                                                                                             tAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAAtt                          <  2:2810194/150‑1 (MQ=255)
                                                                                                                               acacCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAAtttt                        <  2:4391814/150‑1 (MQ=255)
                                                                                                                               acacCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAAtttt                        <  1:4353424/150‑1 (MQ=255)
                                                                                                                                cacCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTc                       >  1:3255246/1‑150 (MQ=255)
                                                                                                                                        aCCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCATGTAATCCGTTTGCTGATAACCAAAAATGCAGGTTAATTTTTTTTTTAATAGCGAAATTTTCTattcatt               >  1:1117438/1‑150 (MQ=255)
                                                                                                                                        aCCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTattcatt               >  2:2573683/1‑150 (MQ=255)
                                                                                                                                         ccccATTGATAATCTTTCTTTATCCTCTTTTAAACCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTattcattc              <  1:2142999/150‑1 (MQ=255)
                                                                                                                                         ccccATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTattcattc              >  2:1993304/1‑150 (MQ=255)
                                                                                                                                         ccccATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTattcattc              >  1:4432462/1‑150 (MQ=255)
                                                                                                                                         ccccATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTattcattc              >  2:3860552/1‑150 (MQ=255)
                                                                                                                                           ccATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGt            <  2:3927176/150‑1 (MQ=255)
                                                                                                                                           ccATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGt            <  1:3247448/150‑1 (MQ=255)
                                                                                                                                               tGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATc        >  1:3133775/1‑150 (MQ=255)
                                                                                                                                               tGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATc        >  1:1380516/1‑150 (MQ=255)
                                                                                                                                                gATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCa       >  2:3841632/1‑150 (MQ=255)
                                                                                                                                                 aTAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCaa      >  2:3816009/1‑150 (MQ=255)
                                                                                                                                                  tAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAg     <  1:1085214/150‑1 (MQ=255)
                                                                                                                                                  tAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAg     >  1:83956/1‑150 (MQ=255)
                                                                                                                                                  tAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAg     >  1:4199918/1‑150 (MQ=255)
                                                                                                                                                    aCCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt   >  1:2688368/1‑150 (MQ=255)
                                                                                                                                                     ccTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATc  >  1:2894028/1‑150 (MQ=255)
                                                                                                                                                     ccTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATc  >  1:3343749/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCCACCATTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATC  >  CP009273/1706301‑1706599

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: