Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,698,154 A→C 7.3% K263T (AAG→ACG)  hslU → protease HslVU ATPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,698,1540AC 7.3% 211.4 / 7.8 70K263T (AAG→ACG) hslUprotease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base C (2/3);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

AAGCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACC  >  NC_002947/5698009‑5698300
                                                                                                                                                 |                                                                                                                                                  
aaGCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCTCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCg                                                                                                                                                 <  7:280624/149‑1 (MQ=255)
      aaGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCaa                                                                                                                                           >  8:6012/1‑149 (MQ=255)
          tGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGtt                                                                                                                                       >  3:205836/1‑149 (MQ=255)
           gAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCACGCTTGGCAACCTTg                                                                                                                                      <  3:280663/149‑1 (MQ=255)
            aaaGAAGCGTTTAAGATTGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTg                                                                                                                                      <  6:348158/148‑1 (MQ=255)
             aagaagCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGCCAACGTTGCCACGCTTGGCAACCTTgtc                                                                                                                                    <  6:227378/149‑3 (MQ=255)
              agaagCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAGGCGTGGCAACGTTggcg                                                                                                                                   >  2:134970/1‑149 (MQ=255)
                 agCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGgcggcg                                                                                                                                <  4:205836/149‑1 (MQ=255)
                      tGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAACTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGAc                                                                                                                           >  1:199087/1‑149 (MQ=255)
                             ggTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTccc                                                                                                                    >  6:349817/1‑149 (MQ=255)
                              gTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCg                                                                                                                   <  6:106984/149‑1 (MQ=255)
                                  gcgATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGaa                                                                                                               <  7:300951/149‑1 (MQ=255)
                                   cgATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAg                                                                                                              <  7:234250/149‑1 (MQ=255)
                                          gaagCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTa                                                                                                         >  1:82217/1‑147 (MQ=255)
                                                 gccgccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGAc                                                                                                <  7:4870/149‑1 (MQ=255)
                                                   cgccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACct                                                                                              <  4:372353/149‑1 (MQ=255)
                                                    gccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCCACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACctg                                                                                             <  1:127777/149‑1 (MQ=255)
                                                    gccTGGCCAACGAGCAAGAGCTCAAGGCGAAGGCCCTGGAAGAAGTCGATCAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGTCAACGTTTGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCTACctg                                                                                             <  6:338824/149‑1 (MQ=255)
                                                     ccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAACGTTGCCACGCTTGGCAACCTTgtc                                                                                                                                    <  2:195987/109‑3 (MQ=39)
                                                     ccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAACGTTGCCACGCTTGGCAACCTTgtc                                                                                                                                    >  1:195987/1‑107 (MQ=39)
                                                       tGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACgt                                                                                                                         >  7:206646/1‑118 (MQ=255)
                                                       tGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACgt                                                                                                                         <  8:206646/118‑1 (MQ=255)
                                                            aaCGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGc                                                                                      <  8:341763/148‑1 (MQ=255)
                                                               gAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGAt                                                                                  >  1:311322/1‑149 (MQ=255)
                                                                aGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATc                                                                                 >  8:88228/1‑149 (MQ=255)
                                                                   aaGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGaa                                                                              <  2:129238/149‑1 (MQ=255)
                                                                    agagCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTgg                                                                                                                                              >  2:359995/1‑84 (MQ=255)
                                                                    agagCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAg                                                                             >  8:20683/1‑149 (MQ=255)
                                                                      agCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGc                                                                                      <  6:264457/138‑1 (MQ=255)
                                                                      agCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGc                                                                                      >  5:264457/1‑138 (MQ=255)
                                                                         tCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTg                                                                         >  1:323767/1‑148 (MQ=255)
                                                                           aaggccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGa                                                                                   >  1:326440/1‑136 (MQ=255)
                                                                           aaggccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGa                                                                                   <  2:326440/136‑1 (MQ=255)
                                                                           aaggccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGCTCGaaa                                                                             <  4:389981/142‑2 (MQ=255)
                                                                           aaggccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGCTCGaaa                                                                             >  3:389981/1‑141 (MQ=255)
                                                                           aaggccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGc                                                                        <  3:110036/147‑1 (MQ=255)
                                                                             ggccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCg                                                                    >  7:379337/1‑149 (MQ=255)
                                                                                 aaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCaa                                                                >  5:301655/1‑149 (MQ=255)
                                                                                  aggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAACTCGACAAGGTTGCCCAGCGTGGCAACGGTGGCGGCGCCGACGTGTCCCGTGACGGCGTACAGCGGGCCCTGCTGCCGCGGGTCGACGGCTGCACCGGCacc                                                               >  4:293942/1‑147 (MQ=255)
                                                                                  aggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAACGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAac                                                               <  5:240499/149‑1 (MQ=255)
                                                                                    gccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAacac                                                             <  6:54442/149‑1 (MQ=255)
                                                                                      ccTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCg                                                                                       <  4:220183/121‑1 (MQ=255)
                                                                                      ccTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCg                                                                                       >  3:220183/1‑121 (MQ=255)
                                                                                      ccTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACAcc                                                            >  8:358399/1‑148 (MQ=255)
                                                                                           aaGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTg                                                      <  2:323767/149‑1 (MQ=255)
                                                                                           aaGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTg                                                      >  5:171842/1‑149 (MQ=255)
                                                                                            aGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTgg                                                     >  2:40566/1‑149 (MQ=255)
                                                                                            aGCAGTCGAACAGAAAGCTATCGTGTTCATCGACGAAATCGATAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTgg                                                     <  8:147348/149‑1 (MQ=255)
                                                                                              cAGGCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGc                                                   <  4:292189/149‑1 (MQ=255)
                                                                                                 tCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATg                                                >  6:96005/1‑149 (MQ=255)
                                                                                                      cAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCaa                                           <  1:317578/149‑1 (MQ=255)
                                                                                                      cAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCaa                                           <  4:27038/149‑1 (MQ=255)
                                                                                                          aCGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAgacc                                       <  7:20683/149‑1 (MQ=255)
                                                                                                          aCGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAgacc                                       <  6:171842/149‑1 (MQ=255)
                                                                                                          aCGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAgacc                                       <  7:88228/149‑1 (MQ=255)
                                                                                                                     ttCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCt                            >  2:309619/1‑149 (MQ=255)
                                                                                                                        aTCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGt                          >  8:149257/1‑148 (MQ=255)
                                                                                                                           gacgaAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCa                       <  2:82217/148‑1 (MQ=255)
                                                                                                                            acgaAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGACAACACCGACCACATTCTGTTCAtt                     >  5:207444/1‑149 (MQ=255)
                                                                                                                             cgaAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTg                    >  7:373886/1‑149 (MQ=255)
                                                                                                                                aaTCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCt                 >  6:88751/1‑149 (MQ=255)
                                                                                                                                    gACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGtt                         >  3:253223/1‑137 (MQ=255)
                                                                                                                                    gACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGtt                         <  4:253223/137‑1 (MQ=255)
                                                                                                                                    gACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCgg              <  2:199087/148‑1 (MQ=255)
                                                                                                                                     aCAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGt            >  5:218180/1‑149 (MQ=255)
                                                                                                                                     aCAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGt            <  6:301655/149‑1 (MQ=255)
                                                                                                                                        acgTTGCCACGCTTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCa         >  4:394903/3‑149 (MQ=255)
                                                                                                                                          gTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCAtt       >  7:125702/1‑149 (MQ=255)
                                                                                                                                               cAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCAcc  >  8:332522/1‑149 (MQ=255)
                                                                                                                                               cAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCAcc  >  4:372077/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AAGCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACC  >  NC_002947/5698009‑5698300

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: ‑