Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,333,354 G→A 8.7% P151L (CCG→CTG)  PP_2931 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,333,3540GA 8.7% 196.7 / 6.4 69P151L (CCG→CTG) PP_2931hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base A (2/4);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

GCCGCGCTCCTGCATGACGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTAA  >  NC_002947/3333206‑3333500
                                                                                                                                                    |                                                                                                                                                  
gacgcgCTCCTGCATGACGCGCACGGCCTCGACCTTGGCCTGATTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGCCCTGGATCATCTCGCCCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCCAAGGTATGGCCCCCAGCg                                                                                                                                                    >  8:114184/3‑149 (MQ=255)
 ccGCGCTCCTGCATGACGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCCCATCAGCAGCTGCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCgg                                                                                                                                                   <  7:114184/149‑1 (MQ=255)
             aTGACGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGctgctg                                                                                                                                       <  8:268889/149‑1 (MQ=255)
               gACGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGctgctg                                                                                                                                       >  8:215266/1‑147 (MQ=255)
                aCGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCAGCTGCAGCCGCTgggg                                                                                                                                    <  1:17572/149‑3 (MQ=255)
                aCGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCAGCTGCAGCCGCTggg                                                                                                                                     <  8:69700/148‑2 (MQ=255)
                  gcgcACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCAGCTGCAGCCGCTgggg                                                                                                                                    <  3:241616/147‑3 (MQ=255)
                  gcgcACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCAGCTGCAGCCGCTgggg                                                                                                                                    >  4:241616/1‑145 (MQ=255)
                           cTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATgg                                                                                                                         >  8:15915/1‑149 (MQ=255)
                              gACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGgca                                                                                                                       <  7:215266/148‑1 (MQ=255)
                                 cTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCAt                                                                                                                           >  7:233193/1‑141 (MQ=255)
                                 cTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCAt                                                                                                                           <  8:233193/141‑1 (MQ=255)
                                 cTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCAt                                                                                                                   >  5:79232/1‑149 (MQ=255)
                                    ggCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATc                                                                                                                  <  5:199481/147‑1 (MQ=255)
                                     gCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAgctgc                                                                                                                                         <  6:282734/123‑1 (MQ=255)
                                     gCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAg                                                                                                               <  2:239282/149‑1 (MQ=255)
                                     gCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAg                                                                                                               <  2:157432/149‑1 (MQ=255)
                                     gCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAg                                                                                                               <  2:136584/149‑1 (MQ=255)
                                                 aTCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGt                                                                                                    >  1:220027/1‑148 (MQ=255)
                                                 aTCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGt                                                                                                    <  2:220027/148‑1 (MQ=255)
                                                  tCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCg                                                                                                  >  3:59303/1‑149 (MQ=255)
                                                  tCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGACGGTGacgacg                                                                                                  >  8:345983/1‑149 (MQ=255)
                                                    cGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGgt                                                                                                <  7:16336/149‑1 (MQ=255)
                                                     gCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGt                                                                                                         >  6:354918/1‑139 (MQ=255)
                                                     gCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGt                                                                                                         <  5:354918/139‑1 (MQ=255)
                                                        ccTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGtt                                                                                             >  3:120578/1‑148 (MQ=255)
                                                          tGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGAc                                                                                          >  3:342435/1‑149 (MQ=255)
                                                              cgctgcgctgCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGcac                                                                                      <  7:137916/149‑1 (MQ=255)
                                                                ctgcgctgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGc                                                                                                                                                  <  4:18800/87‑1 (MQ=255)
                                                                ctgcgctgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGc                                                                                                                                                  >  3:18800/1‑87 (MQ=255)
                                                                   cgctgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCt                                                                                 >  6:38254/1‑149 (MQ=255)
                                                                      tgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGAc                                                                                          <  1:313884/137‑1 (MQ=255)
                                                                      tgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGAc                                                                                          >  2:313884/1‑137 (MQ=255)
                                                                      tgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGc                                                                              >  7:189113/1‑149 (MQ=255)
                                                                      tgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGc                                                                              <  4:2483/149‑1 (MQ=255)
                                                                           ggcggcCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACg                                                                         >  3:222450/1‑149 (MQ=255)
                                                                              ggcCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGt                                                                      >  3:181884/1‑149 (MQ=255)
                                                                               gcCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTg                                                                     <  4:342435/149‑1 (MQ=255)
                                                                                ccTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTgg                                                                    >  5:242272/1‑149 (MQ=255)
                                                                                 cTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGgtg                                                                                               >  2:354966/1‑121 (MQ=255)
                                                                                 cTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGgtg                                                                                               <  1:354966/121‑1 (MQ=255)
                                                                                      tgatCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCgctg                                                              <  7:15915/149‑1 (MQ=255)
                                                                                         tctcGCTCATCAGCAGGTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGctgctg                                                           >  3:317320/1‑149 (MQ=255)
                                                                                         tctcGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGctgctg                                                           >  6:288754/1‑149 (MQ=255)
                                                                                           tcGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGt                                                          >  8:296435/1‑148 (MQ=255)
                                                                                             gCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCAc                                                       <  4:222450/149‑1 (MQ=255)
                                                                                                     gcagcTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGAtt                                               >  6:378011/1‑149 (MQ=255)
                                                                                                         cTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCg                                           >  5:55941/1‑149 (MQ=255)
                                                                                                          tCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGt                                          >  1:328646/1‑149 (MQ=255)
                                                                                                           ccATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTg                                         <  4:59303/149‑1 (MQ=255)
                                                                                                           ccATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTg                                         <  7:296435/149‑1 (MQ=255)
                                                                                                           ccATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGAAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGACGCTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCg                                                                 >  2:316816/1‑125 (MQ=255)
                                                                                                            cATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTg                                         <  4:120578/148‑1 (MQ=255)
                                                                                                                   gCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCg                                 >  4:358069/1‑149 (MQ=255)
                                                                                                                   gCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCg                                 >  6:333799/1‑149 (MQ=255)
                                                                                                                       gTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGc                             >  5:144889/1‑149 (MQ=255)
                                                                                                                       gTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGc                             >  3:48348/1‑149 (MQ=255)
                                                                                                                            ttGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCg                        >  7:314561/1‑149 (MQ=255)
                                                                                                                               aCCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGc                      <  6:55941/148‑1 (MQ=255)
                                                                                                                                 cGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCtt                   >  7:182588/1‑149 (MQ=255)
                                                                                                                                    aGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGa                <  4:317320/149‑1 (MQ=255)
                                                                                                                                    aGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGa                <  8:189113/149‑1 (MQ=255)
                                                                                                                                      gTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGc              >  7:140531/1‑149 (MQ=255)
                                                                                                                                            cccccAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTg        >  2:194380/1‑149 (MQ=255)
                                                                                                                                            cccccAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTg        <  4:181884/149‑1 (MQ=255)
                                                                                                                                              cccAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTggg      <  2:328646/149‑1 (MQ=255)
                                                                                                                                              cccAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTggg      <  6:41754/149‑1 (MQ=255)
                                                                                                                                              cccAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTggg      <  5:288754/149‑1 (MQ=255)
                                                                                                                                               ccAGCAGCTGCAGCCGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGc                      >  4:28061/1‑132 (MQ=255)
                                                                                                                                               ccAGCAGCTGCAGCCGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGc                      <  3:28061/132‑1 (MQ=255)
                                                                                                                                                  gCGGCTGCAGCTGCTGGAGAAAGTTTATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCGCTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTaa  >  8:103567/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCCGCGCTCCTGCATGACGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTAA  >  NC_002947/3333206‑3333500

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