Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,706,450 C→T 14.2% 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.2% 850.9 / 100.7 289intergenic (+35/‑576)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base C (109/139);  new base T (21/20);  total (130/159)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 291 total aligned reads displayed.

CTGCCACCATTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATC  >  CP009273/1706301‑1706599
                                                                                                                                                     |                                                                                                                                                     
   ccaccaTTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACctt                                                                                                                                                    <  2:1226361/150‑1 (MQ=255)
    caccaTTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACcttt                                                                                                                                                   <  2:3963051/150‑1 (MQ=255)
     accaTTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACctttc                                                                                                                                                  <  1:40456/150‑1 (MQ=255)
      ccaTTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATAAAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACctttct                                                                                                                                                 >  2:4021476/1‑150 (MQ=255)
       caTTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACctttctt                                                                                                                                                <  2:3908377/150‑1 (MQ=255)
        aTTGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACctttcttt                                                                                                                                               <  2:3931524/150‑1 (MQ=255)
          tGGTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTAt                                                                                                                                             <  2:894346/150‑1 (MQ=255)
            gTTGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATcc                                                                                                                                           >  2:1412411/1‑150 (MQ=255)
              tGATCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCctc                                                                                                                                         <  2:151447/150‑1 (MQ=255)
                aTCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCtt                                                                                                                                       <  2:1696420/150‑1 (MQ=255)
                aTCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCtt                                                                                                                                       <  1:2663625/150‑1 (MQ=255)
                aTCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCtt                                                                                                                                       <  2:738102/150‑1 (MQ=255)
                 tCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCttt                                                                                                                                      >  2:3447317/1‑150 (MQ=255)
                 tCCTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCttt                                                                                                                                      <  2:1049343/150‑1 (MQ=255)
                 tCCTGCACGGGCGTTATAACTGCTTTGCGCGCTAGCCACGCTGTGGGTATTGTTTTCTTGAAGTTCTTTGTGAATACTAAGAGAATGTTGACATCTGAAGTAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCttt                                                                                                                                      <  1:4597902/150‑1 (MQ=255)
                   cTGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGGATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTa                                                                                                                                    >  2:544530/1‑150 (MQ=255)
                    tGCACGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTaa                                                                                                                                   <  2:1854086/150‑1 (MQ=255)
                       aCGGGCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAAc                                                                                                                                <  1:1216137/150‑1 (MQ=255)
                         gggCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACAt                                                                                                                              <  2:1323517/150‑1 (MQ=255)
                         gggCGTTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACAt                                                                                                                              <  2:3980658/150‑1 (MQ=255)
                               tataCCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCttttt                                                                                                                        >  1:3680478/1‑150 (MQ=255)
                                ataCCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTa                                                                                                                       <  1:1352457/150‑1 (MQ=255)
                                  aCCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAAc                                                                                                                     <  1:3851447/150‑1 (MQ=255)
                                    cTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGt                                                                                                                   <  1:69536/150‑1 (MQ=255)
                                      gCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCa                                                                                                                 <  1:2267001/150‑1 (MQ=255)
                                        aTTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAat                                                                                                               <  1:4584603/150‑1 (MQ=255)
                                         ttGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAatg                                                                                                              <  1:3333809/150‑1 (MQ=255)
                                         ttGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAatg                                                                                                              <  1:3220541/150‑1 (MQ=255)
                                          tGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAatga                                                                                                             >  1:2757733/1‑150 (MQ=255)
                                            cccGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAatgatg                                                                                                           <  1:2452102/150‑1 (MQ=255)
                                                 aaGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCAtt                                                                                                      <  1:2564661/150‑1 (MQ=255)
                                                   gCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGc                                                                                                    <  2:4400755/150‑1 (MQ=255)
                                                    ccccGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCt                                                                                                   <  1:822118/150‑1 (MQ=255)
                                                     cccGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCtt                                                                                                  <  2:1487639/150‑1 (MQ=255)
                                                      ccGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACTAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTa                                                                                                 <  1:3001963/150‑1 (MQ=255)
                                                       cGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGc                                                                                               >  2:675387/1‑150 (MQ=255)
                                                        gCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGc                                                                                               >  2:130354/1‑150 (MQ=255)
                                                           gtgGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGtt                                                                                            >  1:259126/1‑150 (MQ=255)
                                                           gtgGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGtt                                                                                            >  2:2707173/1‑150 (MQ=255)
                                                              gCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTatc                                                                                         <  2:3248471/150‑1 (MQ=255)
                                                              gCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTatc                                                                                         >  1:18596/1‑150 (MQ=255)
                                                               ctctTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTatca                                                                                        <  1:1583691/150‑1 (MQ=255)
                                                                 ctTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTatcatc                                                                                      <  1:152919/150‑1 (MQ=255)
                                                                   tGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAg                                                                                    >  1:3531212/1‑150 (MQ=255)
                                                                     tattatTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGt                                                                                  >  2:2952695/1‑150 (MQ=255)
                                                                      attattGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTa                                                                                 >  2:2908963/1‑150 (MQ=255)
                                                                           tGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCg                                                                            <  1:303154/150‑1 (MQ=255)
                                                                           tGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCg                                                                            >  2:1532908/1‑150 (MQ=255)
                                                                            gAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGt                                                                           >  1:2945041/1‑150 (MQ=255)
                                                                            gAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGt                                                                           >  1:930127/1‑150 (MQ=255)
                                                                            gAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGt                                                                           >  2:3377255/1‑150 (MQ=255)
                                                                              aGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGttt                                                                         >  1:1750349/1‑150 (MQ=255)
                                                                               gATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTg                                                                        >  2:1833703/1‑150 (MQ=255)
                                                                                aTCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGc                                                                       <  1:1617629/150‑1 (MQ=255)
                                                                                     tgtgAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATa                                                                  >  2:2373872/1‑150 (MQ=255)
                                                                                     tgtgAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATa                                                                  <  2:403031/150‑1 (MQ=255)
                                                                                        gAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAAcc                                                               <  1:1978463/150‑1 (MQ=255)
                                                                                         aaTACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCa                                                              <  2:1380041/150‑1 (MQ=255)
                                                                                         aaTACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCa                                                              <  1:604474/150‑1 (MQ=255)
                                                                                             cAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCaaaaa                                                          >  1:4151629/1‑150 (MQ=255)
                                                                                              aaaGAGAAAGTTGACATCTGAAGAAAAGGGGTAATACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAAt                                                         <  1:3133145/150‑1 (MQ=255)
                                                                                              aaaGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAAt                                                         <  1:2815132/150‑1 (MQ=255)
                                                                                              aaaGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAAt                                                         <  1:1008290/150‑1 (MQ=255)
                                                                                               aaGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATg                                                        >  1:3905652/1‑150 (MQ=255)
                                                                                               aaGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATg                                                        >  1:2406021/1‑150 (MQ=255)
                                                                                                  agaAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAg                                                     >  2:4584082/1‑150 (MQ=255)
                                                                                                   gaAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAgg                                                    <  2:1122519/150‑1 (MQ=255)
                                                                                                   gaAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAgg                                                    >  1:2224582/1‑150 (MQ=255)
                                                                                                   gaAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAgg                                                    >  1:1412991/1‑150 (MQ=255)
                                                                                                    aaaGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGt                                                   <  2:3724902/150‑1 (MQ=255)
                                                                                                     aaGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGtt                                                  >  2:3140794/1‑150 (MQ=255)
                                                                                                           aCATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGt                                            <  2:693210/150‑1 (MQ=255)
                                                                                                             aTCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCATCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttt                                          >  1:2243461/1‑150 (MQ=255)
                                                                                                              tCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttt                                         >  1:1869121/1‑150 (MQ=255)
                                                                                                                tGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttt                                       <  2:929487/150‑1 (MQ=255)
                                                                                                                tGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttt                                       <  1:1126378/150‑1 (MQ=255)
                                                                                                                 gaagaaAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttttt                                      <  2:1628728/150‑1 (MQ=255)
                                                                                                                     aaaaGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATa                                  >  2:3595433/1‑150 (MQ=255)
                                                                                                                     aaaaGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATa                                  <  1:2051344/150‑1 (MQ=255)
                                                                                                                     aaaaGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATa                                  <  2:2864361/150‑1 (MQ=255)
                                                                                                                       aaGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGc                                >  1:759708/1‑150 (MQ=255)
                                                                                                                        aGGGGTAACACCGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCg                               >  2:3024196/1‑150 (MQ=255)
                                                                                                                            gTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAAt                           >  1:1316499/1‑150 (MQ=255)
                                                                                                                                 acCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCt                      >  2:3509435/1‑150 (MQ=255)
                                                                                                                                 acCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCt                      <  2:2207332/150‑1 (MQ=255)
                                                                                                                                   cGATTACCCCATTGATAATCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTat                    <  2:4520797/150‑1 (MQ=255)
                                                                                                                                   cGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTat                    >  1:1087924/1‑150 (MQ=255)
                                                                                                                                     aTTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTattc                  <  2:3062426/150‑1 (MQ=255)
                                                                                                                                         ccacATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCTGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTattcattc              <  1:4582202/147‑1 (MQ=255)
                                                                                                                                          cccATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCg             >  1:436677/1‑150 (MQ=255)
                                                                                                                                          cccATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCg             <  2:3520229/150‑1 (MQ=255)
                                                                                                                                           ccATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGt            <  2:4087882/150‑1 (MQ=255)
                                                                                                                                           ccATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGt            >  2:2971657/1‑150 (MQ=255)
                                                                                                                                            cATTGATAACCTTACTTTATCCTCTTTTAAAACATCATTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCAGTTTGCGGATAACCAAAAATGCAGTTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTg           <  2:4114716/150‑1 (MQ=255)
                                                                                                                                             aTTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGa          >  2:2531907/1‑150 (MQ=255)
                                                                                                                                              ttGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGAt         >  2:492645/1‑150 (MQ=255)
                                                                                                                                                gATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCa       >  1:4445239/1‑150 (MQ=255)
                                                                                                                                                 aTAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCaa      >  2:3215385/1‑150 (MQ=255)
                                                                                                                                                   aaCCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGa    >  1:2306322/1‑150 (MQ=255)
                                                                                                                                                    aCCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt   >  2:2670958/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: