Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 961,946 T→A 8.6% F724Y (TTT→TAT)  ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273961,9460TA 8.6% 1098.1 / 73.6 312F724Y (TTT→TAT) ycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (141/143);  new base A (12/15);  total (153/159)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 323 total aligned reads displayed.

ATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGA  >  CP009273/961798‑962093
                                                                                                                                                    |                                                                                                                                                   
aTTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAAttt                                                                                                                                                     <  2:3432426/149‑1 (MQ=255)
 ttCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAAtttt                                                                                                                                                    <  1:782289/149‑1 (MQ=255)
   cAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTc                                                                                                                                                  <  1:589236/149‑1 (MQ=255)
       gAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAgcg                                                                                                                                              >  2:2615813/1‑149 (MQ=255)
             aTGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGc                                                                                                                                        >  1:1167252/1‑149 (MQ=255)
               ggAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTg                                                                                                                                      >  1:2344100/1‑149 (MQ=255)
                gAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTgg                                                                                                                                     <  1:1830203/149‑1 (MQ=255)
                   aaGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCgg                                                                                                                                  <  2:4171814/149‑1 (MQ=255)
                   aaGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCgg                                                                                                                                  >  2:152926/1‑149 (MQ=255)
                     gTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGAt                                                                                                                                <  1:4431073/149‑1 (MQ=255)
                      tGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGAtt                                                                                                                               <  2:2425289/149‑1 (MQ=255)
                      tGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGAtt                                                                                                                               <  1:1213525/149‑1 (MQ=255)
                        tcACTCGCATCGGCATCGCTCTATAACGCATGGCGACTGCCCTAGACCATAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAATGTCAAACAACGGTCAATTTTTCAGCGCTAGGCTGGCGGATTag                                                                                                                             <  2:1750772/148‑1 (MQ=255)
                        gCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTag                                                                                                                             >  1:2725936/1‑149 (MQ=255)
                                  cggcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAgg                                                                                                                   >  2:1426431/1‑149 (MQ=255)
                                   ggcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAggg                                                                                                                  <  1:3696755/149‑1 (MQ=255)
                                    gcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGa                                                                                                                 >  1:4479964/1‑149 (MQ=255)
                                    gcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGa                                                                                                                 >  1:3649954/1‑149 (MQ=255)
                                      atcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGc                                                                                                               >  1:4068793/1‑149 (MQ=255)
                                       tcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCa                                                                                                              >  1:3380743/1‑149 (MQ=255)
                                         gcgcTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCaaa                                                                                                            >  2:2979871/1‑149 (MQ=255)
                                            cTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAAttt                                                                                                         >  2:2516389/1‑149 (MQ=255)
                                             tataACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAAtttt                                                                                                        >  2:1302161/1‑149 (MQ=255)
                                             tataACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAAtttt                                                                                                        >  2:2149172/1‑149 (MQ=255)
                                               taACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGCTGAGCAGCCTCAGGGAGCAAATTTTAc                                                                                                      >  2:3274435/1‑149 (MQ=255)
                                                aaCGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACACCGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTAcc                                                                                                     <  2:3097886/149‑1 (MQ=255)
                                                        tcGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGTTTGATACTCCACATCAAGGTCAAACAACGGTAAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGt                                                                                             <  2:635886/148‑1 (MQ=255)
                                                          gACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTat                                                                                           >  1:314153/1‑149 (MQ=255)
                                                          gACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTat                                                                                           >  1:32332/1‑149 (MQ=255)
                                                          gACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTat                                                                                           >  1:2040913/1‑149 (MQ=255)
                                                           aCTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTatc                                                                                          >  1:2000722/1‑149 (MQ=255)
                                                               cccTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTAtcatca                                                                                      >  2:848544/1‑149 (MQ=255)
                                                               cccTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTAtcatca                                                                                      <  2:1501235/149‑1 (MQ=255)
                                                                ccTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAg                                                                                     >  1:3567413/1‑149 (MQ=255)
                                                                ccTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAg                                                                                     <  2:2903790/149‑1 (MQ=255)
                                                                 ctctAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGt                                                                                    >  2:1892284/1‑149 (MQ=255)
                                                                     aaCAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGtt                                                                                <  2:4595962/149‑1 (MQ=255)
                                                                      aCAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGttt                                                                               >  1:4339149/1‑149 (MQ=255)
                                                                      aCAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGttt                                                                               <  2:2746625/149‑1 (MQ=255)
                                                                       cAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTg                                                                              <  2:863728/149‑1 (MQ=255)
                                                                         aaGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGc                                                                            >  2:822150/1‑149 (MQ=255)
                                                                              aaGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGcc                                                                       >  1:1893854/1‑149 (MQ=255)
                                                                                gCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAg                                                                     <  1:328497/149‑1 (MQ=255)
                                                                                    cGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGa                                                                 <  1:1272999/149‑1 (MQ=255)
                                                                                     gCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGAt                                                                <  2:3660396/149‑1 (MQ=255)
                                                                                     gCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGAt                                                                >  1:2856856/1‑149 (MQ=255)
                                                                                      cTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATa                                                               >  1:3527722/1‑149 (MQ=255)
                                                                                        aTCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCACACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAc                                                             <  2:535901/149‑1 (MQ=255)
                                                                                        aTCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAc                                                             <  2:2319477/149‑1 (MQ=255)
                                                                                        aTCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAc                                                             >  2:3223328/1‑149 (MQ=255)
                                                                                        aTCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAc                                                             >  2:949913/1‑149 (MQ=255)
                                                                                           aaCTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACggg                                                          >  1:2887188/1‑149 (MQ=255)
                                                                                              tGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAg                                                       >  1:1705154/1‑149 (MQ=255)
                                                                                               gCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGa                                                      <  1:4222823/149‑1 (MQ=255)
                                                                                                cAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGaa                                                     >  1:3437337/1‑149 (MQ=255)
                                                                                                cAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGaa                                                     <  1:2535972/149‑1 (MQ=255)
                                                                                                  aGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAata                                                   >  1:3549209/1‑149 (MQ=255)
                                                                                                  aGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAata                                                   <  1:1063494/149‑1 (MQ=255)
                                                                                                   ggATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAatat                                                  <  2:3160417/149‑1 (MQ=255)
                                                                                                   ggATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAatat                                                  <  1:2402432/149‑1 (MQ=255)
                                                                                                    gATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATg                                                 <  1:140668/149‑1 (MQ=255)
                                                                                                     atatCAATGGATTGATACTCCACATCCAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGc                                                <  2:503379/149‑1 (MQ=255)
                                                                                                     atatCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGc                                                <  2:2626238/149‑1 (MQ=255)
                                                                                                     atatCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGc                                                <  1:112387/149‑1 (MQ=255)
                                                                                                     atatCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGc                                                <  1:968959/149‑1 (MQ=255)
                                                                                                      tatCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCg                                               >  1:2109197/1‑149 (MQ=255)
                                                                                                         cAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCt                                            <  2:1964987/149‑1 (MQ=255)
                                                                                                         cAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCt                                            <  2:1634204/149‑1 (MQ=255)
                                                                                                          aaTGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGGGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTa                                           >  2:2036059/1‑149 (MQ=255)
                                                                                                           aTGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAt                                          <  2:218053/149‑1 (MQ=255)
                                                                                                           aTGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAt                                          >  2:208254/1‑149 (MQ=255)
                                                                                                            tGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAtt                                         <  2:1073928/149‑1 (MQ=255)
                                                                                                             ggATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAttt                                        <  2:1598841/149‑1 (MQ=255)
                                                                                                             ggATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAttt                                        <  2:4068793/149‑1 (MQ=255)
                                                                                                              gATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTc                                       >  2:1392185/1‑149 (MQ=255)
                                                                                                                ttGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaa                                     >  2:4387131/1‑149 (MQ=255)
                                                                                                                ttGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaa                                     <  2:2256711/149‑1 (MQ=255)
                                                                                                                 tGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaac                                    >  1:4125918/1‑149 (MQ=255)
                                                                                                                   aTACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaacaa                                  >  2:3903793/1‑149 (MQ=255)
                                                                                                                      cTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATg                               <  1:4628273/149‑1 (MQ=255)
                                                                                                                        ccACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCt                             <  1:1504069/149‑1 (MQ=255)
                                                                                                                        ccACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCt                             <  1:3766342/149‑1 (MQ=255)
                                                                                                                         cacaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTg                            <  2:37433/149‑1 (MQ=255)
                                                                                                                         cacaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTg                            <  1:3940385/149‑1 (MQ=255)
                                                                                                                              cAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGtttt                       <  2:1662631/149‑1 (MQ=255)
                                                                                                                                 ggTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTg                    <  2:3266022/149‑1 (MQ=255)
                                                                                                                                   tCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACTGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGaa                  >  1:1639899/1‑149 (MQ=255)
                                                                                                                                   tCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGaa                  <  2:4437317/149‑1 (MQ=255)
                                                                                                                                   tCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGaa                  >  2:1484216/1‑149 (MQ=255)
                                                                                                                                      aacaacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGc               >  1:637470/1‑149 (MQ=255)
                                                                                                                                      aacaacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGc               <  2:2850875/149‑1 (MQ=255)
                                                                                                                                       acaacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCa              >  2:2660122/1‑149 (MQ=255)
                                                                                                                                       acaacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCa              <  1:299376/149‑1 (MQ=255)
                                                                                                                                         aacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATa            <  2:3354964/149‑1 (MQ=255)
                                                                                                                                         aacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATa            >  2:4295842/1‑149 (MQ=255)
                                                                                                                                          acGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATaa           >  1:552446/1‑149 (MQ=255)
                                                                                                                                          acGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATaa           <  2:3836220/149‑1 (MQ=255)
                                                                                                                                                aaTTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTCTTTCAACAAATGCTGGTTTTTTGAATGCATAACGACac     >  1:829538/1‑148 (MQ=255)
                                                                                                                                                aaTTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACaa     >  2:4028515/1‑149 (MQ=255)
                                                                                                                                                 aTTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACaaa    >  2:496071/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGA  >  CP009273/961798‑962093

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: