Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 725,276 A→T 10.0% S80C (AGT→TGT)  rhsC → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273725,2760AT10.0% 293.2 / 6.5 132S80C (AGT→TGT) rhsCRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base A (63/51);  new base T (2/15);  total (65/67)
Fisher's exact test for biased strand distribution p-value = 1.18e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

CCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTAT  >  CP009273/725129‑725404
                                                                                                                                                   |                                                                                                                                
cccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagca                                                                                                                                  >  2:3104277/1‑148 (MQ=37)
cccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagca                                                                                                                                  <  1:1705279/148‑1 (MQ=37)
cccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACATCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagca                                                                                                                                  <  2:972367/148‑1 (MQ=37)
  cACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGt                                                                                                                                <  1:2423298/148‑1 (MQ=37)
  cACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGt                                                                                                                                <  1:2423056/148‑1 (MQ=37)
         gTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGa                                                                                                                         <  1:2159541/148‑1 (MQ=37)
          tGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGAc                                                                                                                        >  2:2893685/1‑148 (MQ=37)
          tGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGAc                                                                                                                        >  2:2893715/1‑148 (MQ=37)
           ggCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACa                                                                                                                       <  1:299319/148‑1 (MQ=37)
             ccTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACaaa                                                                                                                     <  1:1702557/148‑1 (MQ=37)
              cTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACaaaa                                                                                                                    <  2:2329492/148‑1 (MQ=37)
              cTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACaaaa                                                                                                                    <  1:932905/148‑1 (MQ=37)
              cTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACaaaa                                                                                                                    <  1:2700843/148‑1 (MQ=37)
               tGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACaaaaa                                                                                                                   >  2:1389528/1‑148 (MQ=37)
               tGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACaaaaa                                                                                                                   >  2:1389507/1‑148 (MQ=37)
               tGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATACTCTCCCGCACCTACAGCAGTTACCGGACaaaaa                                                                                                                   >  2:1389758/1‑148 (MQ=21)
                 ttCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACg                                                                                                                 <  2:2860406/148‑1 (MQ=37)
                  tCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCATTTACCGGACAAAAACGc                                                                                                                <  2:3044779/148‑1 (MQ=37)
                  tCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACCAAAACGc                                                                                                                <  2:3044719/148‑1 (MQ=37)
                     gtgtGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTACAGCAGTTACCGGCCAAAAACGCCcg                                                                                                             >  1:2436282/1‑148 (MQ=255)
                     gtgtGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCcg                                                                                                             >  1:348333/1‑148 (MQ=37)
                     gtgtGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTGCCGTCTCTCCCCCACCTCCAGCAGGTACCGGACACAAACGCCcg                                                                                                             >  1:348218/1‑148 (MQ=11)
                     gtgtGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTGCATTCTCTCCCGCACCGACAGCAGTTACCCGACCACAACCCCcg                                                                                                             >  1:646620/1‑148 (MQ=11)
                          ccccGGCGGGGTGACGTCCGGCCATCGGGTAAATCCGGTGCTGGGTGCAAAGGTCCTTCCGGGTGAAACCGACATCGCCCTGCCGGGCCCGCTGCCGTTCTTTCTCTCCCGAACCTACCGCCGTTACCGGACAAAAACGCCCGCGCCg                                                                                                        <  2:2817237/148‑1 (MQ=11)
                                  gggTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAg                                                                                                >  2:222104/1‑148 (MQ=35)
                                         gTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTGCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGc                                                                                         >  1:1660608/1‑148 (MQ=21)
                                          tCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACCGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcc                                                                                        <  2:3242837/148‑1 (MQ=21)
                                          tCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCGGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcc                                                                                        <  2:2069578/148‑1 (MQ=21)
                                              gCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCgg                                                                                    >  1:2834116/1‑148 (MQ=34)
                                                actCCGGTCAATCCGCTGCTCGGTGCAACGGCCCTTCCCGGTGAAACCGACAGCGCCGTGCCCGGCCCGCTGCCGTTCAGTCTCTCCCCCACAGCCCGCTGTTCCCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCt                                                                                  <  2:2520042/146‑1 (MQ=11)
                                                 aTCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCGGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTg                                                                                 <  2:2082893/148‑1 (MQ=21)
                                                  tCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTgg                                                                                >  2:2577377/1‑148 (MQ=34)
                                                   ccGGTCAATCCGCTGCTGGGTGAAATGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCGGCGGTTCATTCTCGCCCGCACCTACCGCAGTTCCCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGa                                                                               <  1:3209342/148‑1 (MQ=11)
                                                     ggTCAATCCGCTGCTCGGTGCAAAGGTCCTGCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGaaa                                                                             <  2:3449114/148‑1 (MQ=21)
                                                        cAATCCGCTGCTCGGTGCAATGGTCCTTCCCGGTGAACCCGACATCGCCCGGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTCCCGCTGTTACCGCACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATg                                                                          <  2:664837/148‑1 (MQ=11)
                                                        cAAGCCGCTGCTCGGTGCAACGGTCCTTCCCGGTGAAACCGACAGCGCCCGGCCCGGCCCGCTGCCTTTCATTCACTCCCGCACCTACCGCTGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATg                                                                          <  2:667771/148‑1 (MQ=11)
                                                           tCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCt                                                                       >  1:2147253/1‑148 (MQ=34)
                                                             cgcgGCTCGGGGCAACGGTCCTGCCCGGTGAAACCGACAGCGCCCGGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCCACCGCAGTTCCCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGc                                                                     <  2:2730825/148‑1 (MQ=11)
                                                             cGCTGCGCGGTGCAACGGTCCTGCCCGGTGAAACCGACATCGCCCGGCCCGGCCCGGTGCCGTTCATTCTCTCCCGCACCTACCGCAGTTCCCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGc                                                                     <  2:2730826/148‑1 (MQ=11)
                                                              gctgctCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCg                                                                    >  1:714943/1‑148 (MQ=32)
                                                               gtgctCGGTGCAACGGGCCTTCCCGGTGAAACCGACAGCGCCCTGCCCGGCCCGCGGCCGTTCATTCTCTCCCCCACCTACCGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCgg                                                                   <  2:2550470/147‑1 (MQ=11)
                                                               ctgctCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCgg                                                                   >  1:2472604/1‑148 (MQ=32)
                                                               ctgctCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCgg                                                                   >  1:2472242/1‑148 (MQ=32)
                                                               ctgctCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCAACACCTGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCgg                                                                   <  2:2861427/148‑1 (MQ=12)
                                                                tgctCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACCGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGa                                                                  <  2:1139470/148‑1 (MQ=17)
                                                                        gCAAGGGTCCTTCCCGGTCAAACCGACACCGCCCGGCCCGGCCCGCTGCCGTACATTCTCACCCGCACCAACAGCAGTTCCCGCAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCt                                                          <  1:1742810/148‑1 (MQ=11)
                                                                          aaaGTTACTGCCGGGTCAAACGGAAATCGCCCGGCCCGGCCCGGTGCTGTTATTTATCTCCCGCACCTACTGATTTTACCGTAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTa                                                        <  1:2093543/148‑1 (MQ=2)
                                                                          aaaGGTCCTTCCCGGTGAAACCCACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTCCCGCACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTa                                                        <  1:2093538/148‑1 (MQ=12)
                                                                           aaGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTAc                                                       >  1:823817/1‑148 (MQ=32)
                                                                            aGGTCCTTCCCGGTGAAACCGACATCGCCCTGCACGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACa                                                      >  2:3295761/1‑148 (MQ=17)
                                                                              ggccTGCCCGGTGAACCCGACACCGCCCGCCCCGGCCCGCTGCCTTTCATTCGCGCCCGCACCTCCAGCAGTAACCGGAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGc                                                    <  1:1968766/146‑1 (MQ=11)
                                                                              gTCCTTCCCGGTTAAACCGCCATCGCCCGCCCCGGCCGGCTGCGGTAAATTCTCCCCCGCACCAACCGCAGTTCCCGCACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGc                                                    <  1:1968638/148‑1 (MQ=11)
                                                                                 cTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACCAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTgc                                                 >  1:3334319/1‑148 (MQ=17)
                                                                                    cccGTTGAAACCGACATCGCCCGGCCCGGCCCGCTGCCTTTATTTCTCTCCCCCCCCTCCCGCAGTTCCCGCAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGcgcg                                              <  2:1207764/148‑1 (MQ=11)
                                                                                    cccGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGcgcg                                              <  2:3464952/148‑1 (MQ=32)
                                                                                    cccGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGcgcg                                              >  1:2879377/1‑148 (MQ=32)
                                                                                      cGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGAt                                            >  2:712784/1‑148 (MQ=32)
                                                                                         tGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAac                                         >  1:1338880/1‑148 (MQ=32)
                                                                                          gAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAaca                                        >  2:1341971/1‑148 (MQ=32)
                                                                                            aaCCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCGTACAGCTGCGCGATAacaca                                      >  1:1076948/1‑148 (MQ=17)
                                                                                             aCCTACATCGCCCGGCCCGGCCCGCTGCCGTTCTTTCTCCCCCGCACCTACCGCTGTTACCGGAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAacacac                                     <  1:1179265/148‑1 (MQ=11)
                                                                                               cGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTACAGCAGTTACCGTACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTg                                   <  2:1021599/148‑1 (MQ=17)
                                                                                                gACATCGCCCTGCCCGGCCCGCTGCGGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGa                                  <  2:521170/148‑1 (MQ=17)
                                                                                                gACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGa                                  <  2:3494498/148‑1 (MQ=32)
                                                                                                 aCATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAAGGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGAt                                 >  2:1736683/1‑148 (MQ=17)
                                                                                                  caccGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATa                                <  2:926584/145‑1 (MQ=17)
                                                                                                  cATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATa                                >  1:3183841/1‑148 (MQ=32)
                                                                                                    tCGCCCTGCCGGGCCGGCTGCGGTTCATTCTCACCCGCACCTACAGCTGTTACCGGAAAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACt                              <  2:646620/148‑1 (MQ=11)
                                                                                                    tCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCATTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACt                              <  2:348333/148‑1 (MQ=17)
                                                                                                    gcGCCGTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTCCAGCAGTTCCCTGACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACt                              <  2:348218/147‑1 (MQ=11)
                                                                                                     cGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTc                             >  1:2934645/1‑148 (MQ=32)
                                                                                                       cccggcccggCCCGCTGCCTTTCTTTCTCTCCCGCACCTACAGCTGTTACCGAACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAg                           <  1:3529254/148‑1 (MQ=11)
                                                                                                        ccTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGt                          >  1:3526869/1‑148 (MQ=32)
                                                                                                         gggcccggccggCTGCGTTTCTTTCTCGCCCGCACATACACCTTTAACCGCAAAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTg                         <  1:2244256/146‑1 (MQ=11)
                                                                                                         cTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGCGCCGCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTg                         >  1:967885/1‑148 (MQ=17)
                                                                                                         cTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTg                         >  1:2705017/1‑148 (MQ=32)
                                                                                                         cTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTg                         >  2:2753432/1‑148 (MQ=32)
                                                                                                         cTGCCCGCCCCGCTGCGTTTCTTTCTCCCCCGCACCTACAGCTGTTACCGCAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTg                         <  1:764940/148‑1 (MQ=11)
                                                                                                          tGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGa                        >  1:1599858/1‑148 (MQ=32)
                                                                                                            cccggcccgCTGCCTTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATa                      <  2:3512827/148‑1 (MQ=17)
                                                                                                            cccggcccgCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTCCCGTACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATa                      <  2:1336967/148‑1 (MQ=17)
                                                                                                            cccggcccgCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATa                      <  2:3513032/148‑1 (MQ=32)
                                                                                                            cccggcccgCTGCCGTTCATTCTCTCCCCCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATa                      <  2:1818353/148‑1 (MQ=17)
                                                                                                              cggcccgCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAc                    >  2:2818256/1‑148 (MQ=32)
                                                                                                              cggcccgCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAc                    >  2:2817597/1‑148 (MQ=32)
                                                                                                                   cgCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAAAGCCCCCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcggcgg               >  1:2182095/1‑148 (MQ=17)
                                                                                                                       gCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGa           >  1:1195456/1‑148 (MQ=32)
                                                                                                                          gTTCATTCTCTCCCGCACCTACAGCTGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGc        >  1:607364/1‑148 (MQ=17)
                                                                                                                          gTTCATTCTCTCCCGCACCTACAGCTGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGc        >  1:607690/1‑148 (MQ=17)
                                                                                                                           ttCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGcc       >  2:2335356/1‑148 (MQ=32)
                                                                                                                           ttCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGcc       >  2:3425989/1‑148 (MQ=32)
                                                                                                                           ttCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGcc       >  2:695847/1‑148 (MQ=32)
                                                                                                                           ttCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGcc       >  2:695870/1‑148 (MQ=32)
                                                                                                                            tCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCt      >  2:3217582/1‑148 (MQ=32)
                                                                                                                              aTTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGt    >  1:2872218/1‑148 (MQ=32)
                                                                                                                              aTTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGt    >  2:451991/1‑148 (MQ=32)
                                                                                                                              aTTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGt    >  1:931652/1‑148 (MQ=32)
                                                                                                                                tctctcCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTAt  >  2:285290/1‑148 (MQ=11)
                                                                                                                                tctctcCCCCACCACCACCAGTTACCGCACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTAt  <  1:3482775/148‑1 (MQ=11)
                                                                                                                                tctcgcCCGCAAATCCAGCAGTTACCGGACAACAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTAt  <  1:2649697/148‑1 (MQ=11)
                                                                                                                                                   |                                                                                                                                
CCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTAT  >  CP009273/725129‑725404

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