Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,931,041 A→T 7.9% N45Y (AAT→TAT)  yebK → putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,931,0410AT 7.9% 625.8 / 35.5 203N45Y (AAT→TAT) yebKputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (88/99);  new base T (6/10);  total (94/109)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 204 total aligned reads displayed.

GAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCA  >  CP009273/1930894‑1931187
                                                                                                                                                   |                                                                                                                                                  
gAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCa                                                                                                                                                    >  1:2925813/1‑148 (MQ=255)
   tctcTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATg                                                                                                                                                 <  1:1308931/148‑1 (MQ=255)
    ctctTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGt                                                                                                                                                >  2:2530526/1‑148 (MQ=255)
    ctctTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGt                                                                                                                                                >  1:1172748/1‑148 (MQ=255)
    ctctTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGt                                                                                                                                                >  1:2276405/1‑148 (MQ=255)
      ctTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTa                                                                                                                                              >  1:127152/1‑148 (MQ=255)
        ttACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGc                                                                                                                                            >  2:3268919/1‑148 (MQ=255)
          aCATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGa                                                                                                                                          >  1:2930875/1‑148 (MQ=255)
            atcatGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAAc                                                                                                                                        <  2:1569101/148‑1 (MQ=255)
                 gAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACg                                                                                                                                   <  1:1659393/148‑1 (MQ=255)
                 gAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACg                                                                                                                                   >  1:745571/1‑148 (MQ=255)
                 gAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACg                                                                                                                                   <  1:1659414/148‑1 (MQ=255)
                  aaTATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACgg                                                                                                                                  <  2:751464/148‑1 (MQ=255)
                  aaTATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACgg                                                                                                                                  <  2:2396249/148‑1 (MQ=255)
                         tGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCg                                                                                                                           >  2:2047887/1‑148 (MQ=255)
                          ggAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGt                                                                                                                          <  1:2379375/148‑1 (MQ=255)
                           gAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGtt                                                                                                                         >  1:1895667/1‑148 (MQ=255)
                           gAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGtt                                                                                                                         >  2:3311024/1‑148 (MQ=255)
                            aaaaaaTCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGttt                                                                                                                        <  1:1937838/148‑1 (MQ=255)
                               aaaTCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTg                                                                                                                     >  1:1051262/1‑148 (MQ=255)
                                  tCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCg                                                                                                                  <  2:2466820/148‑1 (MQ=255)
                                      gTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgcagc                                                                                                              <  1:3038298/148‑1 (MQ=255)
                                       tctcAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgcagca                                                                                                             <  1:70710/148‑1 (MQ=255)
                                       tctcAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgcagca                                                                                                             <  1:2437303/148‑1 (MQ=255)
                                       tctcAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgcagca                                                                                                             >  1:1019032/1‑148 (MQ=255)
                                       tctcAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgcagca                                                                                                             >  2:2283474/1‑148 (MQ=255)
                                        ctcAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCAt                                                                                                            <  1:379056/148‑1 (MQ=255)
                                        atcATCTGTAAAATTTTAGCAAATCAGATCGCAATGTTGCCGAGGTCTTTGTGGATTCGCCGTATAACGCGATCAATTAGAGTTTTGCTGATTTGGCACTGGATGCCATTGTTATCTAACCGACGGTGAATCGTTTCTGTCGCAGCAt                                                                                                            <  1:378917/147‑1 (MQ=255)
                                          cAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATgg                                                                                                          <  2:2313250/148‑1 (MQ=255)
                                            gCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGac                                                                                                        <  1:1217980/148‑1 (MQ=255)
                                            gCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGac                                                                                                        >  2:744628/1‑148 (MQ=255)
                                            gCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGac                                                                                                        >  1:2544754/1‑148 (MQ=255)
                                             cTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGaca                                                                                                       >  2:1851022/1‑148 (MQ=255)
                                               tgAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcg                                                                                                     >  1:3025608/2‑148 (MQ=255)
                                                gAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgc                                                                                                    >  2:674871/1‑148 (MQ=255)
                                                     tttGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGt                                                                                               >  1:1306210/1‑148 (MQ=255)
                                                     tttGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGt                                                                                               <  2:1225645/148‑1 (MQ=255)
                                                      ttGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGtt                                                                                              >  1:1511178/1‑148 (MQ=255)
                                                        gAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGtttt                                                                                            <  2:2644475/148‑1 (MQ=255)
                                                          gCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCTATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTcc                                                                                          >  1:3034841/1‑148 (MQ=255)
                                                              aTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTg                                                                                                                     <  1:2030986/117‑1 (MQ=255)
                                                               tCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGAtt                                                                                     >  2:316767/1‑148 (MQ=255)
                                                                cAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCACTGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACCCGCGCGGTTTTCCTGAttt                                                                                    >  1:2631846/1‑148 (MQ=255)
                                                                 agagCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGAtttt                                                                                   >  2:2210732/1‑148 (MQ=255)
                                                                   agCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATATCGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCTATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTACTGATTTTaa                                                                                 <  1:2978049/148‑1 (MQ=255)
                                                                     cgcAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAAc                                                                               >  2:1087157/1‑148 (MQ=255)
                                                                       cAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACtt                                                                             <  2:1581435/148‑1 (MQ=255)
                                                                          aGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCAt                                                                          <  1:519778/148‑1 (MQ=255)
                                                                          aGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCAt                                                                          <  1:2182357/148‑1 (MQ=255)
                                                                          aGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCAt                                                                          >  2:1443482/1‑148 (MQ=255)
                                                                            ttGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCt                                                                        <  2:305353/148‑1 (MQ=255)
                                                                               ccGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGc                                                                     >  2:1321468/1‑148 (MQ=255)
                                                                                  aGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGcaca                                                                  >  1:3176225/1‑148 (MQ=255)
                                                                                  aGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGcaca                                                                  >  1:1614114/1‑148 (MQ=255)
                                                                                   ggTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACag                                                                 >  1:1177778/1‑148 (MQ=255)
                                                                                   ggTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACag                                                                 <  1:1654961/148‑1 (MQ=255)
                                                                                    gTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACaga                                                                <  2:127152/148‑1 (MQ=255)
                                                                                    gTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACaga                                                                >  1:111650/1‑148 (MQ=255)
                                                                                     tCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACagag                                                               >  2:2742566/1‑148 (MQ=255)
                                                                                       aTTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTc                                                             <  2:216506/148‑1 (MQ=255)
                                                                                         tCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTg                                                           <  1:2047887/148‑1 (MQ=255)
                                                                                             gCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCTATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGa                                                       <  1:1395421/148‑1 (MQ=255)
                                                                                             gCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCTATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGa                                                       <  2:3034841/148‑1 (MQ=255)
                                                                                             gCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGa                                                       >  2:1960068/1‑148 (MQ=255)
                                                                                               ttCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAAt                                                     <  1:1490690/148‑1 (MQ=255)
                                                                                                tCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCTAATg                                                    >  2:1280765/1‑148 (MQ=255)
                                                                                                tCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATg                                                    <  2:294536/148‑1 (MQ=255)
                                                                                                   cccGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCa                                                 <  1:1227337/148‑1 (MQ=255)
                                                                                                   cccGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCa                                                 <  2:2001913/148‑1 (MQ=255)
                                                                                                     cGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACt                                               >  1:2714406/1‑148 (MQ=255)
                                                                                                       aTAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTcc                                             <  2:223817/148‑1 (MQ=255)
                                                                                                       aTAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTcc                                             >  1:1803019/1‑148 (MQ=255)
                                                                                                        tAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTccc                                            >  1:1891703/1‑148 (MQ=255)
                                                                                                          aCGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTa                                          <  2:2034889/148‑1 (MQ=255)
                                                                                                          aCGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTa                                          <  2:3122656/148‑1 (MQ=255)
                                                                                                          aCGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTa                                          <  2:2930776/148‑1 (MQ=255)
                                                                                                          aCGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTa                                          >  2:2881603/1‑148 (MQ=255)
                                                                                                          aCGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTa                                          >  1:2361471/1‑148 (MQ=255)
                                                                                                          aCGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTa                                          <  2:1895681/148‑1 (MQ=255)
                                                                                                           cgcgATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTAt                                         <  2:2276405/148‑1 (MQ=255)
                                                                                                              gATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGtt                                      >  2:596554/1‑148 (MQ=255)
                                                                                                               aTCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGACTCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTa                                     >  2:3261331/1‑148 (MQ=255)
                                                                                                               aTCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTGCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTa                                     >  2:1658156/1‑148 (MQ=255)
                                                                                                               aTCCATTAGATTATTGCTGCTTTGGCACTTGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTa                                     <  2:2277729/148‑1 (MQ=255)
                                                                                                                tCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCTATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTaa                                    <  2:463269/148‑1 (MQ=255)
                                                                                                                tCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTaa                                    >  1:819110/1‑148 (MQ=255)
                                                                                                                 ccATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAAt                                   <  2:2085385/148‑1 (MQ=255)
                                                                                                                    ttCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGc                                >  2:689388/1‑148 (MQ=255)
                                                                                                                     tCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCa                               >  2:940451/1‑148 (MQ=255)
                                                                                                                      cGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCaa                              <  1:1087157/148‑1 (MQ=255)
                                                                                                                         gTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGt                           <  1:2661457/148‑1 (MQ=255)
                                                                                                                                 gctATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAg                   >  2:1894612/1‑148 (MQ=255)
                                                                                                                                  ctATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGa                  <  1:1334/148‑1 (MQ=255)
                                                                                                                                   tATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGAt                 >  1:2515859/1‑148 (MQ=255)
                                                                                                                                     tGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGa               <  2:2502285/148‑1 (MQ=255)
                                                                                                                                        cACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAg            >  1:960493/1‑148 (MQ=255)
                                                                                                                                        cACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAg            <  1:603971/148‑1 (MQ=255)
                                                                                                                                          cTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCg          <  2:1019032/148‑1 (MQ=255)
                                                                                                                                            ggAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGtt        >  1:415955/1‑148 (MQ=255)
                                                                                                                                                  cAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCa  <  2:767176/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACAGAGTCTGGCGAATGGCACTCCCTATGTTAATCGCAATGTCAATGAAGATGACAGCGTTGAATCA  >  CP009273/1930894‑1931187

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: