Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,416,101 T→G 7.1% R157R (CGT→CGG) ‡ PP_1238 → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,416,1010TG 7.1% 211.2 / 7.1 71R157R (CGT→CGG) ‡PP_1238lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base G (2/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

GCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTCCCATCCAC  >  NC_002947/1415963‑1416247
                                                                                                                                          |                                                                                                                                                  
gCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCACGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGACGCTGGCCaag                                                                                                                                           <  5:5220/148‑3 (MQ=255)
      cGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGc                                                                                                                                      >  1:159214/1‑147 (MQ=255)
      cGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGc                                                                                                                                      >  4:108039/1‑147 (MQ=255)
       gAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCcag                                                                                                                                   <  2:253528/149‑1 (MQ=255)
               gCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAgcggc                                                                                                                            >  5:66902/1‑148 (MQ=255)
                        ccGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAgcggc                                                                                                                            >  5:218161/1‑139 (MQ=255)
                        ccGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAgcggc                                                                                                                            <  6:218161/139‑1 (MQ=255)
                        ccGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAgg                                                                                                                  >  8:12059/1‑149 (MQ=255)
                        ccGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAgg                                                                                                                  >  6:233346/1‑149 (MQ=255)
                          ggCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTg                                                                                                                >  5:269702/1‑149 (MQ=255)
                           gCAGTGCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTgc                                                                                                               <  8:192309/149‑1 (MQ=255)
                               ttcttcGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGt                                                                                                           <  5:233346/149‑1 (MQ=255)
                                tcttcGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGACGCTGGCCaa                                                                                                                                            >  7:133958/1‑114 (MQ=255)
                                tcttcGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGACGCTGGCCaa                                                                                                                                            <  8:133958/115‑2 (MQ=255)
                                  ttcGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGc                                                                                                                      >  8:51006/1‑135 (MQ=255)
                                  ttcGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGc                                                                                                                      <  7:51006/135‑1 (MQ=255)
                                    cGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTAc                                                                                                         <  5:101062/146‑1 (MQ=255)
                                     gAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACgtgt                                                                                                     <  5:325460/149‑1 (MQ=255)
                                     gAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACgtg                                                                                                      <  3:110294/148‑1 (MQ=255)
                                       ggACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTg                                                                                                                                                 >  6:230921/1‑103 (MQ=255)
                                       ggACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTg                                                                                                                                                 <  5:230921/103‑1 (MQ=255)
                                           aaCGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCAt                                                                                               >  7:300213/1‑149 (MQ=255)
                                            aCGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAgcggc                                                                                                                            >  6:280312/1‑119 (MQ=255)
                                            aCGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATg                                                                                              <  7:12059/149‑1 (MQ=255)
                                             cGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCACACGTACGCGCacct                                                                                             <  3:108039/149‑4 (MQ=255)
                                              ggCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGa                                                                                            <  6:258994/149‑1 (MQ=255)
                                                 ttCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGcc                                                                                         >  8:18211/1‑149 (MQ=255)
                                                    cGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTg                                                                                       <  2:159214/148‑1 (MQ=255)
                                                     gCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGACGCTGGCCaa                                                                                                                                            >  3:66253/1‑93 (MQ=255)
                                                     gCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGACGCTGGCCaa                                                                                                                                            <  4:66253/94‑2 (MQ=255)
                                                       aTTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATCCCGCCTGGCAGCGGTTCGATGAACTGTCGGCCTCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGc                                                                                     <  6:303370/147‑1 (MQ=255)
                                                        ttGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTg                                                                                  <  6:66902/149‑1 (MQ=255)
                                                           cTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAg                                                                               >  6:190180/1‑149 (MQ=255)
                                                             gaggagCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAgcg                                                                             >  4:263374/1‑149 (MQ=255)
                                                              aggagCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAgcg                                                                             <  1:244473/148‑1 (MQ=255)
                                                                gagCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAgc                                                                              <  5:140752/145‑1 (MQ=255)
                                                                gagCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAgc                                                                              >  6:140752/1‑145 (MQ=255)
                                                                                 tGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTAtg                                                         >  8:185250/1‑149 (MQ=255)
                                                                                   aGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTAtgtg                                                       >  5:268975/1‑149 (MQ=255)
                                                                                      tCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGt                                                     >  8:33249/1‑148 (MQ=255)
                                                                                       cAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGTACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGt                                                     >  7:249605/1‑147 (MQ=255)
                                                                                       cAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGt                                                     >  1:248142/1‑147 (MQ=255)
                                                                                        aaTACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCa                                                   <  6:268975/148‑1 (MQ=255)
                                                                                             cGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTc                                              >  3:76915/1‑148 (MQ=255)
                                                                                             cGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTc                                              >  3:264651/1‑148 (MQ=255)
                                                                                              gCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGa                                            <  7:18211/149‑1 (MQ=255)
                                                                                               ccTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAg                                           <  7:33249/149‑1 (MQ=255)
                                                                                                         tttCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCcggcgg                                 <  2:248142/149‑1 (MQ=255)
                                                                                                         tttCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCcggcgg                                 <  5:190180/149‑1 (MQ=255)
                                                                                                           tCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCggg                                <  2:93752/148‑1 (MQ=255)
                                                                                                             gATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGgcag                             <  1:239818/149‑1 (MQ=255)
                                                                                                              aTGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGgcag                             >  1:55938/1‑148 (MQ=255)
                                                                                                              aTGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGgcag                             >  8:140510/1‑148 (MQ=255)
                                                                                                               tGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGgcggca                           <  3:263374/149‑1 (MQ=255)
                                                                                                                      tCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGcc                      >  2:26302/1‑147 (MQ=255)
                                                                                                                        ggCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGc                  >  4:95825/1‑149 (MQ=255)
                                                                                                                         gCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGt                                                     <  5:3514/113‑1 (MQ=255)
                                                                                                                         gCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGt                                                     >  6:3514/1‑113 (MQ=255)
                                                                                                                           aTCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTg               <  4:264651/149‑1 (MQ=255)
                                                                                                                               cTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGtt           >  2:262195/1‑149 (MQ=255)
                                                                                                                                ttGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTc          >  8:265713/1‑149 (MQ=255)
                                                                                                                                    ccAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTCCCAt      >  2:149538/1‑149 (MQ=255)
                                                                                                                                    ccAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTCCCAt      >  8:228498/1‑149 (MQ=255)
                                                                                                                                     cAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTc          >  7:157597/1‑144 (MQ=255)
                                                                                                                                     cAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTc          <  8:157597/144‑1 (MQ=255)
                                                                                                                                     cAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTCCCATc     <  7:228498/149‑1 (MQ=255)
                                                                                                                                     cAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTCCCATc     <  2:99542/149‑1 (MQ=255)
                                                                                                                                         gTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGaa                                                                <  5:303706/86‑1 (MQ=255)
                                                                                                                                         gTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGaa                                                                >  6:303706/1‑86 (MQ=255)
                                                                                                                                         gTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTCCCATCCAc  >  5:246640/1‑148 (MQ=255)
                                                                                                                                         gTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTCCCATCCAc  <  6:246640/148‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
GCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCTGGCAGCGTTTCGATGAACTGTCGGCATCGCTTGGCCAGCGTCTGGCCAGCACTGCCAGCAGCGGCGACAGCGAGGTGCGCGTACGTGTGCGCATGGAGCCTGGCGTGCAGCGCAACACCTCCGAAGTGTATGTGGTCAGCGTCGAGCGCCCGGCGGGCAGCACTGCCGAGCCTGCGTTCCCATCCAC  >  NC_002947/1415963‑1416247

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