Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,407,899 T→G 6.7% I129L (ATC→CTC)  PP_3890 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,407,8990TG 6.7% 187.4 / 5.2 60I129L (ATC→CTC) PP_3890hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base G (2/2);  total (31/29)
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.50e-01

ACGCGGGTCACGCAGGTCGCGCACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCGATGTACAGGCCGAAGTACCCACCCCGCGCCATAGGCAAC  >  NC_002947/4407754‑4408047
                                                                                                                                                 |                                                                                                                                                    
aCGCGGGTCACGCAGGTCGCGCACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGAtgct                                                                                                                                                   <  6:421995/149‑1 (MQ=255)
    gggTCACGCAGGTCGCGCACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGccc                                                                                                                                               <  3:108369/149‑1 (MQ=255)
         aCGCAGGTCGCGCACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGc                                                                                                                                           <  3:299029/148‑1 (MQ=255)
          cGCAGGTCGCGCACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGcac                                                                                                                                         <  1:167258/149‑1 (MQ=255)
           gCAGGTCGCGCACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGAt                                                                                                                                                      <  1:164325/135‑1 (MQ=255)
           gCAGGTCGCGCACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGcaca                                                                                                                                        >  5:288805/1‑149 (MQ=255)
               gTCGCGCACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGAt                                                                                                                                    <  6:1841/149‑1 (MQ=255)
                 cgcgCACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATgg                                                                                                                                  >  3:247063/1‑149 (MQ=255)
                 cgcgCACCGACTTGTCGGTCGCCGGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATgg                                                                                                                                  <  8:166782/149‑1 (MQ=255)
                  gcgcACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGc                                                                                                                                 >  8:16226/1‑149 (MQ=255)
                  gcgcACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGc                                                                                                                                 <  6:335034/149‑1 (MQ=255)
                  gcgcACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCCGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACCGCGTCGATGCTGCCCCGGCACACGATGGc                                                                                                                                 >  6:42515/1‑149 (MQ=255)
                    gcACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCa                                                                                                                               <  8:4812/149‑1 (MQ=255)
                         gACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTaa                                                                                                                          <  2:293228/149‑1 (MQ=255)
                            ttGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATtc                                                                                                                       >  3:254896/1‑149 (MQ=255)
                                 ggtcgCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCtcgtt                                                                                                                  <  4:254896/149‑1 (MQ=255)
                                     gCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCtcgttcgt                                                                                                               >  3:98432/1‑148 (MQ=255)
                                      ccTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCa                                                                                                             <  2:43467/149‑1 (MQ=255)
                                        tGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGGAATTCTCGTTCGTCAgc                                                                                                           >  4:28037/1‑149 (MQ=255)
                                            gccgcACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATTCGCCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCAAGTACGTCAgctgct                                                                                                       >  3:208039/1‑149 (MQ=255)
                                            gccgcACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAgctgct                                                                                                       >  3:408474/1‑149 (MQ=255)
                                             ccgcACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGGGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCGGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGctgctg                                                                                                      <  8:323262/149‑1 (MQ=255)
                                                      ggCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAg                                                                                               <  4:98432/147‑1 (MQ=255)
                                                      ggCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAg                                                                                               <  7:16226/147‑1 (MQ=255)
                                                        ccTCGCTCTACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGTTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGctt                                                                                           <  4:208039/149‑3 (MQ=255)
                                                            gCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCctg                                                                                       >  6:332287/1‑149 (MQ=255)
                                                                aaCTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCctggctg                                                                                   <  4:247063/149‑1 (MQ=255)
                                                                aaCTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCctggc                                                                                     <  5:332287/147‑1 (MQ=255)
                                                                     cacttcGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGc                                                                                <  4:408474/147‑1 (MQ=255)
                                                                       ctgcgcgGGTCGCGCGCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCgaga                                                                            <  6:288805/146‑1 (MQ=255)
                                                                         tcGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGAcc                                                                          >  6:262509/1‑149 (MQ=255)
                                                                             gggTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGgcg                                                                      >  6:189795/1‑149 (MQ=255)
                                                                                    gcTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGaa                                                               >  3:359992/1‑149 (MQ=255)
                                                                                      tCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAAcg                                                             >  8:232388/1‑149 (MQ=255)
                                                                                        aTGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAAcggc                                                           >  4:37922/1‑149 (MQ=255)
                                                                                        aTGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAAcggc                                                           >  3:339699/1‑149 (MQ=255)
                                                                                            ggCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCt                                                       >  2:356546/1‑149 (MQ=255)
                                                                                              ccACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTg                                                     <  5:42515/149‑1 (MQ=255)
                                                                                               cACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTgg                                                    >  2:24736/1‑149 (MQ=255)
                                                                                                      gCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGa                                              >  4:284112/1‑148 (MQ=255)
                                                                                                              gaggaggTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCGa                                        <  3:385600/146‑1 (MQ=255)
                                                                                                              gaggaggTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCGa                                        >  4:385600/1‑146 (MQ=255)
                                                                                                              gaggaggTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCGATg                                      >  1:50729/1‑148 (MQ=255)
                                                                                                              gaggaggTAGGCACGGATGGCCTCGACGGCGTCGATGCTACCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAg                                                                                               >  5:265114/1‑91 (MQ=255)
                                                                                                              gaggaggTAGGCACGGATGGCCTCGACGGCGTCGATGCTACCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAg                                                                                               <  6:265114/91‑1 (MQ=255)
                                                                                                                    gTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAAcggcg                                                          >  8:178102/1‑122 (MQ=255)
                                                                                                                            ggATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAgctg                                                                                                         >  4:389409/1‑67 (MQ=255)
                                                                                                                            ggATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAgctg                                                                                                         <  3:389409/67‑1 (MQ=255)
                                                                                                                            ggATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAAcggc                                                           >  7:206086/1‑113 (MQ=255)
                                                                                                                            ggATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAAcggc                                                           <  8:206086/113‑1 (MQ=255)
                                                                                                                                 gcatCGACGCCGTCGAGGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCctggctg                                                                                   >  3:183978/4‑84 (MQ=255)
                                                                                                                                 gcatCGACGCCGTCGAGGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCctggctg                                                                                   >  3:177560/4‑84 (MQ=255)
                                                                                                                                 gcatCGACGCCGTCGAGGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCctggctg                                                                                   <  4:177560/81‑1 (MQ=255)
                                                                                                                                 gcatCGACGCCGTCGAGGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCctggctg                                                                                   <  4:183978/81‑1 (MQ=255)
                                                                                                                                    tCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCGATGTACAGGCCGAAGTACCCACCCcg               >  4:327361/1‑149 (MQ=255)
                                                                                                                                      gACTGCGTCGATGATGCCCCTGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCTAGACCTGCGCGTAGAACGGCGGCTTGTCCTTGAACTCGATGTACATGCCGAAGTACCCACCAcgcg             >  4:298827/1‑149 (MQ=255)
                                                                                                                                        cGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCCATGTACAGGCCGAAGTACCCCCCCCgcgc            >  7:387405/1‑148 (MQ=255)
                                                                                                                                         ggCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGc                                                   <  4:87541/108‑1 (MQ=255)
                                                                                                                                         ggCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGc                                                   >  3:87541/1‑108 (MQ=255)
                                                                                                                                           cGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTAAGCTGCTGCAGGTAGGCGTCCTGGCTTGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCGATGTACAGGCCGAAGTACCCACCCCGCGCCATa        <  3:298827/149‑1 (MQ=255)
                                                                                                                                             tCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCGATGTACAGGCCGAAGTACCCACCCCGCGCCATAgg      <  3:327361/149‑1 (MQ=255)
                                                                                                                                               gATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCGATGTACAGGCCGAAGTACCCACCCCGCGCCATAGGc     <  7:232388/148‑1 (MQ=255)
                                                                                                                                                 tgctgcCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCGATGTACAGGCCGAAGTACCCACCCCGCGCCATAGGCAac  <  3:37922/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACGCGGGTCACGCAGGTCGCGCACCGACTTGTCGGTCGCCTGACGCCGCACCTCGGCCTCGCTCAACTTCACTTCGCGGGTCGCGCTCATGCGGCCACCGTGGCGGGCAGGAGGAGGTAGGCACGGATGGCCTCGACGGCGTCGATGCTGCCCCGGCACACGATGGCCAGGTAATTCTCGTTCGTCAGCTGCTGCAGGTAGGCGTCCTGGCTGGGCGAGACCGGCGCGTCGAACGGCGGCTTGGCCTTGAACTCGATGTACAGGCCGAAGTACCCACCCCGCGCCATAGGCAAC  >  NC_002947/4407754‑4408047

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