Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,131,521 C→A 100% Q168K (CAG→AAG)  flgH → flagellar protein of basal‑body outer‑membrane L ring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,131,5210CA95.2% 618.4 / 20.6 209Q168K (CAG→AAG) flgHflagellar protein of basal‑body outer‑membrane L ring
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base A (95/103);  total (99/110)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 215 total aligned reads displayed.

ACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAACGAAGC  >  CP009273/1131375‑1131666
                                                                                                                                                  |                                                                                                                                                 
  gCTCGTGCCGATGTCGACGCCTCCGGTGGTAACACGTAAAACGGAACGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTCCGGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGAt                                                                                                                                               <  1:1765495/149‑1 (MQ=255)
  gCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGAt                                                                                                                                               >  1:3904955/1‑149 (MQ=255)
  gCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCAGGTGGTGGGTGAAAAAAAGAt                                                                                                                                               >  1:894022/1‑149 (MQ=255)
  gCTCGTGCCGATGTCCAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTCGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGAt                                                                                                                                               <  1:703292/149‑1 (MQ=255)
     cGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGc                                                                                                                                            <  2:1229541/149‑1 (MQ=255)
     cGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGc                                                                                                                                            <  2:2087488/149‑1 (MQ=255)
       tGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCa                                                                                                                                          <  2:1562859/149‑1 (MQ=255)
        gCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGACAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGAGTGCCAt                                                                                                                                         <  2:1990861/149‑1 (MQ=255)
         cggATGTCGAAGCCTCCGGTGGTAACACTTTCAACGGAACGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTTACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTTCATGTGGTTGGTGAAAAAAAGATTGCCAtt                                                                                                                                        <  1:2317921/147‑1 (MQ=255)
         ccGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCAtt                                                                                                                                        <  1:1350064/149‑1 (MQ=255)
         ccGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCAtt                                                                                                                                        <  2:4102456/149‑1 (MQ=255)
         ccGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCAtt                                                                                                                                        <  1:3408460/149‑1 (MQ=255)
          cGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTa                                                                                                                                       >  2:924538/1‑149 (MQ=255)
           gATGTCGAAGCCTCCGGTGTTAACACGTTCAACGGAACGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTaa                                                                                                                                      <  1:188663/149‑1 (MQ=255)
           gATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTaa                                                                                                                                      <  2:3256418/149‑1 (MQ=255)
              gTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCa                                                                                                                                   <  2:2886606/149‑1 (MQ=255)
               tCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAg                                                                                                                                  >  1:1061876/1‑149 (MQ=255)
                 gAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAggg                                                                                                                                <  2:1221154/149‑1 (MQ=255)
                 gAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAggg                                                                                                                                >  2:3821200/1‑149 (MQ=255)
                  aaGCCTCCGGTGTTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGt                                                                                                                               <  1:2614901/149‑1 (MQ=255)
                   aGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTa                                                                                                                              >  2:437834/1‑149 (MQ=255)
                         cGGTGGTAACACGTTCACCGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAAt                                                                                                                        <  2:3862478/149‑1 (MQ=255)
                          ggtggtAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTCGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAAtt                                                                                                                       <  2:3797329/149‑1 (MQ=255)
                           gtggtAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAAttt                                                                                                                      >  1:188524/1‑149 (MQ=255)
                             ggtAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCACCCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTAt                                                                                                                    >  1:3658476/1‑149 (MQ=255)
                              gtAACACGTTCAACGTAAAGGGCGGGGCAAATGCCAGAAATACATTTCGCGGCACGTTGAGGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTAtt                                                                                                                   <  1:24864/149‑1 (MQ=255)
                              gtAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTAtt                                                                                                                   >  1:1642279/1‑149 (MQ=255)
                               tAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTc                                                                                                                  >  1:3796057/1‑149 (MQ=255)
                                  cacGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCt                                                                                                               >  2:3903633/1‑149 (MQ=255)
                                  cacGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCt                                                                                                               >  1:1727788/1‑149 (MQ=255)
                                   acGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCtt                                                                                                              >  1:3351724/1‑149 (MQ=255)
                                   acGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCtt                                                                                                              >  2:1930281/1‑149 (MQ=255)
                                   acGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCtt                                                                                                              <  1:3941456/149‑1 (MQ=255)
                                      ttCATCGGAAAGGGCGGGGCCAATGCGAGCAATACCTTTAGCGGCACGTTGACGTTGACTTTTGACCAGGTTCTGGTAAACGGCAAACTGCTTGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGTGTACAGAATTTATTCGCTtctc                                                                                                           <  2:3559216/149‑1 (MQ=255)
                                      ttCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTtctc                                                                                                           <  2:2407834/149‑1 (MQ=255)
                                      ttCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTtctc                                                                                                           <  2:3226500/149‑1 (MQ=255)
                                           cGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCg                                                                                                      >  2:2144528/1‑149 (MQ=255)
                                             gAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTg                                                                                                    >  1:430755/1‑149 (MQ=255)
                                             gAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTg                                                                                                    >  1:1791067/1‑149 (MQ=255)
                                             gAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTg                                                                                                    >  1:198479/1‑149 (MQ=255)
                                             gAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACATGGTAATGGTCATCGGCAACCTGCATGTGGTGGGTGTAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTg                                                                                                    <  1:2077296/149‑1 (MQ=255)
                                              aaaGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTgg                                                                                                   <  2:66741/149‑1 (MQ=255)
                                                       ggCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCAc                                                                                          <  2:3286373/149‑1 (MQ=255)
                                                       ggCCAATGCCAGCAAGACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCAc                                                                                          <  2:2400474/149‑1 (MQ=255)
                                                          cactGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCa                                                                                       <  2:2104660/146‑1 (MQ=255)
                                                           aaTGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCAc                                                                                      >  1:2552137/1‑149 (MQ=255)
                                                             tGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTa                                                                                    <  1:3503700/149‑1 (MQ=255)
                                                               ccAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATc                                                                                  <  2:1424383/149‑1 (MQ=255)
                                                               ccAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATc                                                                                  <  2:3870993/149‑1 (MQ=255)
                                                                cAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTGCTCTGGCGTGGTTAATCCACGCACTATCa                                                                                 >  1:3920275/1‑149 (MQ=255)
                                                                 aGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAg                                                                                >  1:514519/1‑149 (MQ=255)
                                                                 aGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAg                                                                                >  1:2236050/1‑149 (MQ=255)
                                                                 aGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAg                                                                                <  1:3630255/149‑1 (MQ=255)
                                                                 aGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAg                                                                                <  1:3630259/149‑1 (MQ=255)
                                                                     aTACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgc                                                                            >  2:1070380/1‑149 (MQ=255)
                                                                             aGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtac                                                                    <  1:548206/149‑1 (MQ=255)
                                                                              gCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtacc                                                                   >  1:1310149/1‑149 (MQ=255)
                                                                                ggCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccgt                                                                 >  1:2001298/1‑149 (MQ=255)
                                                                                   aCGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccgtacc                                                              >  2:2211162/1‑149 (MQ=255)
                                                                                      ttGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTc                                                           >  1:505002/1‑149 (MQ=255)
                                                                                        gacggtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTa                                                         <  1:2047615/149‑1 (MQ=255)
                                                                                         acggtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTAc                                                        >  1:1262483/1‑149 (MQ=255)
                                                                                         acggtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTAc                                                        <  1:3416741/149‑1 (MQ=255)
                                                                                         acggtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTAc                                                        <  1:4098958/149‑1 (MQ=255)
                                                                                            gtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCa                                                     >  2:138821/1‑149 (MQ=255)
                                                                                            gtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCa                                                     >  2:2607399/1‑149 (MQ=255)
                                                                                            gtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCa                                                     <  2:3689564/149‑1 (MQ=255)
                                                                                            gtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCa                                                     >  1:3999529/1‑149 (MQ=255)
                                                                                            gtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCa                                                     <  1:3753868/149‑1 (MQ=255)
                                                                                            gtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCa                                                     <  1:298752/149‑1 (MQ=255)
                                                                                            gtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCa                                                     <  1:1708687/149‑1 (MQ=255)
                                                                                                cggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTg                                                 <  2:3415671/149‑1 (MQ=255)
                                                                                                  gtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGc                                               >  2:3102577/1‑149 (MQ=255)
                                                                                                     gACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGa                                            >  1:188970/1‑149 (MQ=255)
                                                                                                      aCCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGAt                                           <  2:3730058/149‑1 (MQ=255)
                                                                                                      aCCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGAt                                           <  2:1685939/149‑1 (MQ=255)
                                                                                                               tGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCAtt                                  >  2:2467292/1‑149 (MQ=255)
                                                                                                                   cAACGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAAt                              >  2:3285752/1‑149 (MQ=255)
                                                                                                                    aaCGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATa                             <  1:1511800/149‑1 (MQ=255)
                                                                                                                    aaCGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATa                             <  1:1417883/149‑1 (MQ=255)
                                                                                                                      cGGCAACCTGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACg                           <  1:980931/149‑1 (MQ=255)
                                                                                                                        gCAACCTGCATGTGGTGGGTGAAAAACAGAGTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTa                         <  2:1185753/149‑1 (MQ=255)
                                                                                                                              tGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAAt                   >  2:3163847/1‑149 (MQ=255)
                                                                                                                              tGCATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAAt                   >  1:480124/1‑149 (MQ=255)
                                                                                                                                cATGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATgg                 <  1:1049941/149‑1 (MQ=255)
                                                                                                                                 aTGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGc                >  2:437452/1‑149 (MQ=255)
                                                                                                                                 aTGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGc                <  1:3871553/149‑1 (MQ=255)
                                                                                                                                 aTGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGc                >  1:1540006/1‑149 (MQ=255)
                                                                                                                                 aTGTGGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGc                >  2:2764832/1‑149 (MQ=255)
                                                                                                                                  tgtgGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCt               <  1:3821200/149‑1 (MQ=255)
                                                                                                                                  tgtgGTGGGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCt               <  2:2533224/149‑1 (MQ=255)
                                                                                                                                    tggtggGTGAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTAc             <  2:198479/149‑1 (MQ=255)
                                                                                                                                      gtgggtgAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACAt           <  2:3795792/149‑1 (MQ=255)
                                                                                                                                        gggtgAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTa         >  1:3612383/1‑149 (MQ=255)
                                                                                                                                          gtgAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAAc       <  1:1726239/149‑1 (MQ=255)
                                                                                                                                          gtgAAAAAAAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAAc       <  1:343364/149‑1 (MQ=255)
                                                                                                                                              aaaaaaGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAACGAAg   <  1:3447969/149‑1 (MQ=255)
                                                                                                                                              aaaaaaGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAACGAAg   <  2:1633419/149‑1 (MQ=255)
                                                                                                                                               aaaaAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACTCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAACGAAgc  >  2:690623/1‑149 (MQ=255)
                                                                                                                                               aaaaAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAACGAAgc  <  2:2965606/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCATTGAATACGTAGGCAATGGCTACATTAACGAAGC  >  CP009273/1131375‑1131666

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