Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 158,054 A→C 7.6% N22T (AAC→ACC)  PP_0148 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947158,0540AC 7.6% 188.2 / 3.1 62N22T (AAC→ACC) PP_0148hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/35);  new base C (3/2);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTATGTGGCGCCTGCCAGTTGCAGCCGC  >  NC_002947/157906‑158202
                                                                                                                                                    |                                                                                                                                                     
ccGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCaa                                                                                                                                                       >  8:236292/1‑149 (MQ=255)
ccGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCaa                                                                                                                                                       <  4:308083/149‑1 (MQ=255)
ccGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCaa                                                                                                                                                       <  7:236292/149‑1 (MQ=255)
  ggTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAAcc                                                                                                                                                     <  7:82040/149‑1 (MQ=255)
   gTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCt                                                                                                                                                    <  2:85008/149‑1 (MQ=255)
     cTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCACCCTGa                                                                                                                                                  >  2:234223/1‑149 (MQ=255)
     cTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCACCCTGa                                                                                                                                                  >  2:227371/1‑149 (MQ=255)
       tgtgCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGAc                                                                                                                                                 <  2:221890/148‑1 (MQ=255)
        gtgCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACg                                                                                                                                                >  1:226693/1‑148 (MQ=255)
           ccGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGAc                                                                                                                                             <  3:65554/148‑1 (MQ=255)
                 gTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGcc                                                                                                                                       <  1:2325/148‑1 (MQ=255)
                         tATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCAt                                                                                                                             >  1:166732/1‑149 (MQ=255)
                           tCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCACCCTGACGGACATTGCCCACTGTCAtt                                                                                                                            <  1:227371/148‑1 (MQ=255)
                           tCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCACCCTGACGGACATTGCCCACTGTCAtt                                                                                                                            <  1:234223/148‑1 (MQ=255)
                             cTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGtt                                                                                                                         >  6:159290/1‑149 (MQ=255)
                                aaaGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGa                                                                                                                      <  5:159290/149‑1 (MQ=255)
                                 aaGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGAt                                                                                                                     <  4:119342/149‑1 (MQ=255)
                                 aaGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACAAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGAt                                                                                                                     <  7:300/149‑1 (MQ=255)
                                   gggCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCt                                                                                                                   >  2:233179/1‑149 (MQ=255)
                                    gggaaGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGTTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTg                                                                                                                  <  4:220638/146‑1 (MQ=255)
                                    ggCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTg                                                                                                                  <  2:166732/149‑1 (MQ=255)
                                     gCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGc                                                                                                                 <  5:405338/149‑1 (MQ=255)
                                     gCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGc                                                                                                                 <  5:230691/149‑1 (MQ=255)
                                            ttgagccagcTTTCGCTTAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCgg                                                                                                           <  5:52633/145‑1 (MQ=255)
                                            ttCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGTGGTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCggg                                                                                                          >  3:43703/1‑149 (MQ=255)
                                               agccagcTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTc                                                                                                        <  6:208084/148‑1 (MQ=255)
                                               agccagcTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCa                                                                                                       <  8:81336/149‑1 (MQ=255)
                                               agccagcTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATC‑‑‑‑GGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCgg                                                                                                   <  3:283408/149‑1 (MQ=255)
                                                gccagcTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCACCCTGACGGCCATTGCCGACTGGCATTGTGCGATCTGCAGGCGGGTCAg                                                                                                      >  4:8720/1‑149 (MQ=255)
                                                     cTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTg                                                                                                 <  7:186744/149‑1 (MQ=255)
                                                         cGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGAc                                                                                              <  6:285340/148‑1 (MQ=255)
                                                            tGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAACGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACtggt                                                                                          >  1:305909/1‑149 (MQ=255)
                                                              aaGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCCGCGGTGGACtggtg                                                                                         <  7:368230/148‑1 (MQ=255)
                                                              aaGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACtggtg                                                                                         <  2:226693/148‑1 (MQ=255)
                                                              aaGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACtggtg                                                                                         <  4:328100/148‑1 (MQ=255)
                                                                gCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCCGCGGTGGACTggtggt                                                                                       >  4:367163/1‑148 (MQ=255)
                                                                      gTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTgg                                                                                <  8:37303/149‑1 (MQ=255)
                                                                       tGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGc                                                                               <  5:379319/149‑1 (MQ=255)
                                                                                gtaggtGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGc                                                                      >  1:344369/1‑149 (MQ=255)
                                                                                 taggtGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGc                                                                      <  3:270305/148‑1 (MQ=255)
                                                                                   ggtGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGcccc                                                                   >  7:105896/1‑149 (MQ=255)
                                                                                     tGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGATGGAGGGCAACCTGACGGACATTGACCACTGTCATTTTTAGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCGCGt                                                                 >  7:279211/1‑149 (MQ=255)
                                                                                           aGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGTGTTATCGGCTGGAGGGCAACGTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTgg                                                            <  4:43703/148‑1 (MQ=255)
                                                                                                   aCCGAAGGGGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGt                                                    <  3:8720/148‑1 (MQ=255)
                                                                                                     cGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTg                                                   <  5:359008/147‑1 (MQ=255)
                                                                                                     cGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTg                                                   >  6:359008/1‑147 (MQ=255)
                                                                                                             gCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACggg                                          <  1:233179/148‑1 (MQ=255)
                                                                                                               tGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGc                                         >  2:187555/1‑147 (MQ=255)
                                                                                                                gCCATTGTGGAGCCCTGCGTTAGCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGc                                         >  7:66169/1‑147 (MQ=255)
                                                                                                                  cATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTgg                                                  >  1:19989/1‑135 (MQ=255)
                                                                                                                  cATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTgg                                                  <  2:19989/135‑1 (MQ=255)
                                                                                                                   aTTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGcgagcc                                    >  6:318349/1‑148 (MQ=255)
                                                                                                                   aTTGTGGAGCCCTGCGGTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACAGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGcgagccg                                   <  8:351630/149‑1 (MQ=255)
                                                                                                                       tgGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAgct                               >  3:51738/1‑149 (MQ=255)
                                                                                                                          aGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAgctgc                             >  5:204082/1‑148 (MQ=255)
                                                                                                                            cccTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGcgagccga                                  >  7:364236/1‑141 (MQ=255)
                                                                                                                            cccTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGcgagccga                                  <  8:364236/141‑1 (MQ=255)
                                                                                                                              cTGCGTTATCGGCTGGAGGGCAACCTGACGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTAtgtg                       >  6:242737/1‑149 (MQ=255)
                                                                                                                                gCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGAGGGTCAGCGGTGGACTGGTGGTCAACTGGCTTACACTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTATGTgg                      >  7:145989/1‑149 (MQ=255)
                                                                                                                                 cGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTATGTGgc                     <  2:305909/149‑1 (MQ=255)
                                                                                                                                         gCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTATGTGGCGCCTGCCa             <  2:344369/149‑1 (MQ=255)
                                                                                                                                               gggCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTATGTGGCGCCTGCCAGTTg         >  5:263163/1‑147 (MQ=255)
                                                                                                                                                ggCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTATGTGGCGCCTGCCAGTTGCAg      >  7:210234/1‑149 (MQ=255)
                                                                                                                                                ggCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTATGTGGCGCCTGCCAGTTGCAg      <  5:318349/149‑1 (MQ=255)
                                                                                                                                                  cAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTATGTGGCGCCTGCCAGTTGCAGcc    <  6:107452/149‑1 (MQ=255)
                                                                                                                                                    aCCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTATGTGGCGCCTGCCAGTTGCAGCCGc  <  8:105896/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGGTCTTGTGCCGGTAGTGGGCGCTATCCTAAAAGGGCAAGCTTTCAGCCAGCTTTCGCTGAAGCAAAGGTGACGAGAGGTAGGTGATGCAGGATGTAACCGAAGGTGGCTGCCATTGTGGAGCCCTGCGTTATCGGCTGGAGGGCAACCTGACGGACATTGCCCACTGTCATTGTTCGATCTGCAGGCGGGTCAGCGGTGGACTGGTGGTCACCTGGCTTACCCTGCCCCGTTCTGGGTTCCAGTGGTTAACGGGCGAGCCGAGCTGCTATGTGGCGCCTGCCAGTTGCAGCCGC  >  NC_002947/157906‑158202

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: