Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 828,374 G→T 92.4% V11V (GTC→GTA)  ybiA ← DUF1768 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273828,3740GT92.4% 739.1 / 48.2 236V11V (GTC→GTA) ybiADUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base G (7/11);  new base T (117/101);  total (124/112)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 236 total aligned reads displayed.

TTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTT  >  CP009273/828227‑828520
                                                                                                                                                   |                                                                                                                                                  
tttttCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTAc                                                                                                                                                  <  1:243570/150‑1 (MQ=255)
tttttCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTAc                                                                                                                                                  <  1:1712398/150‑1 (MQ=255)
 ttttCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACg                                                                                                                                                 >  1:2286128/1‑150 (MQ=255)
     catcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGtt                                                                                                                                             >  1:991189/1‑150 (MQ=255)
        caAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGa                                                                                                                                          <  2:88141/150‑1 (MQ=255)
          aGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATc                                                                                                                                        >  2:2885888/1‑150 (MQ=255)
             aaTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCtt                                                                                                                                     <  1:2684249/150‑1 (MQ=255)
               tttttGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCATTtg                                                                                                                                   >  2:619489/1‑150 (MQ=255)
                ttttGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTtgt                                                                                                                                  <  1:1424297/150‑1 (MQ=255)
                ttttGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTtgt                                                                                                                                  >  1:1503250/1‑150 (MQ=255)
                     tgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCg                                                                                                                             >  1:4160818/1‑150 (MQ=255)
                     tgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCg                                                                                                                             >  1:1229527/1‑150 (MQ=255)
                      gCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGa                                                                                                                            >  1:608568/1‑150 (MQ=255)
                       ccTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGaa                                                                                                                           >  1:1897141/1‑150 (MQ=255)
                        cTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAAc                                                                                                                          <  2:2380224/150‑1 (MQ=255)
                         tGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACg                                                                                                                         >  2:3197180/1‑150 (MQ=255)
                         tGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACg                                                                                                                         <  1:24970/150‑1 (MQ=255)
                             aaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCa                                                                                                                     >  2:2840061/1‑150 (MQ=255)
                              aTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCAt                                                                                                                    <  2:1840907/150‑1 (MQ=255)
                               tAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCAtt                                                                                                                   <  2:2413260/150‑1 (MQ=255)
                                aGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTa                                                                                                                  <  2:2911154/150‑1 (MQ=255)
                                 gtgtTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTaa                                                                                                                 >  2:1238921/1‑150 (MQ=255)
                                  tgtTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTaaa                                                                                                                <  2:4288279/150‑1 (MQ=255)
                                   gtTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAAc                                                                                                               >  1:4424959/1‑150 (MQ=255)
                                     tCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCg                                                                                                             >  1:886215/1‑150 (MQ=255)
                                     tCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCg                                                                                                             <  1:1035714/150‑1 (MQ=255)
                                     tCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTGGTGCTCGAACGGGCATTAAACCg                                                                                                             <  1:867408/150‑1 (MQ=255)
                                           gTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTAt                                                                                                       <  2:1700967/150‑1 (MQ=255)
                                            tAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATg                                                                                                      >  2:3557870/1‑150 (MQ=255)
                                             aGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGt                                                                                                     >  1:3680205/1‑150 (MQ=255)
                                               gCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtg                                                                                                   <  2:3853220/150‑1 (MQ=255)
                                                ccAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtgg                                                                                                  >  1:2114747/1‑150 (MQ=255)
                                                 cAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggt                                                                                                 <  1:4231543/150‑1 (MQ=255)
                                                 cAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggt                                                                                                 >  2:973028/1‑150 (MQ=255)
                                                  aGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtg                                                                                                <  2:2377934/150‑1 (MQ=255)
                                                   ggTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtgg                                                                                               >  1:1157659/1‑150 (MQ=255)
                                                   ggTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtgg                                                                                               >  2:2556663/1‑150 (MQ=255)
                                                   ggTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtgg                                                                                               >  2:3351562/1‑150 (MQ=255)
                                                     ttttCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGtt                                                                                             <  2:1854064/150‑1 (MQ=255)
                                                       ttCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAt                                                                                           >  2:112199/1‑150 (MQ=255)
                                                          ccGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGc                                                                                        >  2:2027551/1‑150 (MQ=255)
                                                                  cTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCt                                                                                <  2:1300021/150‑1 (MQ=255)
                                                                   ttGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCtt                                                                               >  1:819940/1‑150 (MQ=255)
                                                                        gggCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcag                                                                          >  1:4059530/1‑150 (MQ=255)
                                                                            cAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaac                                                                      >  1:3059605/1‑150 (MQ=255)
                                                                            cAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaac                                                                      <  1:2173593/150‑1 (MQ=255)
                                                                             aGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaaca                                                                     <  2:311319/150‑1 (MQ=255)
                                                                             aGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaaca                                                                     >  2:68009/1‑150 (MQ=255)
                                                                                cggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacagaa                                                                  <  1:3258231/150‑1 (MQ=255)
                                                                                cggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacagaa                                                                  <  1:3383524/150‑1 (MQ=255)
                                                                                cggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacagaa                                                                  >  1:3949644/1‑150 (MQ=255)
                                                                                      aGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCa                                                            >  1:3370196/1‑150 (MQ=255)
                                                                                      aGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCa                                                            >  1:2591539/1‑150 (MQ=255)
                                                                                      aGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCa                                                            >  1:3440290/1‑150 (MQ=255)
                                                                                        ttAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAtt                                                          <  2:3381501/150‑1 (MQ=255)
                                                                                              aaaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCt                                                    <  2:3082659/150‑1 (MQ=255)
                                                                                               aaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTg                                                   >  1:1929565/1‑150 (MQ=255)
                                                                                               aaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTg                                                   <  2:3920172/150‑1 (MQ=255)
                                                                                               aaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTg                                                   >  2:504964/1‑150 (MQ=255)
                                                                                               aaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTg                                                   <  2:1622958/150‑1 (MQ=255)
                                                                                                 tCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGACGCTTTCCAGAAACAGAAAAACCATTACCCCTGaa                                                 >  2:2159453/1‑150 (MQ=255)
                                                                                                   aCCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGaaaa                                               >  1:2849358/1‑150 (MQ=255)
                                                                                                   aCCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGaaaa                                               >  1:630821/1‑150 (MQ=255)
                                                                                                   aCCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGaaaa                                               >  2:322313/1‑150 (MQ=255)
                                                                                                    ccGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAAc                                              <  2:2690777/150‑1 (MQ=255)
                                                                                                           cATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGaaaaa                                       >  1:1954630/1‑150 (MQ=255)
                                                                                                            aTCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAAt                                      >  1:637688/1‑150 (MQ=255)
                                                                                                            aTCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGaaaaaa                                      >  1:3884271/1‑149 (MQ=255)
                                                                                                             tCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATg                                     <  1:4028278/150‑1 (MQ=255)
                                                                                                              cGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGc                                    >  1:3605218/1‑150 (MQ=255)
                                                                                                              cGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGc                                    <  2:108314/150‑1 (MQ=255)
                                                                                                               gCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGcc                                   >  1:289020/1‑150 (MQ=255)
                                                                                                               gCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGcc                                   >  2:1564092/1‑150 (MQ=255)
                                                                                                               gCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGcc                                   >  1:3573138/1‑150 (MQ=255)
                                                                                                                 tGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCcac                                 >  2:2447730/1‑150 (MQ=255)
                                                                                                                 tGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCcac                                 >  2:1412483/1‑150 (MQ=255)
                                                                                                                    tGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAat                              <  2:2617387/150‑1 (MQ=255)
                                                                                                                     gCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAata                             >  1:2099693/1‑150 (MQ=255)
                                                                                                                      cTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAatat                            >  2:383936/1‑150 (MQ=255)
                                                                                                                       tGTAGAACTTTAGGATGGTGTCTTGCATTACGTGTGGGCTCCTGTGAGCTCGAACGGGCACTAACCCGCATGATGATGGTGGTTAGTGCGAGCCGCTTTCCAGAAACAGAAAAACCCTTCCCCCCGAAAACCGCCAAAGGCCACCATAtt                           <  2:714246/150‑1 (MQ=255)
                                                                                                                       tGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATAtt                           <  2:4194778/150‑1 (MQ=255)
                                                                                                                        gTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTg                          <  1:1370760/150‑1 (MQ=255)
                                                                                                                         tAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTgg                         <  1:3552874/150‑1 (MQ=255)
                                                                                                                            aaaTTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTg                      <  1:1659518/150‑1 (MQ=255)
                                                                                                                            aaaTTTATGATGGTGTCTTGCATTACGTGTTGGATCCTGTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTg                      <  1:3365177/150‑1 (MQ=255)
                                                                                                                                ttATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTta                  >  1:2724584/1‑150 (MQ=255)
                                                                                                                                  atgatgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtata                >  2:4509159/1‑150 (MQ=255)
                                                                                                                                  atgatgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtata                >  2:1971717/1‑150 (MQ=255)
                                                                                                                                      tgGTGTCTTGCATTACGTGTTGGATCCTGTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGt            <  1:3135888/150‑1 (MQ=255)
                                                                                                                                       ggTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTa           <  2:4424959/150‑1 (MQ=255)
                                                                                                                                       ggTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTa           >  2:2061503/1‑150 (MQ=255)
                                                                                                                                       ggTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTa           >  1:2395598/1‑150 (MQ=255)
                                                                                                                                         tgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAtt         <  1:4381630/150‑1 (MQ=255)
                                                                                                                                         tgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAtt         <  1:505643/150‑1 (MQ=255)
                                                                                                                                         tgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGGTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAtt         <  2:3458427/150‑1 (MQ=255)
                                                                                                                                          gtCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAttt        <  1:1600550/150‑1 (MQ=255)
                                                                                                                                            cTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCa      >  2:1504269/1‑150 (MQ=255)
                                                                                                                                            cTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCa      >  1:2232087/1‑150 (MQ=255)
                                                                                                                                              tGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAgg    <  2:4051447/150‑1 (MQ=255)
                                                                                                                                              tGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAgg    <  1:3336006/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTT  >  CP009273/828227‑828520

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: ‑