Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,539,013 T→G 8.0% K658T (AAG→ACG)  PP_4872 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,539,0130TG 8.0% 191.4 / 6.3 62K658T (AAG→ACG) PP_4872hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base G (2/3);  total (27/35)
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.89e-01

CGCGCTCTGATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCG  >  NC_002947/5538868‑5539161
                                                                                                                                                 |                                                                                                                                                    
cgcgcTCTGATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCGGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTtgt                                                                                                                                                   <  2:249431/149‑1 (MQ=255)
cgcgcTCTGATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTtgt                                                                                                                                                   <  5:2529/149‑1 (MQ=255)
cgcgcTCTGATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGACCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTtgt                                                                                                                                                   <  8:140189/149‑1 (MQ=255)
  cgcTCTGATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGc                                                                                                                                                 <  6:72632/149‑1 (MQ=255)
    ctctGATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAg                                                                                                                                               <  2:101626/149‑1 (MQ=255)
     tctGATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGa                                                                                                                                              >  7:31044/1‑149 (MQ=255)
          tAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCg                                                                                                                                         >  2:237255/1‑149 (MQ=255)
          tAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCg                                                                                                                                         <  1:361653/149‑1 (MQ=255)
             gTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTg                                                                                                                                      >  4:115622/1‑149 (MQ=255)
             gTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCGTCTGCACAAGGCGttc                                                                                                                                      <  7:339348/149‑2 (MQ=255)
             gTATCGGGCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCAGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTg                                                                                                                                      <  2:356874/149‑1 (MQ=255)
                              gTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGAc                                                                                                                       <  5:39426/147‑1 (MQ=255)
                                 cAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGAtct                                                                                                                  >  7:28813/1‑149 (MQ=255)
                                  aGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCGTCTGCACAAGGCGttctt                                                                                                                                    >  6:239928/1‑127 (MQ=255)
                                  aGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCGTCTGCACAAGGCGttctt                                                                                                                                    <  5:239928/130‑4 (MQ=255)
                                       gATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCAccc                                                                                                            <  8:31044/149‑1 (MQ=255)
                                         tCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCt                                                                                                          >  1:17144/1‑149 (MQ=255)
                                            gTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCg                                                                                                       >  6:385156/1‑149 (MQ=255)
                                            gTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCg                                                                                                       >  5:419558/1‑149 (MQ=255)
                                                gTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGAtc                                                                                                                   <  2:157215/133‑1 (MQ=255)
                                                gTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctct                                                                                                   >  6:210311/1‑149 (MQ=255)
                                                gTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctc                                                                                                    <  5:348259/148‑1 (MQ=255)
                                                gTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctc                                                                                                    <  6:419558/148‑1 (MQ=255)
                                                gTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctc                                                                                                    >  6:348259/1‑148 (MQ=255)
                                                 tGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCtt                                                                                                  <  2:155663/149‑1 (MQ=255)
                                                 tGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCtt                                                                                                  <  5:210311/149‑1 (MQ=255)
                                                      gcccggcccgATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAg                                                                                              <  3:300033/148‑1 (MQ=255)
                                                        ccggcccgATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGt                                                                                            >  6:274620/1‑148 (MQ=255)
                                                                     gTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTa                                                                                               >  3:252704/1‑132 (MQ=255)
                                                                     gTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTa                                                                                               <  4:252704/132‑1 (MQ=255)
                                                                          ctcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTc                                                                         <  7:174628/149‑1 (MQ=255)
                                                                                    aCGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGt                                                               >  3:270947/1‑149 (MQ=255)
                                                                                     cGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGt                                                               <  3:115622/148‑1 (MQ=255)
                                                                                      gCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAg                                                                >  3:242588/1‑146 (MQ=255)
                                                                                      gCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAg                                                                <  4:242588/146‑1 (MQ=255)
                                                                                             cgcTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtgg                                                      >  4:162737/1‑149 (MQ=255)
                                                                                                   aCGCAGCATCGGCCAGCCCTTCTGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtgg                                                >  8:249299/1‑149 (MQ=255)
                                                                                                    cGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtggc                                               >  2:246727/1‑149 (MQ=255)
                                                                                                      cagcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtggccg                                             >  7:90539/1‑149 (MQ=255)
                                                                                                        gcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAg                                           <  8:90539/149‑1 (MQ=255)
                                                                                                         caTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAgg                                          <  3:162737/149‑1 (MQ=255)
                                                                                                         caTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAgg                                          <  5:385156/149‑1 (MQ=255)
                                                                                                         caTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAgg                                          <  1:237255/149‑1 (MQ=255)
                                                                                                              gCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTg                                                                                                 >  1:170101/1‑89 (MQ=255)
                                                                                                              gCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTg                                                                                                 <  2:170101/89‑1 (MQ=255)
                                                                                                                      ttCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACgg                             >  5:131046/1‑149 (MQ=255)
                                                                                                                          ggTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACggg                            >  8:131831/1‑146 (MQ=255)
                                                                                                                          ggTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACggg                            >  2:253573/1‑146 (MQ=255)
                                                                                                                           gTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTa                          >  4:392430/1‑147 (MQ=255)
                                                                                                                           gTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTa                          <  3:392430/147‑1 (MQ=255)
                                                                                                                              gTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctc                                                                                                    <  2:57474/70‑1 (MQ=255)
                                                                                                                              gTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctc                                                                                                    >  1:57474/1‑70 (MQ=255)
                                                                                                                                    gcgcAAGAACGCCTTGTTCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCGACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGt                <  3:150115/148‑1 (MQ=255)
                                                                                                                                     cgcAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTAc              <  5:274620/149‑1 (MQ=255)
                                                                                                                                        agcaaCGCCGTCTGCACAAGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtggcct                                             <  7:154164/112‑2 (MQ=255)
                                                                                                                                        agcaaCGCCGTCTGCACAAGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtggccg                                             >  8:154164/4‑115 (MQ=255)
                                                                                                                                               ccTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcg    <  4:270947/149‑1 (MQ=255)
                                                                                                                                                cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTAc              >  2:106865/1‑138 (MQ=255)
                                                                                                                                                cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcgc   <  2:235786/149‑1 (MQ=255)
                                                                                                                                                cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcgc   >  4:265906/1‑149 (MQ=255)
                                                                                                                                                cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcgc   <  1:246727/149‑1 (MQ=255)
                                                                                                                                                cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCCGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTAc              <  1:106865/138‑1 (MQ=255)
                                                                                                                                                 ttGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGcgcg  >  7:227759/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGCGCTCTGATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCG  >  NC_002947/5538868‑5539161

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