Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,237,541 T→C 6.5% M1M (TTG→CTG) † PP_4614 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,237,5410TC 6.5% 172.8 / 2.4 61M1M (TTG→CTG) †PP_4614hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base C (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

CGCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGG  >  NC_002947/5237402‑5237688
                                                                                                                                           |                                                                                                                                                    
cGCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCgg                                                                                                                                              <  4:279779/148‑1 (MQ=255)
 gCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGc                                                                                                                                            >  3:195820/1‑149 (MQ=255)
 gCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGc                                                                                                                                            >  3:233452/1‑149 (MQ=255)
       ccAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCTTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTc                                                                                                                                      <  6:281169/149‑1 (MQ=255)
       ccAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTc                                                                                                                                      >  5:4092/1‑149 (MQ=255)
        cAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCTCAAGGCCGGGCGtt                                                                                                                                         >  5:143360/1‑145 (MQ=255)
        cAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCTCAAGGCCGGGCGtt                                                                                                                                         <  6:143360/145‑1 (MQ=255)
        cAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCg                                                                                                                                     >  5:316305/1‑149 (MQ=255)
        cAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCg                                                                                                                                     >  4:150842/1‑149 (MQ=255)
             aGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTtgt                                                                                                                                >  2:135797/1‑149 (MQ=255)
              gTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTtgtg                                                                                                                               <  4:233452/149‑1 (MQ=255)
                  aCGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGa                                                                                                                           >  2:305819/1‑149 (MQ=255)
                        ccgccACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAAt                                                                                                                      <  5:127141/148‑1 (MQ=255)
                         cgccACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATgg                                                                                                                    >  7:158489/1‑149 (MQ=255)
                          gccACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGgc                                                                                                                   <  3:150842/149‑1 (MQ=255)
                           ccaccaAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGgcg                                                                                                                  >  7:289697/1‑149 (MQ=255)
                               caAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTg                                                                                                               >  4:72755/1‑148 (MQ=255)
                                aaGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGaa                                                                                                             <  5:80461/149‑1 (MQ=255)
                                    tCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCCCTGAACCaa                                                                                                         >  7:159085/1‑149 (MQ=255)
                                          aCGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGc                                                                                                    <  6:4092/148‑1 (MQ=255)
                                           cgcgCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGa                                                                                                  <  8:289697/149‑1 (MQ=255)
                                                  gAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGTAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGAGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCg                                                                                           <  1:269734/149‑1 (MQ=255)
                                                   aGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCgg                                                                                          >  2:17858/1‑149 (MQ=255)
                                                    gCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGt                                                                                         >  2:105301/1‑149 (MQ=255)
                                                             cGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACa                                                                                                       <  2:299933/126‑1 (MQ=255)
                                                             cGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACa                                                                                                       >  1:299933/1‑126 (MQ=255)
                                                             cGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTg                                                                                <  4:195820/149‑1 (MQ=255)
                                                              ggCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTg                                                                                >  1:314747/1‑148 (MQ=255)
                                                              ggCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGc                                                                               <  5:42709/149‑1 (MQ=255)
                                                                        gcagcaAGCGCGGCGCCTTGGAACTGAGTCACCCGCACGCCTGATCGGAAATGCAAAACGCCCGGCCTTGCGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCAcc                                                                     <  2:314747/149‑1 (MQ=255)
                                                                         cagcaAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAAGGTCCGGTCAACCTCTGCATCACTCACCg                                                                   <  1:182157/149‑1 (MQ=255)
                                                                         cagcaAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCg                                                                   >  2:203491/1‑149 (MQ=255)
                                                                              aGCGCGGCTCCTTTGATCTGAGTGTCGCGGACGCCTGATCGGCAATGCAAAACGCCCGGCGTTGAGCCGGGCGTTTTCGTTTGTGCCGTTCAATGGCGCTGAACCATGATGCGAACGTCCGGTCAACCTCTGCATCAGTCACGGCAGaa                                                              <  6:292228/149‑1 (MQ=255)
                                                                                 gcggcgCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGc                                                            >  4:113317/1‑148 (MQ=255)
                                                                                   ggcgCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAg                                                         <  7:68718/149‑1 (MQ=255)
                                                                                    gcgcCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGa                                                        <  6:316305/149‑1 (MQ=255)
                                                                                        cTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCAt                                                    <  8:159085/149‑1 (MQ=255)
                                                                                           ggAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGaa                                                 >  7:140164/1‑149 (MQ=255)
                                                                                               cTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGaacaac                                             <  8:140164/149‑1 (MQ=255)
                                                                                                tGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAAcc                                            >  1:254852/1‑149 (MQ=255)
                                                                                                   gTGACCCGGACGCCGGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCAt                                         <  6:4448/149‑1 (MQ=255)
                                                                                                    tGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCAtt                                        <  3:72755/149‑1 (MQ=255)
                                                                                                         cGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTcgaac                                   <  1:17858/149‑1 (MQ=255)
                                                                                                          ggACGCCTGTTCGGCAATTCAAAACGCCCGGGCTTGATGCGGGGGTTTTAGTTTTTGCCGATCAATGGCGCTTAACCAACATGCGAACCTCCGGTCAACCTCTGCATCACTCACCTGCAGAAGCAAGACCATGAACAACCCATTcgaac                                   <  1:203491/149‑1 (MQ=255)
                                                                                                                cTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAg                            >  2:55246/1‑149 (MQ=255)
                                                                                                                 tGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGt                           >  3:314885/1‑149 (MQ=255)
                                                                                                                   aTCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGtg                          <  1:227934/148‑1 (MQ=255)
                                                                                                                    tCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGtgct                        <  1:305819/149‑1 (MQ=255)
                                                                                                                      ggCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCATCATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGc                                                            <  1:68155/111‑1 (MQ=255)
                                                                                                                      ggCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGc                                                            >  2:68155/1‑111 (MQ=255)
                                                                                                                      ggCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGtgctg                       >  5:273082/1‑148 (MQ=255)
                                                                                                                         aaTGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCAtt                   >  5:109795/1‑149 (MQ=255)
                                                                                                                           acgaaaaCGCCCGGCTCAAGGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGc                                                            <  8:23028/103‑1 (MQ=38)
                                                                                                                           acgaaaaCGCCCGGCTCAAGGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGc                                                            >  7:23028/4‑106 (MQ=38)
                                                                                                                             cAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGACATTCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGcc               <  1:272395/149‑1 (MQ=255)
                                                                                                                               aaaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCAtt                   <  8:92505/143‑1 (MQ=255)
                                                                                                                               aaaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCAtt                   >  7:92505/1‑143 (MQ=255)
                                                                                                                               aaaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGcccc             <  8:158489/149‑1 (MQ=255)
                                                                                                                                aaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGccccc            <  1:55246/149‑1 (MQ=255)
                                                                                                                                   gCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGACGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGaa         >  1:61890/1‑149 (MQ=255)
                                                                                                                                          cTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCggg  >  7:184852/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CGCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGGACGCCTGATCGGCAATGCAAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGG  >  NC_002947/5237402‑5237688

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: