Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,687,673:1 +G 5.6% coding (1349/2061 nt) pilJ ← twitching motility protein PilJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,687,6731.G 5.6% 438.0 / 18.5 107coding (1349/2061 nt)pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base . (52/49);  new base G (2/4);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

Only 100 of 120 total aligned reads displayed.

TCCGCTTGGCGGTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGACC  >  NC_002947/5687527‑5687822
                                                                                                                                                   |                                                                                                                                                     
tCCGCTTGGCGGTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGtcgatcg                                                                                                                                                    >  1:1161651/1‑150 (MQ=255)
 ccGCTTGGCGGTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGtcgatcga                                                                                                                                                   <  1:1533098/150‑1 (MQ=255)
    cTTGGCGGTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAtcgactc                                                                                                                                                <  1:1888814/150‑1 (MQ=255)
    cTTGGCGGTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAtcgactc                                                                                                                                                <  1:1968670/150‑1 (MQ=255)
          ggTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCAt                                                                                                                                          <  1:1235837/150‑1 (MQ=255)
           gtgtCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTTGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAGTCGATCGAccggg                                                                                                                                         <  2:2179212/150‑4 (MQ=255)
                  ggATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTcaccca                                                                                                                                  >  1:1847316/1‑150 (MQ=255)
                       tGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGcc                                                                                                                             <  2:1527140/150‑1 (MQ=255)
                       tGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGcc                                                                                                                             >  2:834208/1‑150 (MQ=255)
                         cTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTc                                                                                                                           >  2:1194302/1‑150 (MQ=255)
                            gcgAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGa                                                                                                                        <  1:2019109/150‑1 (MQ=255)
                            gcgAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGa                                                                                                                        <  2:1920060/150‑1 (MQ=255)
                              gAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAgg                                                                                                                      <  2:1628871/150‑1 (MQ=255)
                              gAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAgg                                                                                                                      <  2:917203/150‑1 (MQ=255)
                                 tgttgtCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCt                                                                                                                   >  1:1215449/1‑150 (MQ=255)
                                      tCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTg                                                                                                              >  1:328420/1‑150 (MQ=255)
                                      tCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTg                                                                                                              >  2:1079130/1‑150 (MQ=255)
                                         aTGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATc                                                                                                           <  2:684342/150‑1 (MQ=255)
                                             cGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCtgtg                                                                                                       <  1:103960/150‑1 (MQ=255)
                                                 gATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGc                                                                                                   <  1:1079130/150‑1 (MQ=255)
                                                 gATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGc                                                                                                   >  1:915660/1‑150 (MQ=255)
                                                  aTGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAGTCGATCGAccgg                                                                                                                                          <  2:1756378/110‑3 (MQ=255)
                                                  aTGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAGTCGATCGAccgg                                                                                                                                          >  1:1756378/1‑108 (MQ=255)
                                                      tgttgtGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGt                                                                                              >  2:2345303/1‑150 (MQ=255)
                                                        ttgtGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGc                                                                                            >  2:385366/1‑150 (MQ=255)
                                                             caccacCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGa                                                                                       <  2:382055/150‑1 (MQ=255)
                                                              accaccTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCCACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAg                                                                                      <  1:1500205/150‑1 (MQ=255)
                                                               ccaccTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCGTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAgg                                                                                     >  2:1481072/1‑150 (MQ=255)
                                                                caccTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGc                                                                                    >  2:1734582/1‑150 (MQ=255)
                                                                caccTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCCTTCCAGTTTCGCGACCATGTCACCCACGGCCTAGGAGGCACAGCTGATCTGTGCAGCCTGGTGCTCGGAGGc                                                                                    >  1:1853962/1‑150 (MQ=255)
                                                                  ccTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCt                                                                                  >  2:1202814/1‑150 (MQ=255)
                                                                   cTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCtt                                                                                 >  1:684893/1‑150 (MQ=255)
                                                                   cTCGTTGCCCTTGTTGGCAATGGCCACCGAACGATGGGGGAATCTGGACGACGGAGATGCCTTCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCtt                                                                                 >  1:46360/1‑150 (MQ=255)
                                                                    tCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTg                                                                                >  1:685971/1‑150 (MQ=255)
                                                                     cGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCGCCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGTTGCTCGGAGGCCTTgg                                                                               <  2:215918/150‑1 (MQ=255)
                                                                      gtcgCCCTTGTTGGCAATGGCCACCGAACGCGCGGCCACGTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCCACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGc                                                                              <  1:508322/147‑1 (MQ=255)
                                                                      gTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGc                                                                              <  2:1847316/150‑1 (MQ=255)
                                                                        tGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCa                                                                            <  2:1215449/150‑1 (MQ=255)
                                                                           ccTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCa                                                                                                                                           <  2:1046127/84‑1 (MQ=255)
                                                                           ccTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCa                                                                                                                                           >  1:1046127/1‑84 (MQ=255)
                                                                           ccTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGtt                                                                         <  2:107811/150‑1 (MQ=255)
                                                                           ccTTGTTAGAAATGGCCACCGAACGCTTGGCCACCGTGGCCGATTCGGAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCCCCCAACGCCGAAGAGGCCGCACGCATCGACTCAGCCGGGTGCTCGGATGCCTTGGCCAGtt                                                                         >  2:1728124/1‑150 (MQ=255)
                                                                            cTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGtt                                                                         <  2:1079683/149‑1 (MQ=255)
                                                                            cTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTg                                                                        >  2:230591/1‑150 (MQ=255)
                                                                                ttctcAATGGCCACGGAACGCTCGGTCACCTTGGCCGGTTCGTAGACATGCGCCGAAACCCGGTTGATCGCCCCCACCATGTCACCAACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGGGCTCGGAGGCCTTGGCCAGTTGCcgc                                                                    >  2:970984/5‑150 (MQ=255)
                                                                                     aaTGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAgt                                                               >  1:1176367/1‑150 (MQ=255)
                                                                                     aaTGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAgt                                                               <  2:685971/150‑1 (MQ=255)
                                                                                           cACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAANCCGGTCGGTCGACTCGACCATGTCACCCACCGCCTCGGGGGCCTCGCTGATCTGTGCAGCCTGGAGCTCGGGGGCCTTGGCCAGTTGCCGCGCAGTGTTCTg                                                         >  1:1037830/1‑150 (MQ=255)
                                                                                           cACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGtcgatc                                                                                                                                                     >  1:1075132/1‑58 (MQ=255)
                                                                                           cACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGtcgatc                                                                                                                                                     <  2:1075132/58‑1 (MQ=255)
                                                                                           cACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAGTCGATCGAccgg                                                                                                                                          >  1:1095354/1‑67 (MQ=38)
                                                                                           cACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAGTCGATCGAccgg                                                                                                                                          <  2:1095354/69‑3 (MQ=38)
                                                                                           cACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGGAGCCTGGTGCTCGGAGGCCTTGGCCAATTGTCGCGCAGTGTTCTg                                                         >  2:1640401/1‑150 (MQ=255)
                                                                                            aCCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCAt                                                                                                                                          >  1:408557/1‑68 (MQ=255)
                                                                                            aCCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCAt                                                                                                                                          <  2:408557/68‑1 (MQ=255)
                                                                                            aCCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGTCGTGGCCAGTTGCCGCGCAGTGTTCTGc                                                        <  1:1872927/150‑1 (MQ=255)
                                                                                            aCCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGc                                                        >  1:1501238/1‑150 (MQ=255)
                                                                                            aCCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGc                                                        <  1:1194302/150‑1 (MQ=255)
                                                                                              cGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCCACACCCGGTCGAACTGCTCGACCATGTCACCACCAATCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCgt                                                      <  1:2342927/150‑1 (MQ=255)
                                                                                                 aCGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTc                                                   <  1:1734582/150‑1 (MQ=255)
                                                                                                 aCGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTc                                                   <  1:204132/150‑1 (MQ=255)
                                                                                                   gCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGACCTGGCTGATCTGTGGAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCt                                                 >  1:151922/1‑150 (MQ=255)
                                                                                                     tCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGc                                                                                                                <  1:1360264/85‑1 (MQ=255)
                                                                                                       ggCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTGGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCAc                                             >  1:1147834/1‑150 (MQ=255)
                                                                                                       ggCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGcc                                                                                                  >  2:254308/1‑97 (MQ=255)
                                                                                                       ggCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGcc                                                                                                  <  1:254308/97‑1 (MQ=255)
                                                                                                        gCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCAcg                                            <  1:1202814/150‑1 (MQ=255)
                                                                                                          cACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGTAGCCTGGTGCTCgg                                                                                        <  2:291453/104‑1 (MQ=255)
                                                                                                          cACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGTAGCCTGGTGCTCgg                                                                                        >  1:291453/1‑104 (MQ=255)
                                                                                                             cTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCtgtg                                                                                                       >  1:2194229/1‑86 (MQ=255)
                                                                                                             cTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCtgtg                                                                                                       <  2:2194229/86‑1 (MQ=255)
                                                                                                             cTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCa                                                                                                     <  1:604770/88‑1 (MQ=255)
                                                                                                             cTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCa                                                                                                     >  2:604770/1‑88 (MQ=255)
                                                                                                             cTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCAcggcggc                                       >  2:1266739/1‑150 (MQ=255)
                                                                                                                 gCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGgcggcggcg                                   >  1:1736766/1‑150 (MQ=255)
                                                                                                                  ccGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCTTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGa                                  >  2:1579763/1‑150 (MQ=255)
                                                                                                                   cGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGAc                                 >  1:455912/1‑150 (MQ=255)
                                                                                                                          tAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACctgcact                          <  2:1777033/150‑1 (MQ=255)
                                                                                                                             gCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGa                                  <  2:1683994/139‑1 (MQ=255)
                                                                                                                             gCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACctgcactgca                       <  2:1853962/150‑1 (MQ=255)
                                                                                                                               aTGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGAGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTgcactgcact                     >  2:734328/1‑150 (MQ=255)
                                                                                                                                  cgcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACtgtg                  <  1:2345303/150‑1 (MQ=255)
                                                                                                                                  cgcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACtgtg                  <  2:1176367/150‑1 (MQ=255)
                                                                                                                                   gcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCt                                                                                                               <  1:433118/56‑1 (MQ=255)
                                                                                                                                   gcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCt                                                                                                               >  2:433118/1‑56 (MQ=255)
                                                                                                                                   gcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTgg                 >  1:2005777/1‑150 (MQ=255)
                                                                                                                                   gcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTgg                 >  1:688384/1‑150 (MQ=255)
                                                                                                                                    ccGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGt                <  1:230591/150‑1 (MQ=255)
                                                                                                                                    ccGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGt                <  2:1501238/150‑1 (MQ=255)
                                                                                                                                     cGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGTCCAGTTTGTGCGCATTGTTCGTCGTGTCCTTCACTTTGACGGCGACCTGCACTGCACTGTGGtt               >  2:1109330/1‑150 (MQ=255)
                                                                                                                                     cGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACctgcact                          <  2:1428274/139‑1 (MQ=255)
                                                                                                                                     cGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACctgcact                          >  1:1428274/1‑139 (MQ=255)
                                                                                                                                        aaCCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGTAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGAt            <  1:1109330/150‑1 (MQ=255)
                                                                                                                                          cccGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGAtgg          <  1:1748197/150‑1 (MQ=255)
                                                                                                                                             ggTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCCGTGTTCTGCGTGTCCTACAGCGCGTGGGCGACCGTCACTGCACGGTGGTTGAtggtgg       <  2:1835452/150‑1 (MQ=255)
                                                                                                                                                cgatcgaCTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGa    <  1:385366/150‑1 (MQ=255)
                                                                                                                                                 gatcgaCTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACtgtg                  >  1:861818/1‑135 (MQ=255)
                                                                                                                                                 gatcgaCTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACtgtg                  <  2:861818/135‑1 (MQ=255)
                                                                                                                                                 gatcgaCTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACCGCACTGTGGTTGATGGTGGCGAc   >  1:815960/1‑150 (MQ=255)
                                                                                                                                                  atcgaCTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGAcc  >  1:1743869/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCCGCTTGGCGGTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGACC  >  NC_002947/5687527‑5687822

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