Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,246,077 G→T 6.5% R862R (CGG→AGG)  PP_4623 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,246,0770GT 6.5% 193.4 / 2.9 61R862R (CGG→AGG) PP_4623hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base T (2/2);  total (28/33)
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.44e-01

CAGCACCACCCCAATGCCATCGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGTCCATCCACAACGGATAAAGGTTGCACGCCGCGTCGATCAAGGTCCCC  >  NC_002947/5245931‑5246225
                                                                                                                                                  |                                                                                                                                                    
cAGCACCACCCCAATGCCATCGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCg                                                                                                                                                      >  8:43077/1‑147 (MQ=255)
  gCACCACCCCAATGCCATCGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGc                                                                                                                                                  >  5:79612/1‑149 (MQ=255)
    accaccCGAATGCCTTCGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAAGGGCTCGCTGTCTTGGAAATAACCTGCCTAGTCGTGCTCAGGCAAAGTCCTTGCCAGCTCGCTTAGGGTCCGCAGCGt                                                                                                                                                <  2:111279/149‑1 (MQ=255)
     ccaccCCAATGCCATCGCCGGCTCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTg                                                                                                                                               <  6:360215/149‑1 (MQ=255)
     ccaccCCAATGCCATCGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGCGGACCCTCAGCgag                                                                                                                                               >  4:91945/1‑147 (MQ=255)
      caccCCAATGCCATCGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGCGGACCCTCAGCgagc                                                                                                                                              <  3:91945/149‑4 (MQ=255)
      caccCCAATGCCATCGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTgg                                                                                                                                              >  6:19726/1‑149 (MQ=255)
       accCCAATGCCATCGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTggg                                                                                                                                             <  5:391033/149‑1 (MQ=255)
                 cATCGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTc                                                                                                                                    >  5:40653/1‑148 (MQ=255)
                    cGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGCGGACCCTCAGCgag                                                                                                                                               >  6:51375/1‑132 (MQ=255)
                    cGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGCGGACCCTCAGCgag                                                                                                                                               <  5:51375/134‑3 (MQ=255)
                      ccGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTc                                                                                                                              <  4:332237/149‑1 (MQ=255)
                         gcgcgcAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTgg                                                                                                                           <  6:79612/149‑1 (MQ=255)
                          cgcgcAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGc                                                                                                                          <  7:43077/149‑1 (MQ=255)
                           gcgcAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTggg                                                                                                                                             >  6:322608/1‑129 (MQ=255)
                           gcgcAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTggg                                                                                                                                             <  5:322608/129‑1 (MQ=255)
                               aGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATa                                                                                                                     <  2:138601/149‑1 (MQ=255)
                                gTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAg                                                                                                                    <  6:245965/149‑1 (MQ=255)
                                 ttGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATTACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGa                                                                                                                   <  5:366791/149‑1 (MQ=255)
                                    ctctTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTgggg                                                                                                                                            >  1:179841/1‑121 (MQ=255)
                                    ctctTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTgggg                                                                                                                                            <  2:179841/121‑1 (MQ=255)
                                       ttCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCaa                                                                                                               <  5:19726/147‑1 (MQ=255)
                                       ttCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAg                                                                                                              >  8:186627/1‑148 (MQ=255)
                                         cGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAg                                                                                                              <  7:186627/146‑1 (MQ=255)
                                         cGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGc                                                                                                           >  4:306323/1‑149 (MQ=255)
                                            tCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTc                                                                                                         >  3:48390/1‑148 (MQ=255)
                                             cAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTc                                                                                                         >  5:160426/1‑147 (MQ=255)
                                                  cGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTcg                                                                                                  <  4:48390/149‑1 (MQ=255)
                                                   gTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCt                                                                                                                             <  8:64092/121‑1 (MQ=255)
                                                   gTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCt                                                                                                                             >  7:64092/1‑121 (MQ=255)
                                                   gTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCa                                                                                                                         >  4:44855/1‑125 (MQ=255)
                                                   gTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCa                                                                                                                         <  3:44855/125‑1 (MQ=255)
                                                        aGCGCGACTGGCTCGATAACGTCACGCGGTCTAGGAAATAACCTGCCAAGTCGTGCTCACGCCAAGACCTTGCCAGCTCGCTGAGGGTCCGCAGCCTCGGGCATTGCTCGATTTCTGGCAGATAGAACAAGGTCTCAAACTCGCGGTCa                                                                                            >  4:115174/1‑149 (MQ=255)
                                                         gcgcGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCAc                                                                                           <  6:40653/149‑1 (MQ=255)
                                                                      gATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTca                                                                               <  3:306323/148‑1 (MQ=255)
                                                                        tAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTc                                                                                                   <  8:101157/126‑1 (MQ=255)
                                                                        tAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTc                                                                                                   >  7:101157/1‑126 (MQ=255)
                                                                           cGGGACGCGGTCTTGGCAATAACCTGCCATGTCGTGCTCAGGCCAAGTCCTTGCCAGCGCGCTGAGGGTCCGCAGCGTGGGGCGTTGCTGGATTTCTGGGTGATAGACCAAGGGCTCAAGCTCGCGGTCACGTGCATGCTTCAGCAGGt                                                                         <  5:224555/149‑1 (MQ=255)
                                                                           cGGCTCGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCTGTCACGAGCATGCTTCAGCAGGt                                                                         <  3:115174/149‑1 (MQ=255)
                                                                           cGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACTAGCATGCTTCAGCAGGt                                                                         <  5:143937/149‑1 (MQ=255)
                                                                             gCACGCGGTCTTGGCAATAACCTGCCATGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCtt                                                                                 <  5:384894/139‑1 (MQ=255)
                                                                             gCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCtt                                                                                 >  6:384894/1‑139 (MQ=255)
                                                                                cgcgGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCtt                                                                                 <  8:110323/136‑1 (MQ=255)
                                                                                cgcgGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCtt                                                                                 >  7:110323/1‑136 (MQ=255)
                                                                                           cAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACggg                                                         <  6:251296/149‑1 (MQ=255)
                                                                                            aaTAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACgggg                                                        <  6:160426/149‑1 (MQ=255)
                                                                                                   tGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGt                                                 >  1:137671/1‑149 (MQ=255)
                                                                                                   tGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGt                                                 >  7:110964/1‑149 (MQ=255)
                                                                                                                 cACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGAGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGTCCATCCACAACGGa                                   >  7:332477/1‑149 (MQ=255)
                                                                                                                  aCGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGTCCATCCACAACGGa                                   >  7:312618/1‑148 (MQ=255)
                                                                                                                          tCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGTCCATCCACAAcc                                     >  6:58587/1‑137 (MQ=255)
                                                                                                                          tCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGTCCATCCACAACg                                     <  5:58587/138‑1 (MQ=255)
                                                                                                                             ttGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAg                                                                                                      <  4:272145/70‑1 (MQ=255)
                                                                                                                             ttGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAg                                                                                                      >  3:272145/1‑70 (MQ=255)
                                                                                                                             ttGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAg                                                                                        <  5:247543/84‑1 (MQ=255)
                                                                                                                             ttGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAg                                                                                        >  6:247543/1‑84 (MQ=255)
                                                                                                                                ccAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGTCCATCCACAACGGATAAAGGTTGCAcgcc                    >  1:347257/1‑149 (MQ=255)
                                                                                                                                 cAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACgggg                                                        <  1:359322/112‑1 (MQ=255)
                                                                                                                                 cAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACgggg                                                        >  2:359322/1‑112 (MQ=255)
                                                                                                                                    cTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCGAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGTCCATCCACAACGGAt                                  <  7:196430/131‑1 (MQ=255)
                                                                                                                                    cTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGTCCATCCACAACGGATAAAGGTTGCACGCCGCGt                >  3:220601/1‑149 (MQ=255)
                                                                                                                                    cTCGCTGAGGGTCCGCAGCGTGGGGCCTAGCTCGATTTCTGGCAGATAGACCAAGGGCTCGAGCTCGCGGTCACGAGCATGCTTCACCAGCTAGGCACCTTGAACGGGGCAATCGTCCATCGACAACGGaa                                  >  8:196430/1‑130 (MQ=255)
                                                                                                                                                  gCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGTCCATCCACAACGGATAAAGGTTGCACGCCGCGTCGATCAAGGTcccc  <  4:220601/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAGCACCACCCCAATGCCATCGCCGGCGCGCAGTTGCTCTTCGATCAATGCGTTCAAGCGCGACTGGCTCGATAACGGCACGCGGTCTTGGCAATAACCTGCCAAGTCGTGCTCACGCCAAGTCCTTGCCAGCTCGCTGAGGGTCCGCAGCGTGGGGCCTTGCTCGATTTCTGGCAGATAGACCAAGGGCTCAAGCTCGCGGTCACGAGCATGCTTCAGCAGGTAGGCACCTTGAACGGGGCAATCGTCCATCCACAACGGATAAAGGTTGCACGCCGCGTCGATCAAGGTCCCC  >  NC_002947/5245931‑5246225

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑