Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,333,668 T→G 7.2% I33S (ATC→AGC)  PP_4695 → two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,333,6680TG 7.2% 207.7 / 6.5 70I33S (ATC→AGC) PP_4695two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base G (3/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TGGCTGATCCTGATCGTCGTCACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTGGGGGTGTCTGGC  >  NC_002947/5333525‑5333810
                                                                                                                                               |                                                                                                                                              
tGGCTGATCCTGATCGTCGTCACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCCCGCGCTCGCTGATcca                                                                                                                                             <  7:118924/147‑2 (MQ=255)
  gCTGATCCTGATCGTCGTCACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGt                                                                                                                                         >  6:253047/1‑149 (MQ=255)
  gCTGATCCTGATCGTCGTCACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGt                                                                                                                                         >  7:272641/1‑149 (MQ=255)
  gCTGATCCTGATCGTCGTCACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATAATCGCCGGGTACCTAATTGTGCTGTTCGGCCTGGCCTGGATCAGCGAGCGCGGTCTCATCCCGCGt                                                                                                                                         >  4:21542/1‑149 (MQ=255)
       tccagatcTTCGTCACCTACATCCTCCAGAGCAAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGAt                                                                                                                                    >  4:92980/1‑149 (MQ=255)
        cctgatcGTCGTCACCTTCATCGTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGAt                                                                                                                                    >  7:61133/1‑148 (MQ=255)
        cctgatcGTCGTCACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAACGCGGGATCAGCCCGCGCTCgctg                                                                                                                                   <  5:117714/149‑4 (MQ=255)
          tgatcGTCGTCACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCAt                                                                                                                                 >  4:381291/1‑149 (MQ=255)
            atcGTCGTCACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCAtt                                                                                                                                <  3:363225/148‑1 (MQ=255)
             tcgtcgTCACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCg                                                                                                                              <  3:21542/149‑1 (MQ=255)
               gtccgcACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGc                                                                                                                             <  8:156138/143‑1 (MQ=255)
                   tcACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGccac                                                                                                                          >  6:334921/1‑147 (MQ=255)
                   tcACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGccac                                                                                                                          <  5:334921/147‑1 (MQ=255)
                     aCCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGccacccac                                                                                                                      <  5:96640/149‑1 (MQ=255)
                     aCCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGccaccca                                                                                                                       <  8:272641/148‑1 (MQ=255)
                            tCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTa                                                                                                               >  7:227615/1‑149 (MQ=255)
                              cTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGAACCTCATCTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGAGATCGATCATTCGCCACCCACGTTCCAacac                                                                                                            >  4:26936/1‑149 (MQ=255)
                                     gCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTc                                                                                                      >  8:73040/1‑149 (MQ=255)
                                      ccAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCg                                                                                                     >  6:2580/1‑149 (MQ=255)
                                          cGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTgg                                                                                                 <  5:151184/149‑1 (MQ=255)
                                           gagaTGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTggg                                                                                                >  8:115360/1‑149 (MQ=255)
                                                 ccGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCt                                                                                           <  3:20057/148‑1 (MQ=255)
                                                   gATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACg                                                                                        >  4:213568/1‑149 (MQ=255)
                                                       aGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCa                                                                                     <  5:2580/148‑1 (MQ=255)
                                                        gCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGTATCAGCGAGCGCGGGCTGATCCCGCGTGCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGAGTCTCCGCCAGt                                                                                   >  4:243004/1‑149 (MQ=255)
                                                            tAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCTTGCGTTCGATCATTCGCCACACAGTTATCTACACCCTGTCGCTGGGCGTCTACGCCAGTGc                                                                                 <  8:61133/147‑1 (MQ=255)
                                                              gCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTgg                                                                             <  5:253047/149‑1 (MQ=255)
                                                                   aCCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCtt                                                                        >  5:181919/1‑149 (MQ=255)
                                                                     ccAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCt                                                                      <  3:92980/149‑1 (MQ=255)
                                                                      cAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTa                                                                     <  7:115360/149‑1 (MQ=255)
                                                                            atCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCt                                                                <  3:381291/148‑1 (MQ=255)
                                                                               cTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGt                                                            >  4:278503/1‑149 (MQ=255)
                                                                                  aTCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGtt                                                           <  1:160576/147‑1 (MQ=255)
                                                                                  aTCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTtgg                                                         >  2:115543/1‑149 (MQ=255)
                                                                                     aGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTAGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTtggcct                                                      >  8:151104/1‑149 (MQ=255)
                                                                                       cgcCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTtggcctgg                                                    <  7:284367/149‑1 (MQ=255)
                                                                                         ccGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTtggcctggcc                                                  >  8:181121/1‑149 (MQ=255)
                                                                                              agagctaTGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCa                                             <  1:348679/143‑1 (MQ=255)
                                                                                                 cTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGta                                          >  6:1922/1‑149 (MQ=255)
                                                                                                   cATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGtacg                                        >  3:307767/1‑149 (MQ=255)
                                                                                                    aTGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGtac                                         <  3:278503/147‑1 (MQ=255)
                                                                                                         gCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGtacggctacg                                  <  3:213568/149‑1 (MQ=255)
                                                                                                           tGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGtacggctacggc                                >  5:41253/1‑149 (MQ=255)
                                                                                                            gTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTc                                                                       <  6:416507/109‑1 (MQ=255)
                                                                                                            gTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCACGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCAGGGCCTTc                                                                       >  5:416507/1‑109 (MQ=255)
                                                                                                                 gCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCtt                          >  4:177623/1‑149 (MQ=255)
                                                                                                                   gTGGCCTGGATCAGCGAGCGCGGGCGGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACTAGTACGGCTACGGCTTCCTTGc                        <  7:61262/149‑1 (MQ=255)
                                                                                                                        cTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGc                        >  6:129107/1‑144 (MQ=255)
                                                                                                                        cTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGc                        <  5:129107/144‑1 (MQ=255)
                                                                                                                        cTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGt                    >  4:292655/1‑148 (MQ=255)
                                                                                                                        cTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGt                    >  2:26391/1‑148 (MQ=255)
                                                                                                                         tGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTAc                                               <  6:10707/120‑1 (MQ=255)
                                                                                                                         tGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTAc                                               >  5:10707/1‑120 (MQ=255)
                                                                                                                                tcgaacgCGGGATCAGCCCGCGCTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTgggg           >  1:262083/2‑149 (MQ=255)
                                                                                                                                tcgaacgCGGGATCAGCCCGCGCTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTgggg           >  1:308832/2‑149 (MQ=255)
                                                                                                                                gcgagcgCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTgggg           <  5:1922/149‑1 (MQ=255)
                                                                                                                                gcgagcgCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTgggg           <  1:238818/149‑1 (MQ=255)
                                                                                                                                gcgagcgCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCCATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTgggg           >  5:236653/1‑149 (MQ=255)
                                                                                                                                gcgagcgCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTAGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTgggg           <  7:151104/149‑1 (MQ=255)
                                                                                                                                 cgagcgCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTggggg          >  7:319412/1‑149 (MQ=255)
                                                                                                                                  gagcgCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTGGGGgt         <  7:73040/149‑1 (MQ=255)
                                                                                                                                      gcgGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGtacggcta                                    <  6:328537/118‑1 (MQ=255)
                                                                                                                                      gcgGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGtacggcta                                    >  8:327447/1‑118 (MQ=255)
                                                                                                                                      gcgGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGtacggcta                                    <  7:327447/118‑1 (MQ=255)
                                                                                                                                      gcgGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGtacggcta                                    >  5:328537/1‑118 (MQ=255)
                                                                                                                                      gcgGGATCAGCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCtt                          >  3:111980/1‑128 (MQ=255)
                                                                                                                                      gcgGGATCAGCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCtt                          <  4:111980/128‑1 (MQ=255)
                                                                                                                                         ggCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTGTCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTGGGGGTGTCTGgc  <  4:53805/149‑1 (MQ=255)
                                                                                                                                          gCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTggg            >  6:97154/1‑138 (MQ=255)
                                                                                                                                          gCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTggg            <  5:97154/138‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGGCTGATCCTGATCGTCGTCACCTTCATCCTGCAGAGCCAACGAGATGCCGATGAGCTTTAGCCTGACCCAGATGATCCTGATCAGCGCCGGGTACCTCATGGTGCTGTTCGGCGTGGCCTGGATCAGCGAGCGCGGGCTGATCCCGCGTTCGATCATTCGCCACCCACTTACCTACACCCTGTCGCTGGGCGTCTACGCCAGTGCCTGGGCCTTCTATGGCTCGGTTGGCCTGGCCTACCAGTACGGCTACGGCTTCCTTGCCTGTTACCTGGGGGTGTCTGGC  >  NC_002947/5333525‑5333810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: