Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 868,500 A→T 90.5% S22S (TCA→TCT yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273868,5000AT90.5% 827.3 / 63.8 275S22S (TCA→TCTyliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (11/15);  new base T (138/111);  total (149/126)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

Only 100 of 277 total aligned reads displayed.

ATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGACTCAC  >  CP009273/868355‑868648
                                                                                                                                                 |                                                                                                                                                    
aTTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTTCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGTCTGATAATTTTATTTTTATTGGCAGCAATGTCtttt                                                                                                                                                   <  1:1511490/149‑1 (MQ=255)
  ttGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTAt                                                                                                                                                 <  2:4001355/149‑1 (MQ=255)
      ttAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGt                                                                                                                                             >  1:2317728/1‑149 (MQ=255)
       tAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTc                                                                                                                                            >  1:4575015/1‑149 (MQ=255)
        aattaGCCGCAGCAATGGACAATTTTTGCAGCGTCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTTTCTTTTATTGGTCt                                                                                                                                           <  2:883199/145‑1 (MQ=255)
        aaaaaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCt                                                                                                                                           <  2:3902071/149‑1 (MQ=255)
         aaaaGCCGCAGCAAGGGACAATTTTTGCAGCGTCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCtt                                                                                                                                          <  1:4679685/149‑1 (MQ=255)
          aaaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttt                                                                                                                                         >  1:3474561/1‑149 (MQ=255)
          aaaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttt                                                                                                                                         >  1:598454/1‑149 (MQ=255)
           aaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCtttt                                                                                                                                        >  2:3247036/1‑149 (MQ=255)
            aGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttt                                                                                                                                       <  1:2160535/149‑1 (MQ=255)
             gCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCtttttt                                                                                                                                      >  2:358533/1‑149 (MQ=255)
             gCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTAGTGGTCtttttt                                                                                                                                      >  1:2906095/1‑149 (MQ=255)
                 cagcaAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTtcat                                                                                                                                  <  1:2769024/149‑1 (MQ=255)
                    caAGGGACAATTTTAGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGTCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTcatcat                                                                                                                               <  2:2317373/149‑1 (MQ=255)
                     aaGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCAtt                                                                                                                              <  2:2162456/149‑1 (MQ=255)
                     aaGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCAtt                                                                                                                              >  2:2497457/1‑149 (MQ=255)
                       gggACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTaa                                                                                                                            >  2:1506162/1‑149 (MQ=255)
                       gggACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTaa                                                                                                                            >  2:535457/1‑149 (MQ=255)
                        ggACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAAc                                                                                                                           <  1:2640673/149‑1 (MQ=255)
                         gACAATTTTTGCAGCGTCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACt                                                                                                                          <  1:4128304/149‑1 (MQ=255)
                         gACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACt                                                                                                                          <  1:1739785/149‑1 (MQ=255)
                                attGCAGCGGCACAGCGTTCAGATAGTTATTTTTTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACt                                                                                                                   <  2:1406091/148‑1 (MQ=255)
                                  tGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTgg                                                                                                                 <  1:4165054/149‑1 (MQ=255)
                                   gCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGt                                                                                                                <  1:1227762/149‑1 (MQ=255)
                                    cAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTa                                                                                                               <  1:535457/149‑1 (MQ=255)
                                    cAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTa                                                                                                               <  1:1110488/149‑1 (MQ=255)
                                         gCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGa                                                                                                          >  1:2358760/1‑149 (MQ=255)
                                                gTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGcgcgcg                                                                                                   >  1:2947145/1‑149 (MQ=255)
                                                gTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGcgcgcg                                                                                                   >  1:3170959/1‑149 (MQ=255)
                                                  tCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGt                                                                                                 <  1:2269486/149‑1 (MQ=255)
                                                    aGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaa                                                                                               >  2:477850/1‑149 (MQ=255)
                                                     gATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTAAACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaa                                                                                              >  1:3697888/1‑149 (MQ=255)
                                                         gTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGt                                                                                          >  2:521623/1‑149 (MQ=255)
                                                         gTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGt                                                                                          >  2:4371037/1‑149 (MQ=255)
                                                           tATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGc                                                                                        >  1:4487014/1‑149 (MQ=255)
                                                           tATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGc                                                                                        <  1:3609973/149‑1 (MQ=255)
                                                                  ttAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGc                                                                                 <  2:1763567/149‑1 (MQ=255)
                                                                  ttAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGc                                                                                 >  1:668449/1‑149 (MQ=255)
                                                                   tAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGcc                                                                                <  2:3071644/149‑1 (MQ=255)
                                                                    aaaTGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCg                                                                               >  1:1146778/1‑149 (MQ=255)
                                                                     aaTGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGc                                                                              >  2:4760086/1‑149 (MQ=255)
                                                                       tGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTATTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCtt                                                                            <  1:3811656/149‑1 (MQ=255)
                                                                         tATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTg                                                                          >  2:617740/1‑149 (MQ=255)
                                                                          aTTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGa                                                                         >  2:1220794/1‑149 (MQ=255)
                                                                          aTTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGa                                                                         >  2:1427872/1‑149 (MQ=255)
                                                                             aaCATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACt                                                                      <  1:1808622/149‑1 (MQ=255)
                                                                             aaCATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACt                                                                      >  2:4676237/1‑149 (MQ=255)
                                                                               cATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTa                                                                    >  1:2707668/1‑149 (MQ=255)
                                                                                 tGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTAtt                                                                  <  2:2228037/149‑1 (MQ=255)
                                                                                 tGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTAtt                                                                  <  1:4745106/149‑1 (MQ=255)
                                                                                  gCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTc                                                                 >  1:3989930/1‑149 (MQ=255)
                                                                                   cTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCa                                                                >  1:2436177/1‑149 (MQ=255)
                                                                                     gAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGa                                                              <  2:2780632/149‑1 (MQ=255)
                                                                                         ttATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaaaaa                                                          <  1:3050445/149‑1 (MQ=255)
                                                                                          tataCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAAc                                                         <  2:1680659/149‑1 (MQ=255)
                                                                                           ataCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACg                                                        >  2:414101/1‑149 (MQ=255)
                                                                                            taCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGt                                                       >  1:855946/1‑149 (MQ=255)
                                                                                            taCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGt                                                       <  1:2812451/149‑1 (MQ=255)
                                                                                              cGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTg                                                     >  2:4193864/1‑149 (MQ=255)
                                                                                               gAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTgg                                                    >  1:3198818/1‑149 (MQ=255)
                                                                                                aaaaGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGc                                                   >  1:1455944/1‑149 (MQ=255)
                                                                                                 aaagataaagatAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCt                                                  <  2:2400614/149‑1 (MQ=255)
                                                                                                  aagataaagataAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTa                                                 <  1:4264667/149‑1 (MQ=255)
                                                                                                   agataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATAGGTCTTTTTTTCATCATTAACAATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACATATTCAGAAAAACGTTGGCTAt                                                >  1:4410526/1‑149 (MQ=255)
                                                                                                     ataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAttt                                              >  2:1033111/1‑149 (MQ=255)
                                                                                                     ataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAttt                                              <  2:1995084/149‑1 (MQ=255)
                                                                                                        aagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCtt                                           >  1:2546144/1‑149 (MQ=255)
                                                                                                         agataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTAAGAAAAACGTTGGCTATTTCttt                                          >  1:4076568/1‑149 (MQ=255)
                                                                                                           ataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaa                                        >  2:1883400/1‑149 (MQ=255)
                                                                                                               ggCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGat                                    <  2:1285444/149‑1 (MQ=255)
                                                                                                               ggCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTCTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGat                                    >  1:2212769/1‑149 (MQ=255)
                                                                                                               ggCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGat                                    <  2:4556822/149‑1 (MQ=255)
                                                                                                                gCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGata                                   <  1:2550451/149‑1 (MQ=255)
                                                                                                                 cTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGatat                                  >  1:2651561/1‑149 (MQ=255)
                                                                                                                 cTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGatat                                  >  2:2967893/1‑149 (MQ=255)
                                                                                                                 cTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGatat                                  >  1:1102283/1‑149 (MQ=255)
                                                                                                                  tGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATAAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTTGCTATTTCTTTAAAGATAtt                                 >  1:1390199/1‑149 (MQ=255)
                                                                                                                   gATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTg                                >  1:705477/1‑149 (MQ=255)
                                                                                                                    aTAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCAATAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGAGCCGATATCCGCTTAGAACTAATTCAGAAAAAAGTTGGCTATTTATTTAAAGATATTGa                               >  1:2601175/1‑149 (MQ=255)
                                                                                                                      aaTTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAAc                             >  2:3013809/1‑149 (MQ=255)
                                                                                                                       attttattttTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACg                            <  1:1777194/149‑1 (MQ=255)
                                                                                                                        ttttatttttatTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGt                           <  1:957574/149‑1 (MQ=255)
                                                                                                                        ttttatttttatTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGt                           <  1:3308918/149‑1 (MQ=255)
                                                                                                                         tttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                          <  2:979277/149‑1 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTc                         <  1:2444511/149‑1 (MQ=255)
                                                                                                                             tttttattGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGc                      >  2:4460072/1‑149 (MQ=255)
                                                                                                                                 tattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTg                  >  1:3044009/1‑149 (MQ=255)
                                                                                                                                  attGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGa                 <  2:3808147/149‑1 (MQ=255)
                                                                                                                                     ggCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACAt              <  2:2011507/149‑1 (MQ=255)
                                                                                                                                        agcaCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTaa           >  2:3467900/1‑149 (MQ=255)
                                                                                                                                          caCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAg         <  1:1671223/149‑1 (MQ=255)
                                                                                                                                          caCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAg         <  2:1704152/149‑1 (MQ=255)
                                                                                                                                            cTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa       >  1:2106972/1‑149 (MQ=255)
                                                                                                                                            cTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa       <  2:2660531/149‑1 (MQ=255)
                                                                                                                                            cTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa       >  2:1018288/1‑149 (MQ=255)
                                                                                                                                            cTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa       >  1:4558452/1‑149 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGAc      <  1:4119721/149‑1 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGAc      >  1:4358391/1‑149 (MQ=255)
                                                                                                                                              gTCTTTTATTGGTCTTTTTTTCATCATTTACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGATTTGAACTTATTCAGAAAAACGTTGGATATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGACt     <  1:3354115/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGACTCAC  >  CP009273/868355‑868648

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: