Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,899,442 C→T 10.5% R129C (CGT→TGT)  yobD → UPF0266 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,899,4420CT10.5% 449.3 / 32.9 143R129C (CGT→TGT) yobDUPF0266 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (66/62);  new base T (3/12);  total (69/74)
Fisher's exact test for biased strand distribution p-value = 2.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 147 total aligned reads displayed.

TGATGGGTTTTTATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCTATATTACCTCTTCC  >  CP009273/1899293‑1899591
                                                                                                                                                     |                                                                                                                                                     
tGATGGGTTTTTATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGc                                                                                                                                                       >  2:975996/1‑150 (MQ=255)
tGATGGGTTTTTATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGc                                                                                                                                                       <  1:1637012/150‑1 (MQ=255)
 gATGGGTTTTTATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCg                                                                                                                                                      <  1:29345/150‑1 (MQ=255)
   tGGGTTTTTATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGTGTc                                                                                                                                                    <  2:3250465/150‑1 (MQ=255)
   tGGGTTTTTATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTc                                                                                                                                                    <  2:1550531/150‑1 (MQ=255)
   tGGGTTTTTATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTc                                                                                                                                                    <  2:4288695/150‑1 (MQ=255)
   tGGGTTTTTATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTc                                                                                                                                                    >  1:1939976/1‑150 (MQ=255)
    gggTTTTTATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCg                                                                                                                                                   >  1:2587032/1‑150 (MQ=255)
       tttttATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCt                                                                                                                                                <  2:1788093/150‑1 (MQ=255)
          ttATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGtt                                                                                                                                             <  2:2383507/150‑1 (MQ=255)
           tatataTTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGTGTCGGCTGTTa                                                                                                                                            <  2:3276604/150‑1 (MQ=255)
           tatataTTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGTGTCGGCTGTTa                                                                                                                                            <  2:3276135/150‑1 (MQ=255)
             tatatTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAAt                                                                                                                                          >  1:2811833/1‑150 (MQ=255)
              atatTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATc                                                                                                                                         >  2:910291/1‑150 (MQ=255)
                  tCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCcgcg                                                                                                                                     >  2:2993901/1‑150 (MQ=255)
                  tCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCcgcg                                                                                                                                     >  2:3241286/1‑150 (MQ=255)
                   cTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGt                                                                                                                                    >  2:1932116/1‑150 (MQ=255)
                     ggATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTc                                                                                                                                  >  2:3190920/1‑150 (MQ=255)
                      gATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCg                                                                                                                                 <  2:4156140/150‑1 (MQ=255)
                      gATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCg                                                                                                                                 <  2:4156171/150‑1 (MQ=255)
                       aTCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGa                                                                                                                                <  2:1326858/150‑1 (MQ=255)
                             gTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATc                                                                                                                          <  1:899571/150‑1 (MQ=255)
                             gTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATc                                                                                                                          >  2:872743/1‑150 (MQ=255)
                             gTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATc                                                                                                                          >  1:1632728/1‑150 (MQ=255)
                             gTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATc                                                                                                                          >  2:872646/1‑150 (MQ=255)
                              ttCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATcg                                                                                                                         <  2:2444468/150‑1 (MQ=255)
                                ccGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATcgac                                                                                                                       <  1:1720045/150‑1 (MQ=255)
                                   aaGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGAt                                                                                                                    >  1:1685666/1‑150 (MQ=255)
                                   aaGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGAt                                                                                                                    >  1:1685878/1‑150 (MQ=255)
                                      atcatcTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGTGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTg                                                                                                                 >  2:717910/1‑150 (MQ=255)
                                      atcatcTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTg                                                                                                                 <  1:2570593/150‑1 (MQ=255)
                                      atcatcTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTg                                                                                                                 >  1:2332385/1‑150 (MQ=255)
                                      atcatcTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTg                                                                                                                 >  1:2329206/1‑150 (MQ=255)
                                       tcatcTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTgg                                                                                                                >  2:1725813/1‑150 (MQ=255)
                                        catcTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGa                                                                                                               >  1:343834/1‑150 (MQ=255)
                                        catcTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGa                                                                                                               >  1:4302297/1‑150 (MQ=255)
                                         atcTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGaa                                                                                                              <  1:1739966/150‑1 (MQ=255)
                                           cTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGaaaa                                                                                                            <  2:4744966/150‑1 (MQ=255)
                                           cTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATAAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGTGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGaaaa                                                                                                            <  1:553100/150‑1 (MQ=255)
                                             ttAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGTGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGaaaaaa                                                                                                          <  1:411745/150‑1 (MQ=255)
                                             ttAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGTGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGaaaaaa                                                                                                          <  1:411877/150‑1 (MQ=255)
                                               aaaCAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGTGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAAtt                                                                                                        <  2:3377745/150‑1 (MQ=255)
                                                aaCAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAAttt                                                                                                       <  2:1819635/150‑1 (MQ=255)
                                                 aCAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTta                                                                                                      <  1:4185738/150‑1 (MQ=255)
                                                   aaaaaGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGTGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTtata                                                                                                    <  2:1017545/150‑1 (MQ=255)
                                                   aaaaaGGTTTGTTCTTCGCCAATGTCGGGATTGAATATAGCCGAGTCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTtata                                                                                                    <  2:2819250/150‑1 (MQ=255)
                                                    aaaaGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATaa                                                                                                   <  1:492028/150‑1 (MQ=255)
                                                     aaaGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAg                                                                                                  <  2:3461509/150‑1 (MQ=255)
                                                     aaaGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAg                                                                                                  >  1:3384648/1‑150 (MQ=255)
                                                     aaaGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAg                                                                                                  >  1:2643018/1‑150 (MQ=255)
                                                      aaGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAgt                                                                                                 <  2:4298180/150‑2 (MQ=255)
                                                      aaGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGc                                                                                                 >  1:2829128/1‑150 (MQ=255)
                                                       aGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCt                                                                                                >  1:2392962/1‑150 (MQ=255)
                                                           tttttCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTtctc                                                                                            <  2:3167799/150‑1 (MQ=255)
                                                             tttCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGt                                                                                          >  1:695123/1‑150 (MQ=255)
                                                              ttcttcGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGtt                                                                                         >  1:3933197/1‑150 (MQ=255)
                                                              ttcttcGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGtt                                                                                         >  1:1543965/1‑150 (MQ=255)
                                                               tcttcGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGttt                                                                                        >  1:2478638/1‑150 (MQ=255)
                                                                  tcGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCaa                                                                                     >  2:2095919/1‑150 (MQ=255)
                                                                       aaTGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCaa                                                                                >  1:1962270/1‑150 (MQ=255)
                                                                        aTGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAAt                                                                               >  1:2692944/1‑150 (MQ=255)
                                                                           tCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAg                                                                            <  2:1056307/150‑1 (MQ=255)
                                                                             tGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCACTTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGtt                                                                          >  1:1064833/1‑150 (MQ=255)
                                                                              ggATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTa                                                                         >  2:182295/1‑150 (MQ=255)
                                                                               gATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTAt                                                                        >  1:3062882/1‑150 (MQ=255)
                                                                               gATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTAt                                                                        >  1:346952/1‑150 (MQ=255)
                                                                               gATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTAt                                                                        >  1:346916/1‑150 (MQ=255)
                                                                                aTTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATg                                                                       >  2:3503752/1‑150 (MQ=255)
                                                                                      tataGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTa                                                                 <  2:4648867/150‑1 (MQ=255)
                                                                                        taGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGc                                                               <  2:1318972/150‑1 (MQ=255)
                                                                                          gCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAAGATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCa                                                             <  2:1544278/150‑1 (MQ=255)
                                                                                                  aaaGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATg                                                     >  2:1150406/1‑150 (MQ=255)
                                                                                                     gCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGttt                                                  >  2:4469144/1‑150 (MQ=255)
                                                                                                     gCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGttt                                                  >  2:4469073/1‑150 (MQ=255)
                                                                                                     gCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGttt                                                  >  2:3212526/1‑150 (MQ=255)
                                                                                                       gATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAg                                                <  1:1945207/150‑1 (MQ=255)
                                                                                                            aCTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTAt                                           >  2:3205227/1‑150 (MQ=255)
                                                                                                               tGTCGGAAGATGGCGTGTTGGTGATGCAATTCGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTAtttt                                        <  2:3857731/150‑1 (MQ=255)
                                                                                                                       gATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATc                                <  1:1060600/150‑1 (MQ=255)
                                                                                                                       gATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATc                                <  1:4001978/150‑1 (MQ=255)
                                                                                                                        aTGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCa                               <  1:342112/150‑1 (MQ=255)
                                                                                                                         tGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAg                              >  2:3678591/1‑150 (MQ=255)
                                                                                                                          ggCGTGTTGGTGATGCAATTAGAACAGTGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGa                             <  1:2517292/150‑1 (MQ=255)
                                                                                                                          ggCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGa                             <  1:2401421/150‑1 (MQ=255)
                                                                                                                            cGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACa                           <  2:3013544/150‑1 (MQ=255)
                                                                                                                            cGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACa                           <  2:3013537/150‑1 (MQ=255)
                                                                                                                            cGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACa                           <  2:1808183/150‑1 (MQ=255)
                                                                                                                              tgtTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATa                         <  2:3681710/150‑1 (MQ=255)
                                                                                                                               gtTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAg                        <  1:1725813/150‑1 (MQ=255)
                                                                                                                                       tGCAATTAGAACAGTGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCt                <  1:1767668/150‑1 (MQ=255)
                                                                                                                                       tGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCt                <  1:3279593/150‑1 (MQ=255)
                                                                                                                                         cAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCtat              >  1:1088684/1‑150 (MQ=255)
                                                                                                                                         cAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCtat              <  2:1139820/150‑1 (MQ=255)
                                                                                                                                         cAATTAGAACAGCGTCGGCTGTTAATCCGCGTTAGAAATATCGACGATCTGGACAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCtat              <  2:724866/150‑1 (MQ=255)
                                                                                                                                          aaTTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCtata             <  1:3190934/150‑1 (MQ=255)
                                                                                                                                            ttAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCTATAtt           >  2:2176241/1‑150 (MQ=255)
                                                                                                                                            ttAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCTATAtt           <  1:1240680/150‑1 (MQ=255)
                                                                                                                                            ttAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCTATAtt           <  2:2643018/150‑1 (MQ=255)
                                                                                                                                                   aGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCTATATTACCTCtt    >  2:1122080/1‑150 (MQ=255)
                                                                                                                                                     cGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCTATATTACCTCTTcc  <  2:3384648/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATGGGTTTTTATATATTCTGGATCCGCGTTCCGAAGATCATCTTTAAACAAAAAGGTTTTTTCTTCGCCAATGTCTGGATTGAATATAGCCGAATCAAAGCGATGAACTTGTCGGAAGATGGCGTGTTGGTGATGCAATTAGAACAGCGTCGGCTGTTAATCCGCGTTCGAAATATCGACGATCTGGAAAAAATTTATAAGCTTCTCGTTTCAACTCAATAAGTTATGAATTTAGCCAAAGCTATGTTTAGTGTATTTTTAATAATCAGACATAGCTTAGGCTATATTACCTCTTCC  >  CP009273/1899293‑1899591

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: