Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,463,547 C→G 5.6% D271H (GAC→CAC)  PP_3078 ← ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,463,5470CG 5.6% 227.9 / 3.8 71D271H (GAC→CAC) PP_3078ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base G (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

TTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTA  >  NC_002947/3463399‑3463693
                                                                                                                                                    |                                                                                                                                                  
ttGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTc                                                                                                                                                    >  7:179987/1‑149 (MQ=255)
   tAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCAc                                                                                                                                                  <  4:132020/148‑1 (MQ=255)
   tAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACg                                                                                                                                                 >  8:83771/1‑149 (MQ=255)
   tAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGAACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCAc                                                                                                                                                  <  1:98321/148‑1 (MQ=255)
     gCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTg                                                                                                                                               >  7:90398/1‑149 (MQ=255)
      cccAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGc                                                                                                                                              >  7:58722/1‑149 (MQ=255)
       ctaaGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGcc                                                                                                                                             >  4:185902/3‑149 (MQ=255)
       ccAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGcc                                                                                                                                             >  8:239387/1‑149 (MQ=255)
        cAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGTGGGCACGTGACCTGccc                                                                                                                                            >  2:139798/1‑149 (MQ=255)
        cAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGcc                                                                                                                                             >  2:43003/1‑148 (MQ=255)
          aGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCccac                                                                                                                                          >  8:84588/1‑149 (MQ=255)
            gtgtAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGTGGGCACGTGACCTGCccgcca                                                                                                                                        <  1:139798/149‑1 (MQ=255)
                 ggCCGTTGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCAc                                                                                                                                                  <  4:353853/134‑1 (MQ=255)
                 ggCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCAc                                                                                                                                                  >  3:353853/1‑134 (MQ=255)
                 ggCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGcc                                                                                                                                    <  3:228288/148‑1 (MQ=255)
                               aCTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGaa                                                                                                                      <  8:90398/148‑1 (MQ=255)
                               aCTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAg                                                                                                                     <  8:58722/149‑1 (MQ=255)
                                cTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGaa                                                                                                                      <  2:181349/147‑1 (MQ=255)
                                cTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGaa                                                                                                                      >  1:181349/1‑147 (MQ=255)
                                 tCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGt                                                                                                                   >  2:169137/1‑149 (MQ=255)
                                 tCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGt                                                                                                                   >  4:57863/1‑149 (MQ=255)
                                   gcgGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACTACGCCGGGTCGCGTTCGAAGGTGa                                                                                                                 <  7:239387/149‑1 (MQ=255)
                                       ttGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGAtgct                                                                                                             <  1:305326/149‑1 (MQ=255)
                                            gTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGc                                                                                                        <  4:125002/149‑1 (MQ=255)
                                             tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGcc                                                                                                       >  3:22947/1‑149 (MQ=255)
                                               gTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGt                                                                                                     >  1:68306/1‑149 (MQ=255)
                                               gTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGt                                                                                                     >  8:90909/1‑149 (MQ=255)
                                                tAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGtt                                                                                                    <  1:169137/149‑1 (MQ=255)
                                                  gCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGt                                                                                                  <  7:253237/149‑1 (MQ=255)
                                                   ccACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTc                                                                                                 >  7:290677/1‑149 (MQ=255)
                                                   ccACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTc                                                                                                 >  5:388592/1‑149 (MQ=255)
                                                    cACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGTGGGCACGTGACCTGcc                                                                                                                                             >  6:24724/1‑104 (MQ=255)
                                                    cACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGTGGGCACGTGACCTGcc                                                                                                                                             <  5:24724/104‑1 (MQ=255)
                                                    cACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCAc                                                                                                                                                  <  6:221245/99‑1 (MQ=255)
                                                    cACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCAc                                                                                                                                                  >  5:221245/1‑99 (MQ=255)
                                                      ccGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGAt                                                                                              <  3:116118/149‑1 (MQ=255)
                                                        gCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGa                                                                                               >  7:332401/1‑146 (MQ=255)
                                                        gCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGa                                                                                               <  8:332401/146‑1 (MQ=255)
                                                             tGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGcc                                                                                         <  7:90909/147‑1 (MQ=255)
                                                              gAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCaa                                                                                      >  4:138416/1‑149 (MQ=255)
                                                                    cTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGt                                                                                <  3:57863/149‑1 (MQ=255)
                                                                    cTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTAACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCgg                                                                                 >  2:99800/1‑148 (MQ=255)
                                                                       gcgcgcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATg                                                                             <  7:83771/149‑1 (MQ=255)
                                                                       gcgcgcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATg                                                                             <  1:99800/149‑1 (MQ=255)
                                                                         gcgcCAGTTCGGTGTCGCCGAAGTACTCGATGCGCTGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCg                                                                                  <  6:50564/142‑1 (MQ=255)
                                                                         gcgcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCggagatcc                                                                           >  5:50564/1‑143 (MQ=255)
                                                                                ttCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGc                                                                    >  8:274055/1‑149 (MQ=255)
                                                                                         gCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCTATGCGGTATGGCCCGCTGCCGACTGGc                                                            <  3:185902/148‑1 (MQ=255)
                                                                                         gCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGc                                                            <  7:84588/148‑1 (MQ=255)
                                                                                         gCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCt                                                           >  4:239969/1‑149 (MQ=255)
                                                                                          ccGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGcc                                                                                                       <  8:326755/104‑1 (MQ=255)
                                                                                          ccGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGcc                                                                                                       >  7:326755/1‑104 (MQ=255)
                                                                                          ccGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCtt                                                          <  4:22947/149‑1 (MQ=255)
                                                                                           cGAAGTTCTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTg                                                         >  7:283937/1‑149 (MQ=255)
                                                                                             aaGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTc                                                       <  1:43003/149‑1 (MQ=255)
                                                                                                  cTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGc                                                  >  3:129284/1‑149 (MQ=255)
                                                                                                         cgcAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCAc                                            <  8:290677/148‑1 (MQ=255)
                                                                                                          gcAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCgg                                               <  5:356730/144‑1 (MQ=255)
                                                                                                          gcAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCgg                                               >  6:356730/1‑144 (MQ=255)
                                                                                                          gcAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCAGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCAACg                                          >  7:328096/1‑149 (MQ=255)
                                                                                                                    gAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAg                                 <  4:129284/148‑1 (MQ=255)
                                                                                                                    gAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGt                                >  4:205354/1‑149 (MQ=255)
                                                                                                                            aGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCACCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCgctgc                        >  3:52642/1‑149 (MQ=255)
                                                                                                                              cGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGctgctg                      <  3:239969/149‑1 (MQ=255)
                                                                                                                               ggCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCtgctgc                     <  2:68306/149‑1 (MQ=255)
                                                                                                                                        tGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGc            >  2:100203/1‑149 (MQ=255)
                                                                                                                                          gCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTc          >  3:225697/1‑149 (MQ=255)
                                                                                                                                            gggCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTGGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCgg        <  8:179987/149‑1 (MQ=255)
                                                                                                                                              gCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCggg       >  7:100623/1‑148 (MQ=255)
                                                                                                                                                aGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCAt    <  6:388592/149‑1 (MQ=255)
                                                                                                                                                aGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTACAGGTCGTGCTCGGGCAt    >  7:349869/1‑149 (MQ=255)
                                                                                                                                                  gTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTa  >  6:367251/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTA  >  NC_002947/3463399‑3463693

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