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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273828,3740GT91.6% 656.5 / 47.5 215V11V (GTC→GTA) ybiADUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base G (11/7);  new base T (109/88);  total (120/95)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 216 total aligned reads displayed.

TTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTTC  >  CP009273/828229‑828523
                                                                                                                                                 |                                                                                                                                                     
       aaGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGAt                                                                                                                                            <  1:2838355/150‑1 (MQ=255)
       aaGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGAt                                                                                                                                            >  1:3865020/1‑150 (MQ=255)
        aGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATc                                                                                                                                           <  1:3087048/150‑1 (MQ=255)
        aGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATc                                                                                                                                           <  1:3751385/150‑1 (MQ=255)
         ggAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATcc                                                                                                                                          <  1:1669405/150‑1 (MQ=255)
         ggAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATcc                                                                                                                                          <  1:3857724/150‑1 (MQ=255)
           aaTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCtt                                                                                                                                        <  2:2358437/150‑1 (MQ=255)
           aaTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCtg                                                                                                                                        >  2:3474580/1‑149 (MQ=255)
            aTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGACAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCttt                                                                                                                                       >  1:1066131/1‑150 (MQ=255)
                  gtgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTc                                                                                                                                 <  1:2256157/150‑1 (MQ=255)
                  gtgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTc                                                                                                                                 >  1:3237590/1‑150 (MQ=255)
                  gtgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTc                                                                                                                                 >  1:3448231/1‑150 (MQ=255)
                    gCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGa                                                                                                                               <  1:3656298/150‑1 (MQ=255)
                    gCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGa                                                                                                                               <  2:2672409/150‑1 (MQ=255)
                      cTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAAc                                                                                                                             >  2:358449/1‑150 (MQ=255)
                       tGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACg                                                                                                                            >  2:3818545/1‑150 (MQ=255)
                         aaaaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACggg                                                                                                                          <  2:1044500/150‑1 (MQ=255)
                         aaaaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACggg                                                                                                                          >  1:3119580/1‑150 (MQ=255)
                          aaaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGc                                                                                                                         >  2:1387800/1‑150 (MQ=255)
                          aaaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGc                                                                                                                         >  2:1374116/1‑150 (MQ=255)
                           aaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCa                                                                                                                        >  1:766152/1‑150 (MQ=255)
                            aTAGTGTTCTGAGGTAGGCCTGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCATAGTTATAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCAt                                                                                                                       <  1:2263115/150‑1 (MQ=255)
                              aGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTa                                                                                                                     >  2:163219/1‑150 (MQ=255)
                                  ttCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAAcc                                                                                                                 <  2:438002/150‑1 (MQ=255)
                                    cTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGc                                                                                                               >  1:1986873/1‑150 (MQ=255)
                                     tGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCa                                                                                                              >  2:960498/1‑150 (MQ=255)
                                     tGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCa                                                                                                              >  2:1673419/1‑150 (MQ=255)
                                      gAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCAt                                                                                                             >  1:1875142/1‑150 (MQ=255)
                                      gAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCAt                                                                                                             >  1:3310203/1‑150 (MQ=255)
                                       aGGTAGGCCAGGTTGTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCAtt                                                                                                            <  2:2460747/150‑1 (MQ=255)
                                        ggTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTa                                                                                                           >  2:1508638/1‑150 (MQ=255)
                                           aGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGt                                                                                                        >  2:4272811/1‑150 (MQ=255)
                                               cAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggt                                                                                                    >  1:1124690/1‑150 (MQ=255)
                                                aGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtg                                                                                                   >  1:714571/1‑150 (MQ=255)
                                                aGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtg                                                                                                   >  1:1226592/1‑150 (MQ=255)
                                                  gTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTggtggt                                                                                                 <  1:3196729/150‑1 (MQ=255)
                                                    tttCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTa                                                                                               >  1:1445054/1‑150 (MQ=255)
                                                    tttCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTa                                                                                               <  1:3510195/150‑1 (MQ=255)
                                                      tCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAtt                                                                                             <  1:2816040/150‑1 (MQ=255)
                                                      tCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAtt                                                                                             <  2:1031543/150‑1 (MQ=255)
                                                         cGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCCTTATGTTGGTGGTTATTGCg                                                                                          >  2:3261812/1‑150 (MQ=255)
                                                          gTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGa                                                                                         <  1:1560410/150‑1 (MQ=255)
                                                           tCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAg                                                                                        <  1:2241568/150‑1 (MQ=255)
                                                             aaCCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGcc                                                                                      >  1:243985/1‑150 (MQ=255)
                                                             aaCCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGcc                                                                                      >  1:3580260/1‑150 (MQ=255)
                                                                 ttGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCtt                                                                                  >  1:3718552/1‑150 (MQ=255)
                                                                  tGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCttt                                                                                 >  1:2737801/1‑150 (MQ=255)
                                                                  tGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCttt                                                                                 <  1:3215955/150‑1 (MQ=255)
                                                                  tGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCttt                                                                                 <  2:2638941/150‑1 (MQ=255)
                                                                      gggCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcag                                                                             <  2:1777144/150‑1 (MQ=255)
                                                                      gggCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcag                                                                             <  1:180288/150‑1 (MQ=255)
                                                                         ccAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaa                                                                          >  2:3093000/1‑150 (MQ=255)
                                                                           aGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaaca                                                                        >  1:1383946/1‑150 (MQ=255)
                                                                           aGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaaca                                                                        <  1:950813/150‑1 (MQ=255)
                                                                           aGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaaca                                                                        >  2:1165914/1‑150 (MQ=255)
                                                                             gcggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacaga                                                                      <  1:1314983/150‑1 (MQ=255)
                                                                              cggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacagaa                                                                     >  2:408501/1‑150 (MQ=255)
                                                                               ggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacagaaa                                                                    >  2:3691458/1‑150 (MQ=255)
                                                                                gcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGaaaa                                                                   >  2:1249492/1‑150 (MQ=255)
                                                                                 cAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGaaaaa                                                                  >  2:4215453/1‑150 (MQ=255)
                                                                                  aaaGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAAc                                                                 <  2:3524231/150‑1 (MQ=255)
                                                                                    aGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTGTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCa                                                               <  2:1416865/150‑1 (MQ=255)
                                                                                     gTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAt                                                              <  1:3187535/150‑1 (MQ=255)
                                                                                      ttAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAtt                                                             <  1:1035036/150‑1 (MQ=255)
                                                                                      ttAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAtt                                                             <  1:4257073/150‑1 (MQ=255)
                                                                                      ttAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAtt                                                             <  1:3918475/150‑1 (MQ=255)
                                                                                        aGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTAc                                                           >  2:2297341/1‑150 (MQ=255)
                                                                                         gAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTAcc                                                          >  2:1973813/1‑150 (MQ=255)
                                                                                             aaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTg                                                      <  2:1149256/150‑1 (MQ=255)
                                                                                                  ccGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAAc                                                 <  2:2946997/150‑1 (MQ=255)
                                                                                                   cGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAAcc                                                >  1:321091/1‑150 (MQ=255)
                                                                                                       gTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGaaa                                            >  1:550393/1‑150 (MQ=255)
                                                                                                        tCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGATGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGaaaa                                           <  1:607174/150‑1 (MQ=255)
                                                                                                          aTCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAAt                                         >  1:532324/1‑150 (MQ=255)
                                                                                                          aTCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAAt                                         >  1:2396295/1‑150 (MQ=255)
                                                                                                            cGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGc                                       <  1:821618/150‑1 (MQ=255)
                                                                                                               tGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCcac                                    >  1:2857898/1‑150 (MQ=255)
                                                                                                                ggTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCcaca                                   <  2:2083482/150‑1 (MQ=255)
                                                                                                                ggTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCcaca                                   <  2:1470481/150‑1 (MQ=255)
                                                                                                                   gCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAata                                >  1:4094664/1‑150 (MQ=255)
                                                                                                                     tGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATAtt                              >  2:495964/1‑150 (MQ=255)
                                                                                                                     tGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATAtt                              <  1:2685147/150‑1 (MQ=255)
                                                                                                                      gTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTg                             >  1:631072/1‑150 (MQ=255)
                                                                                                                      gTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTg                             >  2:2291351/1‑150 (MQ=255)
                                                                                                                      gTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCAGTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTg                             >  2:765802/1‑150 (MQ=255)
                                                                                                                        aGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGc                           >  2:2102661/1‑150 (MQ=255)
                                                                                                                        aGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGc                           >  2:392657/1‑150 (MQ=255)
                                                                                                                         gAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCt                          >  2:1056221/1‑150 (MQ=255)
                                                                                                                            aTTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGtt                       >  1:3909533/1‑150 (MQ=255)
                                                                                                                                atgatgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtata                   >  1:2467574/1‑150 (MQ=255)
                                                                                                                                atgatgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtata                   >  1:349396/1‑150 (MQ=255)
                                                                                                                                 tgatgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATAc                  >  2:4168989/1‑150 (MQ=255)
                                                                                                                                  gatgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACa                 <  2:3712541/150‑1 (MQ=255)
                                                                                                                                  gatgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACa                 >  2:2797874/1‑150 (MQ=255)
                                                                                                                                   atgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAg                <  2:715005/150‑1 (MQ=255)
                                                                                                                                     ggTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTa              >  1:2358429/1‑150 (MQ=255)
                                                                                                                                        gtCTTGCATTACGGGTTGGATCCATTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAttt           <  2:3047563/150‑1 (MQ=255)
                                                                                                                                           ttGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAg        >  1:2865060/1‑150 (MQ=255)
                                                                                                                                           ttGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAg        >  2:2975287/1‑150 (MQ=255)
                                                                                                                                                 gACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTtc  <  2:3865020/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTTC  >  CP009273/828229‑828523

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