Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,018,065 A→T 7.3% Q167L (CAG→CTG)  PP_3545 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,018,0650AT 7.3% 228.3 / 5.5 78Q167L (CAG→CTG) PP_3545two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base A (36/36);  new base T (4/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GACGCACCACGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAGCTGCTTGAA  >  NC_002947/4017920‑4018213
                                                                                                                                                 |                                                                                                                                                    
gacgCACCACGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGgc                                                                                                                                                   <  3:241195/149‑1 (MQ=255)
gacgCACCACGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCCCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGgc                                                                                                                                                   >  7:226013/1‑149 (MQ=255)
   gCACCACGCCACTACGATATTCGCTACCACCCGTGGCTCGATCCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCGGATGCCATCACGCAGCGCGTGACCGAACAGGAACGCCTGCCCCCGGCGGAAGAAGCACTGCGCCAGgcgcc                                                                                                                                                >  2:335861/1‑148 (MQ=255)
    caccacGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCTCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCTGGCGCag                                                                                                                                               <  7:124145/149‑1 (MQ=255)
        acGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCATCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCagaaga                                                                                                                                           <  1:367026/149‑1 (MQ=255)
        acGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTCTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCag                                                                                                                                               >  8:294793/1‑145 (MQ=255)
          gCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATg                                                                                                                                         <  8:37797/149‑1 (MQ=255)
              cTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAgg                                                                                                                                     <  5:104230/149‑1 (MQ=255)
               tACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGc                                                                                                                                    >  2:204410/1‑149 (MQ=255)
                aCGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCg                                                                                                                                   >  7:23459/1‑149 (MQ=255)
                  gAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCCGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGt                                                                                                                                 >  3:214727/1‑149 (MQ=255)
                   aaaTTCGCTACCACCCGTGGCTCGATGCTCAGAAGAGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACAGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCAAGTCTCAGAAGATTGAGGCGGTc                                                                                                                                >  4:35886/1‑149 (MQ=255)
                          cTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGt                                                                                                                         >  2:70091/1‑149 (MQ=255)
                          cTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGt                                                                                                                         >  3:153181/1‑149 (MQ=255)
                           tACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGtt                                                                                                                        >  6:3233/1‑149 (MQ=255)
                            accaccCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGtt                                                                                                                        >  7:25964/1‑148 (MQ=255)
                             ccaccCGTGGCTCGATGCCCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTg                                                                                                                       <  8:226013/148‑1 (MQ=255)
                                 ccGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACcgg                                                                                                                  >  4:218855/1‑149 (MQ=255)
                                  cGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACcggc                                                                                                                 <  6:160808/149‑1 (MQ=255)
                                   gTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACcggcg                                                                                                                >  5:222340/1‑149 (MQ=255)
                                    tGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACcggcgg                                                                                                               <  1:204410/149‑1 (MQ=255)
                                        tCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGAtt                                                                                                           >  8:265633/1‑149 (MQ=255)
                                        tCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGAtt                                                                                                           >  5:335075/1‑149 (MQ=255)
                                          gATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGc                                                                                                         <  8:25964/149‑1 (MQ=255)
                                           aTGCTCAGCAGCGTATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGcc                                                                                                        >  7:249255/1‑149 (MQ=255)
                                               tCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCAt                                                                                                     >  7:229449/1‑148 (MQ=255)
                                                   cagcGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGAttt                                                                                                >  3:117256/1‑149 (MQ=255)
                                                    agcGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTc                                                                                               <  3:256012/149‑1 (MQ=255)
                                                    agcGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTc                                                                                               <  3:91307/149‑1 (MQ=255)
                                                    agcGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTc                                                                                               <  2:184732/149‑1 (MQ=255)
                                                    agcGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTc                                                                                               <  2:99748/149‑1 (MQ=255)
                                                       gAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGAttt                                                                                                >  6:94302/1‑145 (MQ=255)
                                                       gAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGAttt                                                                                                <  5:94302/145‑1 (MQ=255)
                                                       gAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTcaac                                                                                            <  6:335075/149‑1 (MQ=255)
                                                                gcggcTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACctgctg                                                                                   <  3:218855/149‑1 (MQ=255)
                                                                       ccTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCa                                                                            >  5:186374/1‑149 (MQ=255)
                                                                         tGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCACCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGAtt                                                                                                           >  8:140058/1‑116 (MQ=255)
                                                                         tGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCACCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGAtt                                                                                                           <  7:140058/116‑1 (MQ=255)
                                                                          gTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGCCCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATcc                                                                         >  8:255860/1‑149 (MQ=255)
                                                                          gTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATcc                                                                         >  8:190984/1‑149 (MQ=255)
                                                                              gCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTgg                                                                     >  4:142378/1‑149 (MQ=255)
                                                                              gCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTg                                                                      <  6:222340/148‑1 (MQ=255)
                                                                              gCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTg                                                                      <  5:374192/148‑1 (MQ=255)
                                                                                cTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACctgc                                                                                     <  4:150975/131‑1 (MQ=255)
                                                                                cTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACctgc                                                                                     >  3:150975/1‑131 (MQ=255)
                                                                                       aTCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCa                                                                               <  2:355767/130‑1 (MQ=255)
                                                                                       aTCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCa                                                                               >  1:355767/1‑130 (MQ=255)
                                                                                       aTCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTgg                                                             >  6:178551/1‑148 (MQ=255)
                                                                                               gcgcGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCg                                                    <  7:265633/149‑1 (MQ=255)
                                                                                                 tcataggcGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGGCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCgg                                                   <  5:33292/141‑1 (MQ=255)
                                                                                                     gACCGAACAGGAACGCCTGCAGCAGGCCGAAGAACCACTTCGCCAGGCGCAGAAGATGGAGGAGGTAGGCCAGTTGACCTGCGTGATTGCCCATGATTTCAACAACCTTCTGGCCAGCATCCTTGGCGCGCTTGAACTGGCCGGGCAAc                                              <  3:35886/149‑1 (MQ=255)
                                                                                                          aaCAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACAGGCCGGATTCCCCATGATTTCAACAACCTGCTCGCCAGCATCCTTCGCGCGCTGGAACTGGCCGGACAACGTCTg                                         >  8:298264/1‑149 (MQ=255)
                                                                                                           aCAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTgg                                        <  7:255860/149‑1 (MQ=255)
                                                                                                            cAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGt                                       <  4:153181/149‑1 (MQ=255)
                                                                                                            cAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGt                                       <  1:70091/149‑1 (MQ=255)
                                                                                                              ggAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTgc                                     <  8:23459/149‑1 (MQ=255)
                                                                                                               gAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTg                                                              <  5:205971/123‑1 (MQ=255)
                                                                                                               gAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTg                                                              >  6:205971/1‑123 (MQ=255)
                                                                                                                       gcagcaGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCAt                                                                                                     <  6:27289/76‑1 (MQ=255)
                                                                                                                       gcagcaGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGcc                            >  7:182435/1‑149 (MQ=255)
                                                                                                                        cagcagGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCTATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCt                           >  3:188823/1‑149 (MQ=255)
                                                                                                                              gCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGa                     >  6:22290/1‑149 (MQ=255)
                                                                                                                                cGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGaca                   >  2:256741/1‑149 (MQ=255)
                                                                                                                                 gaagaaGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGacac                  <  3:142378/149‑1 (MQ=255)
                                                                                                                                 gaagaaGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGacac                  <  1:256741/149‑1 (MQ=255)
                                                                                                                                 gaagaaGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGacac                  >  5:275604/1‑149 (MQ=255)
                                                                                                                                    gaagCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGAcaccac               <  3:34836/149‑1 (MQ=255)
                                                                                                                                    gaagCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGAcacca                >  4:19187/1‑148 (MQ=255)
                                                                                                                                     aagCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCACCAACCTGCTGGCCAGCATCCTAGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACcaca              <  1:335861/149‑1 (MQ=255)
                                                                                                                                       cttcTGCGCCTGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACcaca              >  8:204205/4‑147 (MQ=255)
                                                                                                                                        ttcTGCGCCTGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGa                     >  5:291100/3‑139 (MQ=255)
                                                                                                                                        ttcTGCGCCTGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGa                     <  6:291100/137‑1 (MQ=255)
                                                                                                                                        ttcTGCGCCTGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAg           >  3:131594/3‑149 (MQ=255)
                                                                                                                                         tcTGCGCCTGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCAGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAgc          >  7:264009/2‑149 (MQ=255)
                                                                                                                                             cgcCAGGCGCAGAAGATGGCGGCGGTCGGCCAGTGGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAgctgct      <  4:117256/149‑1 (MQ=255)
                                                                                                                                             cgcCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAgctgc       <  5:22290/148‑1 (MQ=255)
                                                                                                                                             cgcCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATGGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAgctg        <  5:178551/147‑1 (MQ=255)
                                                                                                                                                cAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAGCTGCTTg    >  1:235887/1‑148 (MQ=255)
                                                                                                                                                 aGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGACGATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAGCTGCTTGaa  >  6:110275/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GACGCACCACGCCACTACGAAATTCGCTACCACCCGTGGCTCGATGCTCAGCAGCGAATTCAGGGCGGCTACCTGTTCGCCTATGACATCACGCAGCGCGTGACCGAACAGGAACGCCTGCAGCAGGCCGAAGAAGCACTGCGCCAGGCGCAGAAGATGGAGGCGGTCGGCCAGTTGACCGGCGGGATTGCCCATGATTTCAACAACCTGCTGGCCAGCATCCTTGGCGCGCTGGAACTGGCCGGGCAACGTCTGGTGCAGCAACGCCTGGCGGACACCACACAGCTGCTTGAA  >  NC_002947/4017920‑4018213

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