Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,980,399 T→G 5.2% intergenic (‑145/+16) flgF ← / ← PP_4387 flagellar basal‑body rod protein FlgF/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,980,3990TG 5.2% 230.1 / 2.8 76intergenic (‑145/+16)flgF/PP_4387flagellar basal‑body rod protein FlgF/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base G (2/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TCCACAGTTTCTCCTCTGCGCGTAGCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGG  >  NC_002947/4980250‑4980548
                                                                                                                                                     |                                                                                                                                                     
tCCACAGTTTCTCCTCTGCGCGTAGCTGCCGCTCTGTCACTGGTCACGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGCCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCt                                                                                                                                                       <  2:1279027/150‑1 (MQ=255)
 ccACAGTTTCTCCTCTGCGCGTAGCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTg                                                                                                                                                      >  1:704795/1‑150 (MQ=255)
  cacaGTTTCTCCTCTGCGCGTAGCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGc                                                                                                                                                     >  2:577153/1‑150 (MQ=255)
   acaGTTTCTCCTCTGCGCGTAGCTTCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCg                                                                                                                                                    >  1:1037805/1‑150 (MQ=255)
   acaGTTTCTCCTCTGCGCGTAGCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCg                                                                                                                                                    >  2:1065137/1‑150 (MQ=255)
        ttCTCCTCTGCGCGTAGCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAg                                                                                                                                               <  2:907023/150‑1 (MQ=255)
        ttCTCATCTTCGCGTAGCTGCCGCGCTGTAACTAGCCTATTACCAATTAATGATTGTACCAAGTCTGAAACAGTCAGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAg                                                                                                                                               <  1:672950/150‑1 (MQ=255)
         tctcCTCTGCGCGTAGCTGCGGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAgg                                                                                                                                              <  2:1037805/150‑1 (MQ=255)
             ctctGCGCGTAGCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCtt                                                                                                                                          <  2:606061/150‑1 (MQ=255)
             ctctGCGCGGAGCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCtt                                                                                                                                          <  1:421770/150‑1 (MQ=255)
                    cgTAGCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAAt                                                                                                                                   <  2:1197521/150‑1 (MQ=255)
                    cgTAGCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAAt                                                                                                                                   <  1:577153/150‑1 (MQ=255)
                        gCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGgc                                                                                                                               <  2:334016/150‑1 (MQ=255)
                             cgcgCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCCTGCTCAATGGGCGCGGc                                                                                                                          >  1:464452/1‑150 (MQ=255)
                             cgagcTGTCACTGGGAAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGc                                                                                                                          >  2:524367/4‑150 (MQ=255)
                              gcgcTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAAt                                                                                                                                   <  2:1135051/140‑1 (MQ=255)
                              gcgcTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAAt                                                                                                                                   >  1:1135051/1‑140 (MQ=255)
                                  tGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAggcg                                                                                                                     >  1:1378652/1‑150 (MQ=255)
                                  tGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAggcg                                                                                                                     >  2:1395598/1‑150 (MQ=255)
                                  tGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACAGTAAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAggcg                                                                                                                     >  2:1260583/1‑150 (MQ=255)
                                   gTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAggcgg                                                                                                                    <  1:1296474/150‑1 (MQ=255)
                                   gTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAggcgg                                                                                                                    >  2:120397/1‑150 (MQ=255)
                                    tCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAggcggc                                                                                                                   >  1:650523/1‑150 (MQ=255)
                                     cACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGgcggcg                                                                                                                  >  2:315895/1‑150 (MQ=255)
                                      aCTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCgt                                                                                                                 >  2:212588/1‑150 (MQ=255)
                                          ggCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCgttgtt                                                                                                             <  1:212588/150‑1 (MQ=255)
                                             aGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTa                                                                                                          >  2:719574/1‑150 (MQ=255)
                                              ggTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAg                                                                                                         <  1:1299904/150‑1 (MQ=255)
                                               gTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTTCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAgg                                                                                                        >  2:214788/1‑150 (MQ=255)
                                               gTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAgg                                                                                                        <  2:384500/150‑1 (MQ=255)
                                                 tCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAgggg                                                                                                      >  2:1098304/1‑150 (MQ=255)
                                                   cAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGcc                                                                                                    >  1:423724/1‑150 (MQ=255)
                                                   cAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGcc                                                                                                    >  2:1120920/1‑150 (MQ=255)
                                                    aaGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCa                                                                                                   >  2:386188/1‑150 (MQ=255)
                                                      gCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCaga                                                                                                 >  1:583089/1‑150 (MQ=255)
                                                       cAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCagag                                                                                                >  1:551963/1‑150 (MQ=255)
                                                           aatTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGgc                                                                                            <  1:1040253/150‑1 (MQ=255)
                                                           aatTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGgc                                                                                            <  1:547440/150‑1 (MQ=255)
                                                               gAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCAc                                                                                        >  2:787149/1‑150 (MQ=255)
                                                                   cAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGc                                                                                                     <  2:1057037/133‑1 (MQ=255)
                                                                   cAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGc                                                                                                     >  1:1057037/1‑133 (MQ=255)
                                                                    aaGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGc                                                                                   <  2:603258/150‑1 (MQ=255)
                                                                        cTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCttt                                                                               <  1:1098304/150‑1 (MQ=255)
                                                                              cacCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCgg                                                                         <  1:719574/150‑1 (MQ=255)
                                                                              cacCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCgg                                                                         <  2:574074/150‑1 (MQ=255)
                                                                               acCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCggg                                                                        <  1:787149/150‑1 (MQ=255)
                                                                                 ccGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCGGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGtt                                                                      >  2:173323/1‑150 (MQ=255)
                                                                                  cGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTtg                                                                     <  1:1395598/150‑1 (MQ=255)
                                                                                      ttACGGGGCTTTCAGCGCGGTAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTCGTCGTCGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTggg                                                                 <  1:889881/150‑1 (MQ=255)
                                                                                        aCGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGcc                                                               <  1:304651/150‑1 (MQ=255)
                                                                                             gCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAaca                                                          <  1:1120920/150‑1 (MQ=255)
                                                                                               tttCAGCGCGGCAAGTGGCAATAACCGTCAGTTTGTCGGCTGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGc                                                                <  1:401204/142‑1 (MQ=255)
                                                                                               tttCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGc                                                                >  2:401204/1‑142 (MQ=255)
                                                                                                     cgcgGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGc                                                  <  2:583089/150‑1 (MQ=255)
                                                                                                          cAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCg                                             <  1:173323/150‑1 (MQ=255)
                                                                                                           aaGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATTGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGAGCACCCCGCCTTCGACCGGGTCATGGGCCGAACACAGGTAGCGGCCgg                                            <  2:1378652/150‑1 (MQ=255)
                                                                                                           aaGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCgg                                            <  1:424667/150‑1 (MQ=255)
                                                                                                                gCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGAcc                                       >  2:1343091/1‑150 (MQ=255)
                                                                                                                              tttGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCGATTg                         >  1:21668/1‑150 (MQ=255)
                                                                                                                               ttGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGa                        <  2:650523/150‑1 (MQ=255)
                                                                                                                                    ggcggATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGaa                   >  1:114323/1‑150 (MQ=255)
                                                                                                                                        gATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCAtt                          >  2:911904/1‑139 (MQ=255)
                                                                                                                                        gATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCAtt                          <  1:911904/139‑1 (MQ=255)
                                                                                                                                        gATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCAt                           <  1:820355/138‑1 (MQ=255)
                                                                                                                                        gATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCAt                           >  2:820355/1‑138 (MQ=255)
                                                                                                                                        gATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAAcggg               >  1:354714/1‑150 (MQ=255)
                                                                                                                                          tGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCgg                                                >  1:893153/1‑115 (MQ=255)
                                                                                                                                          tGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCgg                                                <  2:893153/115‑1 (MQ=255)
                                                                                                                                          tGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAAcggggc             >  2:1199463/1‑150 (MQ=255)
                                                                                                                                           ccttGCCGCAGGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATa                      >  1:1351636/3‑140 (MQ=255)
                                                                                                                                           ccttGCCGCAGGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATa                      <  2:1351636/138‑1 (MQ=255)
                                                                                                                                           ccttGCCGCAGGCGGCAAGGGCTTGCTCAACGGGCGCGCAAACCGGGGTCAACCATGGGTACAGAGTGTGCCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAAcggggcg            >  1:739330/3‑150 (MQ=255)
                                                                                                                                            tttGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAAcggggcgg           <  2:114323/150‑1 (MQ=255)
                                                                                                                                             ttGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAAcggggcggg          <  2:423724/150‑1 (MQ=255)
                                                                                                                                             ttGCCGCAGGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAAcggggcggg          <  2:739330/150‑1 (MQ=255)
                                                                                                                                                   ccTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGgcgggggcagc    >  2:77729/1‑150 (MQ=255)
                                                                                                                                                   ccTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGgcgggggcagc    <  2:551963/150‑1 (MQ=255)
                                                                                                                                                     tGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCgg                                                                         <  1:883664/79‑1 (MQ=255)
                                                                                                                                                     tGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCgg                                                                         >  2:883664/1‑79 (MQ=255)
                                                                                                                                                     tGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGgcgggggcagcgg  <  2:21668/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCACAGTTTCTCCTCTGCGCGTAGCTGCCGCGCTGTCACTGGGCAGGTTCCAAGCAATAATTGAACCAAGTCTGAAACACCCGGTTTACGGGGCTTTCAGCGCGGCAAGGGGCAATAACCGTCAGTTTGTCGGCGGATGTTTGCCGCCTGCGGCAAGGGCTTGCTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGG  >  NC_002947/4980250‑4980548

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