Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,885,611 G→T 80.0% R628S (CGT→AGT)  yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,885,6110GT80.0% 583.1 / 194.7 309R628S (CGT→AGT) yoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base T (140/107);  total (174/135)
Fisher's exact test for biased strand distribution p-value = 8.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 310 total aligned reads displayed.

GTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGT  >  CP009273/1885465‑1885759
                                                                                                                                                  |                                                                                                                                                    
gtcgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGaa                                                                                                                                                    >  1:2102945/1‑149 (MQ=255)
 tcgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAAc                                                                                                                                                   <  2:1198643/149‑1 (MQ=255)
  cgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAAcc                                                                                                                                                  <  1:3261460/149‑1 (MQ=255)
  cgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAAcc                                                                                                                                                  <  1:1167112/149‑1 (MQ=255)
    ttCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGc                                                                                                                                                >  1:4444845/1‑149 (MQ=255)
      ccACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCAc                                                                                                                                              >  1:1784553/1‑149 (MQ=255)
        acaGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACgg                                                                                                                                            <  2:3431352/149‑1 (MQ=255)
        acaGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACgg                                                                                                                                            <  2:3370800/149‑1 (MQ=255)
          aGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGc                                                                                                                                          >  2:2153648/1‑149 (MQ=255)
            ggCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCaa                                                                                                                                        <  2:4263505/149‑1 (MQ=255)
                 cAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCa                                                                                                                                   <  2:434086/149‑1 (MQ=255)
                  agagCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCAc                                                                                                                                  >  1:919169/1‑149 (MQ=255)
                   gagCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACg                                                                                                                                 >  2:1547502/1‑149 (MQ=255)
                     gCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACggg                                                                                                                               <  1:313733/149‑1 (MQ=255)
                     gCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACggg                                                                                                                               <  1:198557/149‑1 (MQ=255)
                      cAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGt                                                                                                                              <  2:4534140/149‑1 (MQ=255)
                        ggCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGc                                                                                                                            >  1:199966/1‑149 (MQ=255)
                        ggCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATGACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGc                                                                                                                            <  1:946604/149‑1 (MQ=255)
                          ccTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGt                                                                                                                          <  2:2068150/149‑1 (MQ=255)
                               gTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAAGTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAAGGTCACGGGTGCGTGgcgc                                                                                                                     <  1:299963/149‑1 (MQ=255)
                                 cgcgGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcgg                                                                                                                   <  1:748594/149‑1 (MQ=255)
                                 cgcgGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCgcgg                                                                                                                   <  1:3493350/149‑1 (MQ=255)
                                  gcgGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcggg                                                                                                                  >  2:2913288/1‑149 (MQ=255)
                                  gcgGTATCGATAGCCAGAAGTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCAGCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcggg                                                                                                                  <  1:1381836/149‑1 (MQ=255)
                                    ggTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcgggcg                                                                                                                >  2:3190479/1‑149 (MQ=255)
                                      tATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGggcggc                                                                                                              >  1:2587416/1‑149 (MQ=255)
                                          gATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCACTGTCACGGGTGCGTGGCGCGGGCGGCAGAc                                                                                                          >  1:1705189/1‑149 (MQ=255)
                                             aGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTgg                                                                                                       >  1:4553755/1‑149 (MQ=255)
                                                  gAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCgaga                                                                                                  >  2:3680035/1‑149 (MQ=255)
                                                  gAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCgaga                                                                                                  >  2:1081010/1‑149 (MQ=255)
                                                  gAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCgaga                                                                                                  >  2:1081053/1‑149 (MQ=255)
                                                   aaTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGaa                                                                                                 <  1:2031386/149‑1 (MQ=255)
                                                    aTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGaaa                                                                                                <  1:2121249/149‑1 (MQ=255)
                                                    aTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGaaa                                                                                                >  1:4378538/1‑149 (MQ=255)
                                                    aTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGaaa                                                                                                <  1:4531193/149‑1 (MQ=255)
                                                     ttCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAAc                                                                                               >  1:372634/1‑149 (MQ=255)
                                                        gCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTc                                                                                            >  2:4005433/1‑149 (MQ=255)
                                                        gCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTc                                                                                            <  2:4639092/149‑1 (MQ=255)
                                                        gCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTc                                                                                            <  2:2242726/149‑1 (MQ=255)
                                                           tGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAAATACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTcgcg                                                                                         >  1:4474388/1‑149 (MQ=255)
                                                            gggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCgcgc                                                                                        >  1:2226578/1‑149 (MQ=255)
                                                            gggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCgcgc                                                                                        >  2:2277470/1‑149 (MQ=255)
                                                             ggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGcc                                                                                       >  2:342477/1‑149 (MQ=255)
                                                                 ttCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCACTGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTaa                                                                                   <  2:4390906/149‑1 (MQ=255)
                                                                 ttCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTaa                                                                                   <  2:546927/149‑1 (MQ=255)
                                                                  tCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAg                                                                                  <  1:100315/149‑1 (MQ=255)
                                                                     tGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGAc                                                                               >  1:3856388/1‑149 (MQ=255)
                                                                     tGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGAc                                                                               >  1:1661808/1‑149 (MQ=255)
                                                                      gCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAATGTCGCGCCGTAAGGACg                                                                              >  2:1738710/1‑149 (MQ=255)
                                                                      gCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACg                                                                              >  2:4179641/1‑149 (MQ=255)
                                                                        ttAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCa                                                                            <  2:201139/149‑1 (MQ=255)
                                                                        ttAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCa                                                                            <  1:3999870/149‑1 (MQ=255)
                                                                          aGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATc                                                                          >  2:1427455/1‑149 (MQ=255)
                                                                           gATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCa                                                                         >  2:3398922/1‑149 (MQ=255)
                                                                            aTCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCAc                                                                        <  2:1156294/149‑1 (MQ=255)
                                                                             tCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCAcc                                                                       <  2:2447850/149‑1 (MQ=255)
                                                                               aaTAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAg                                                                     >  1:79005/1‑149 (MQ=255)
                                                                                aTAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGc                                                                    >  2:1265472/1‑149 (MQ=255)
                                                                                 tAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGcc                                                                   >  1:480644/1‑149 (MQ=255)
                                                                                 gaaaaaGGCGCGCATCATACCAGACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGcc                                                                   <  2:776183/148‑1 (MQ=255)
                                                                                  aaaaaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCg                                                                  <  2:3844595/149‑1 (MQ=255)
                                                                                   aaaaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGa                                                                 >  2:2734483/1‑149 (MQ=255)
                                                                                        gcgcgcATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTc                                                            >  2:543790/1‑149 (MQ=255)
                                                                                             catcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACaaa                                                       >  1:875433/1‑149 (MQ=255)
                                                                                             catcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACaaa                                                       <  1:2405164/149‑1 (MQ=255)
                                                                                              atcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAAt                                                      >  2:936951/1‑149 (MQ=255)
                                                                                               tcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATc                                                     >  1:2740277/1‑149 (MQ=255)
                                                                                               tcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATc                                                     <  1:1037298/149‑1 (MQ=255)
                                                                                                  taccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAcc                                                  >  2:180481/1‑149 (MQ=255)
                                                                                                   accatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCa                                                 >  2:3985856/1‑149 (MQ=255)
                                                                                                      atacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCa                                              >  2:1305207/1‑149 (MQ=255)
                                                                                                       tacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCAc                                             >  1:1879412/1‑149 (MQ=255)
                                                                                                       tacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCAc                                             >  1:4269690/1‑149 (MQ=255)
                                                                                                         cTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGc                                           >  2:4151929/1‑149 (MQ=255)
                                                                                                          tCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGcc                                          >  1:1092734/1‑149 (MQ=255)
                                                                                                              tAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGa                                      <  1:739382/149‑1 (MQ=255)
                                                                                                              tAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGa                                      >  2:2379817/1‑149 (MQ=255)
                                                                                                              tAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGa                                      >  1:680787/1‑149 (MQ=255)
                                                                                                                aCAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtc                                    >  2:3248359/1‑149 (MQ=255)
                                                                                                                aCAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtc                                    >  2:4102425/1‑149 (MQ=255)
                                                                                                                 cAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcg                                   >  2:2140043/1‑149 (MQ=255)
                                                                                                                  aaaTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgt                                  <  1:2652786/149‑1 (MQ=255)
                                                                                                                  aaaTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgt                                  <  2:3040031/149‑1 (MQ=255)
                                                                                                                  aaaTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgt                                  <  2:876347/149‑1 (MQ=255)
                                                                                                                   aaTTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgtc                                 <  1:2277470/149‑1 (MQ=255)
                                                                                                                     ttACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCgg                               >  2:653943/1‑149 (MQ=255)
                                                                                                                      tACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGc                              <  2:371862/149‑1 (MQ=255)
                                                                                                                       aCCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCg                             >  1:1322032/1‑149 (MQ=255)
                                                                                                                        ccTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGt                            <  1:2524676/149‑1 (MQ=255)
                                                                                                                           ggaggaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTcgc                         >  2:3926076/1‑149 (MQ=255)
                                                                                                                              ggaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGac                      >  2:1767731/1‑148 (MQ=255)
                                                                                                                              ggaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGaa                      >  1:1610155/1‑149 (MQ=255)
                                                                                                                                aTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATc                    <  2:3562746/149‑1 (MQ=255)
                                                                                                                                 tGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCa                   <  1:4320273/149‑1 (MQ=255)
                                                                                                                                  ggTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAg                  >  1:4271881/1‑149 (MQ=255)
                                                                                                                                      tCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGc              >  1:2380617/1‑149 (MQ=255)
                                                                                                                                         cAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCt           >  1:283476/1‑149 (MQ=255)
                                                                                                                                           aGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTAc         <  1:557466/149‑1 (MQ=255)
                                                                                                                                           aGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTAc         >  1:3406212/1‑149 (MQ=255)
                                                                                                                                               aaCTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTAccccc     >  2:1017738/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGT  >  CP009273/1885465‑1885759

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: