Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,175,657 G→T 11.8% E191* (GAA→TAA)  cobB → deacetylase of acs and cheY, chemotaxis regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,175,6570GT11.8% 724.0 / 49.1 227E191* (GAA→TAA) cobBdeacetylase of acs and cheY, chemotaxis regulator
Reads supporting (aligned to +/- strand):  ref base G (72/128);  new base T (16/10);  total (89/138)
Fisher's exact test for biased strand distribution p-value = 1.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

Only 100 of 229 total aligned reads displayed.

TTTCCGGCACCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGA  >  CP009273/1175510‑1175803
                                                                                                                                                   |                                                                                                                                                  
        aCCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGaaa                                                                                                                                            <  2:2939059/148‑1 (MQ=255)
        aCCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGaaa                                                                                                                                            <  2:3202597/148‑1 (MQ=255)
        aCCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGaaa                                                                                                                                            <  1:3182344/148‑1 (MQ=255)
        aCCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGaaa                                                                                                                                            <  2:2496071/148‑1 (MQ=255)
         cccTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAAc                                                                                                                                           >  1:988947/1‑148 (MQ=255)
          ccTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACt                                                                                                                                          >  2:2342191/1‑148 (MQ=255)
             tGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCa                                                                                                                                       >  2:866684/1‑148 (MQ=255)
                  ccGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGgcg                                                                                                                                  >  1:2839016/1‑148 (MQ=255)
                        gTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGcaca                                                                                                                            <  1:667129/148‑1 (MQ=255)
                        gTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGcaca                                                                                                                            <  2:3444436/148‑1 (MQ=255)
                                gTTTGGCGAAAGGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGaa                                                                                                                    <  2:2306124/148‑1 (MQ=255)
                                  ttGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACt                                                                                                                  >  1:1709154/1‑148 (MQ=255)
                                  ttGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACt                                                                                                                  <  2:2675030/148‑1 (MQ=255)
                                    ggCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGa                                                                                                                <  1:327425/148‑1 (MQ=255)
                                      cGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAAt                                                                                                              <  2:354413/148‑1 (MQ=255)
                                       gAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATc                                                                                                             <  2:3349414/148‑1 (MQ=255)
                                       gAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATc                                                                                                             <  2:3785553/148‑1 (MQ=255)
                                       gAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATc                                                                                                             <  2:3811933/148‑1 (MQ=255)
                                         aaTGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCtt                                                                                                           >  2:46140/1‑148 (MQ=255)
                                          aTGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTg                                                                                                          >  1:1858907/1‑148 (MQ=255)
                                          aTGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTg                                                                                                          <  2:1922680/148‑1 (MQ=255)
                                          aTGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTg                                                                                                          >  2:1370487/1‑148 (MQ=255)
                                           tGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGa                                                                                                         <  1:310309/148‑1 (MQ=255)
                                            gCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGaa                                                                                                        >  1:2221070/1‑148 (MQ=255)
                                             acACTCGTCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGACTTTGGCACATGCGGGCATGTTTATCCGGGGGCTGGGTTTGTTCACGAAGAGAAACTGCATGGCGCGCACACAGTTGATCTGAATCTTGAAc                                                                                                       <  2:2986618/147‑1 (MQ=255)
                                              cACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAAcc                                                                                                      >  1:114086/1‑148 (MQ=255)
                                                    gCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCa                                                                                                <  1:2533982/148‑1 (MQ=255)
                                                        ggATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGtt                                                                                            >  1:1520243/1‑148 (MQ=255)
                                                         gATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTg                                                                                           >  2:2934016/1‑148 (MQ=255)
                                                         gATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTg                                                                                           >  1:415458/1‑148 (MQ=255)
                                                         gATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTg                                                                                           >  2:376465/1‑148 (MQ=255)
                                                         gATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCAGTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTg                                                                                           <  2:2717312/148‑1 (MQ=255)
                                                          aTGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTgg                                                                                          <  2:3472402/148‑1 (MQ=255)
                                                           tGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGt                                                                                         >  1:1276497/1‑148 (MQ=255)
                                                           tGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGt                                                                                         >  2:3625499/1‑148 (MQ=255)
                                                            gAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTa                                                                                        <  1:438856/148‑1 (MQ=255)
                                                              aaTTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAAt                                                                                      <  1:1236503/148‑1 (MQ=255)
                                                              aaTTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAAt                                                                                      <  1:1236242/148‑1 (MQ=255)
                                                               aTTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATg                                                                                     >  2:3032977/1‑148 (MQ=255)
                                                                tttATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGa                                                                                    >  1:1830444/1‑148 (MQ=255)
                                                                  tataTGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACCTCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAAt                                                                                  <  2:1797201/148‑1 (MQ=255)
                                                                  tataTGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAAt                                                                                  <  1:2979484/148‑1 (MQ=255)
                                                                  tataTGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAAt                                                                                  <  1:1361930/148‑1 (MQ=255)
                                                                  tataTGGCGTTGTCGATGGCCGATATTTTCATTGCCAGTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAAt                                                                                  <  2:1345417/148‑1 (MQ=255)
                                                                           ttGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCgaga                                                                         >  1:367491/1‑148 (MQ=255)
                                                                            tGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGaa                                                                        <  2:1829457/148‑1 (MQ=255)
                                                                              tCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAat                                                                      <  1:2331585/148‑1 (MQ=255)
                                                                               cGATGGCCGATATTTTCATTGCCAGTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGTGTCAGGTTGGTAATGAATTTGCCGAGAAata                                                                     <  2:1177397/148‑1 (MQ=255)
                                                                                gATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGGGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAatat                                                                    <  1:699535/148‑1 (MQ=255)
                                                                                       gATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGcc                                                             <  1:3342911/148‑1 (MQ=255)
                                                                                         tatTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCg                                                           >  1:2006721/1‑148 (MQ=255)
                                                                                         tatTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCg                                                           >  2:2147528/1‑148 (MQ=255)
                                                                                           ttttCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGc                                                         >  1:3142840/1‑148 (MQ=255)
                                                                                           ttttCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGc                                                         >  2:3194774/1‑148 (MQ=255)
                                                                                            tttCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCGCGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCa                                                        <  1:838818/148‑1 (MQ=255)
                                                                                            tttCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGACATATTACGGCCCGGCa                                                        >  1:793364/1‑148 (MQ=255)
                                                                                              tCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAg                                                      <  1:2332708/148‑1 (MQ=255)
                                                                                               cATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGc                                                     <  2:3115960/148‑1 (MQ=255)
                                                                                                aTTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCCTGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGcc                                                    >  1:3564391/1‑148 (MQ=255)
                                                                                                 ttGCCATTGGCACATCCGGGCATGTTTCTCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCa                                                   <  2:335517/148‑1 (MQ=255)
                                                                                                 ttGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCa                                                   >  1:583354/1‑148 (MQ=255)
                                                                                                 ttGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCa                                                   <  2:694529/148‑1 (MQ=255)
                                                                                                 ttGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCa                                                   >  1:2799145/1‑148 (MQ=255)
                                                                                                  tGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAg                                                  >  2:1283946/1‑148 (MQ=255)
                                                                                                   gCCATTGGCACATCCGGGCATGTTTATGCGGCGGCTGGGTTTTTTCACGAAGCGAAACTGCATGGCGCGCACTCCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTTCCGAGAAATATTACGGCCCGGCAAGCCAgg                                                 <  2:1331970/148‑1 (MQ=255)
                                                                                                   gCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAgg                                                 <  2:1704303/148‑1 (MQ=255)
                                                                                                   gCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAgg                                                 <  2:1331892/148‑1 (MQ=255)
                                                                                                   gCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAgg                                                 <  2:3303589/148‑1 (MQ=255)
                                                                                                     cATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAggtg                                               <  1:926833/148‑1 (MQ=255)
                                                                                                      aTTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAggtgg                                              >  2:2446975/1‑148 (MQ=255)
                                                                                                          gCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGcc                                          <  1:632468/148‑1 (MQ=255)
                                                                                                           cacaTCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGACCTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCt                                         >  2:613679/1‑148 (MQ=255)
                                                                                                             caTCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGTGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGa                                       >  1:658435/1‑148 (MQ=255)
                                                                                                               tCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGt                                     <  1:1627254/148‑1 (MQ=255)
                                                                                                                ccGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGtt                                    <  2:3484105/148‑1 (MQ=255)
                                                                                                                ccGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGtt                                    <  1:3476718/148‑1 (MQ=255)
                                                                                                                ccGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGtt                                    <  2:1428643/148‑1 (MQ=255)
                                                                                                                ccGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGtt                                    <  2:1586610/148‑1 (MQ=255)
                                                                                                                   ggCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTttgt                                 <  1:3424434/148‑1 (MQ=255)
                                                                                                                   ggCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTttgt                                 <  2:1811998/148‑1 (MQ=255)
                                                                                                                   ggCATGTTGATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTttgt                                 <  2:93469/148‑1 (MQ=255)
                                                                                                                    gCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTttgtt                                <  1:2412698/148‑1 (MQ=255)
                                                                                                                    gCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTttgtt                                >  1:2143131/1‑148 (MQ=255)
                                                                                                                     cATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTttgttg                               >  1:2599372/1‑148 (MQ=255)
                                                                                                                        gTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGaaa                            <  1:3238467/148‑1 (MQ=255)
                                                                                                                        gTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGaaa                            >  2:1245591/1‑148 (MQ=255)
                                                                                                                          ttATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAg                          <  2:1762706/148‑1 (MQ=255)
                                                                                                                           tATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGt                         <  2:1261150/148‑1 (MQ=255)
                                                                                                                           tATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGt                         <  2:2060102/148‑1 (MQ=255)
                                                                                                                              ccGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGc                      <  2:114101/148‑1 (MQ=255)
                                                                                                                              ccGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGc                      >  1:2392956/1‑148 (MQ=255)
                                                                                                                               cggcggCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCt                     <  1:1116328/148‑1 (MQ=255)
                                                                                                                                   ggcTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAg                 <  2:2670721/148‑1 (MQ=255)
                                                                                                                                       gggTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGAt             >  2:3721251/1‑148 (MQ=255)
                                                                                                                                       gggTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGAt             >  1:2270034/1‑148 (MQ=255)
                                                                                                                                       gggTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGAt             <  1:903163/148‑1 (MQ=255)
                                                                                                                                           ttGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTaaa         >  1:2009358/1‑148 (MQ=255)
                                                                                                                                             gTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAg       >  2:1016897/1‑148 (MQ=255)
                                                                                                                                             gTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAg       <  1:3803684/148‑1 (MQ=255)
                                                                                                                                                  cGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGa  <  2:1986676/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTTCCGGCACCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGA  >  CP009273/1175510‑1175803

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: