Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,413,433 T→A 5.7% N155I (AAC→ATC)  PP_1235 ← AhpC/TSA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,413,4330TA 5.7% 209.4 / 2.5 70N155I (AAC→ATC) PP_1235AhpC/TSA family protein
Reads supporting (aligned to +/- strand):  ref base T (35/31);  new base A (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGC  >  NC_002947/1413286‑1413579
                                                                                                                                                     |                                                                                                                                                  
gggTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGtt                                                                                                                                                   <  5:123709/149‑1 (MQ=255)
         gggATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGAACAAGGCtt                                                                                                                                           >  2:1771/1‑148 (MQ=255)
           gATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTggg                                                                                                                                        >  1:256103/1‑149 (MQ=255)
            atatTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGggc                                                                                                                                       >  5:298192/1‑149 (MQ=255)
            atatTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGAACAAGGCTTgaaa                                                                                                                                       <  1:1771/149‑4 (MQ=255)
                     gCCATCCCGCTGGCGACGCTGATCAAGGCGCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGgccgcc                                                                                                                                >  1:24130/1‑148 (MQ=255)
                        aTCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACGGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATAc                                                                                                                           <  2:95676/149‑1 (MQ=255)
                           ccGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACgg                                                                                                                         >  2:276859/1‑148 (MQ=255)
                             gatGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGTAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCa                                                                                                                       <  4:181666/146‑1 (MQ=255)
                              cTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATc                                                                                                                     <  2:256103/149‑1 (MQ=255)
                                      gCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGcc                                                                                                             >  6:222407/1‑149 (MQ=255)
                                         gATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCgg                                                                                                          >  5:260451/1‑149 (MQ=255)
                                              aGGCCGTGCGGGAGGCGTGGCCCAAGCAGGAGCCTAGCGTTTCCCCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACt                                                                                                     <  7:240685/149‑1 (MQ=255)
                                               ggCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACtt                                                                                                    >  5:118077/1‑149 (MQ=255)
                                               ggCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCAGCCAATACGGCATCCACATGCCAGGTCACtt                                                                                                    >  4:296565/1‑149 (MQ=255)
                                                gCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACttt                                                                                                   >  4:78266/1‑149 (MQ=255)
                                                 ccGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTc                                                                                                  <  5:204276/149‑1 (MQ=255)
                                                       cTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTcacgc                                                                                              >  1:87263/1‑149 (MQ=255)
                                                       cTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGAACAAGGCTTgaa                                                                                                                                        <  5:147899/106‑3 (MQ=255)
                                                       cTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGAACAAGGCTTgaa                                                                                                                                        >  6:147899/1‑104 (MQ=255)
                                                         ggATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTcacgccac                                                                                           <  6:298192/149‑1 (MQ=255)
                                                                   cccAAGCAGGAGCCTAGCGTTTCACCAATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCgg                                                                                                          >  6:317822/1‑124 (MQ=255)
                                                                   cccAAGCAGGAGCCTAGCGTTTCACCAATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCgg                                                                                                          <  5:317822/124‑1 (MQ=255)
                                                                    ccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTgg                                                                                >  3:63402/1‑149 (MQ=255)
                                                                     cAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTc                                                                                    >  1:309337/1‑144 (MQ=255)
                                                                     cAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTc                                                                                    <  2:309337/144‑1 (MQ=255)
                                                                     cAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgc                                                                               <  5:222407/149‑1 (MQ=255)
                                                                     cAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgc                                                                               <  1:276859/149‑1 (MQ=255)
                                                                           ggAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCCAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCgcagca                                                                         >  8:252133/1‑149 (MQ=255)
                                                                            gAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCAGTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTACCGCCACGCCATTCCTGGCGCAGCAc                                                                        >  4:87062/1‑149 (MQ=255)
                                                                                 tAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGTTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTg                                                                                 <  8:96108/135‑1 (MQ=255)
                                                                                 tAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGTTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTg                                                                                 >  7:96108/1‑135 (MQ=255)
                                                                                   gCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTcc                                                                 <  3:296565/149‑1 (MQ=255)
                                                                                       atggccGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGt                                                             <  5:210192/145‑1 (MQ=255)
                                                                                          aCCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTcc                                                                 <  2:53885/142‑1 (MQ=255)
                                                                                          aCCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTcc                                                                 >  1:53885/1‑142 (MQ=255)
                                                                                           ccGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGAt                                                         <  3:78266/149‑1 (MQ=255)
                                                                                            cGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGAt                                                         >  6:143485/1‑148 (MQ=255)
                                                                                                      aaTGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgc                                              >  7:125686/1‑149 (MQ=255)
                                                                                                       aTGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgct                                             <  6:260451/149‑1 (MQ=255)
                                                                                                        tGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgctg                                            >  4:247282/1‑149 (MQ=255)
                                                                                                        tGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgctg                                            <  1:41977/149‑1 (MQ=255)
                                                                                                        tGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgctg                                            >  3:59162/1‑149 (MQ=255)
                                                                                                        tGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgctg                                            <  2:24130/149‑1 (MQ=255)
                                                                                                        aattgaatGGGGCTGCCTTGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgct                                             >  6:60773/8‑149 (MQ=255)
                                                                                                         gTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgctgc                                           >  3:193701/1‑149 (MQ=255)
                                                                                                           aaTTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCgg                                         >  2:261633/1‑149 (MQ=255)
                                                                                                            gaatGGGGCTGCCTTGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCgg                                         >  2:251000/4‑149 (MQ=255)
                                                                                                                  ggCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgct                                             <  3:157266/139‑1 (MQ=255)
                                                                                                                  ggCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACg                                   <  3:247282/149‑1 (MQ=255)
                                                                                                                  ggCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGAACCCCGTCCTTGTCGATCAGGAAGGtgct                                             >  4:157266/1‑139 (MQ=255)
                                                                                                                    cTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGcc                                 >  1:287152/1‑149 (MQ=255)
                                                                                                                       ccaGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGCGCTGCGGTCAACGccc                                <  7:252133/145‑1 (MQ=255)
                                                                                                                          ggCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGt                            >  6:125371/1‑148 (MQ=255)
                                                                                                                                 ccATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTg                    <  8:125686/149‑1 (MQ=255)
                                                                                                                                 ccATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTg                    >  2:36410/1‑149 (MQ=255)
                                                                                                                                    ttCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGcc                  >  4:118901/1‑148 (MQ=255)
                                                                                                                                       aaTTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGc                                  >  1:186773/1‑129 (MQ=255)
                                                                                                                                       aaTTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGc                                  <  2:186773/129‑1 (MQ=255)
                                                                                                                                       aaTTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATa               >  7:211362/1‑148 (MQ=255)
                                                                                                                                        aTTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACa             <  5:125371/149‑1 (MQ=255)
                                                                                                                                         tttCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAg            >  6:260507/1‑149 (MQ=255)
                                                                                                                                             aaGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGctt         <  4:193701/148‑1 (MQ=255)
                                                                                                                                              aGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCgg                                         <  6:163568/115‑1 (MQ=255)
                                                                                                                                              aGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCgg                                         >  5:163568/1‑115 (MQ=255)
                                                                                                                                                ccTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGcttctt      >  4:79327/1‑148 (MQ=255)
                                                                                                                                                ccTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCttcttc     <  3:79327/149‑1 (MQ=255)
                                                                                                                                                 cTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCttcttc     <  3:118901/148‑1 (MQ=255)
                                                                                                                                                    gTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAg   <  8:211362/147‑1 (MQ=255)
                                                                                                                                                    gTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGc  <  1:261633/148‑1 (MQ=255)
                                                                                                                                                    gTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGc  <  1:251000/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGTACTCCTTGCCATACAGCTTCTTCAGC  >  NC_002947/1413286‑1413579

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: