Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,015,788 C→G 10.9% intergenic (‑51/+183) PP_5679 ← / ← PP_4421 hypothetical protein/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,015,7880CG10.9% 170.1 / 12.5 64intergenic (‑51/+183)PP_5679/PP_4421hypothetical protein/aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base G (2/5);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

CCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACG  >  NC_002947/5015645‑5015929
                                                                                                                                                 |                                                                                                                                              
ccGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCgg                                                                                                                                           >  3:127117/1‑149 (MQ=255)
      gCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCgg                                                                                                                                      <  2:22481/148‑1 (MQ=255)
       cAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCggg                                                                                                                                     <  4:19998/148‑1 (MQ=255)
        agagTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCt                                                                                                                                   >  2:269135/1‑149 (MQ=255)
         gagTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTCGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCtt                                                                                                                                  >  1:61035/1‑149 (MQ=255)
           gTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTGTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCtttt                                                                                                                                <  5:300157/149‑1 (MQ=255)
              aCTTCTGGGCACGATGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACt                                                                                                                             <  5:374882/149‑1 (MQ=255)
                tgctgGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACttt                                                                                                                           >  5:147923/1‑149 (MQ=255)
                    gggCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGATCCAGCTCGGTCATCCTATCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGcgcg                                                                                                                                       <  5:267344/133‑1 (MQ=255)
                    gggCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGcgcg                                                                                                                                       >  6:267344/1‑133 (MQ=255)
                     ggCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCg                                                                                                                        <  3:161906/147‑1 (MQ=255)
                       cACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGcc                                                                                                                    >  3:271272/1‑149 (MQ=255)
                          gAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTg                                                                                                                 >  1:261471/1‑149 (MQ=255)
                             gCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGTCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGa                                                                                                               <  1:269135/149‑1 (MQ=255)
                               ttCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCggg                                                                                                                                     <  4:306484/124‑1 (MQ=255)
                               ttCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCggg                                                                                                                                     >  3:306484/1‑124 (MQ=255)
                               ttCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAg                                                                                                            >  4:32736/1‑149 (MQ=255)
                                 ctcGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCATTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAg                                                                                                            <  6:68483/147‑1 (MQ=255)
                                    ggCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTAt                                                                                                         <  3:32736/147‑1 (MQ=255)
                                    ggCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATc                                                                                                        <  5:159851/148‑1 (MQ=255)
                                          ccTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCGGGCtta                                                                                                                                 >  4:314509/1‑116 (MQ=39)
                                          ccTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCGGGCtta                                                                                                                                 <  3:314509/117‑2 (MQ=39)
                                            tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATc                                                                                                        <  2:365663/140‑1 (MQ=255)
                                            tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATc                                                                                                        >  1:365663/1‑140 (MQ=255)
                                                cGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGt                                                                                            <  1:109245/149‑1 (MQ=255)
                                                 gTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTAc                                                                                          <  8:224983/149‑1 (MQ=255)
                                                      ggTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGAccc                                                                                                  >  2:366209/1‑136 (MQ=255)
                                                      ggTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTcc                                                                                      <  3:30501/148‑1 (MQ=255)
                                                         aGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtgt                                                                                  >  7:332961/1‑149 (MQ=255)
                                                           ccACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtgtg                                                                                 >  3:304943/1‑148 (MQ=255)
                                                            cACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGaa                                                                               <  5:99669/149‑1 (MQ=255)
                                                             aCTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAAc                                                                              >  3:203456/1‑149 (MQ=255)
                                                              cTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACt                                                                             >  3:361909/1‑149 (MQ=255)
                                                               tGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtgtg                                                                                 <  1:169832/144‑1 (MQ=255)
                                                               tGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtgtg                                                                                 >  2:169832/1‑144 (MQ=255)
                                                                 cttcttTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtga                                                                           <  2:61035/148‑1 (MQ=255)
                                                                    cttTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgatg                                                                         <  4:127117/147‑1 (MQ=255)
                                                                     tttGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgatga                                                                        >  6:16515/1‑147 (MQ=255)
                                                                       tGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACa                                                                      <  5:3406/147‑1 (MQ=255)
                                                                        gATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgat                                                                          >  7:160848/1‑142 (MQ=255)
                                                                        gATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgat                                                                          <  8:160848/142‑1 (MQ=255)
                                                                            ggAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCtt                                                               <  8:332961/149‑1 (MQ=255)
                                                                                  aGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCttt                                                         >  7:4081/1‑149 (MQ=255)
                                                                                  aGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCttt                                                         >  7:189527/1‑149 (MQ=255)
                                                                                      cGGGCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGAtt                                                     <  4:21899/149‑1 (MQ=255)
                                                                                        gTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATtgt                                                   <  4:203456/149‑1 (MQ=255)
                                                                                               tAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTg                                                        <  2:7344/137‑1 (MQ=255)
                                                                                               tAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTg                                                        >  1:7344/1‑137 (MQ=255)
                                                                                                   gcTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCt                                                           <  2:102426/130‑1 (MQ=255)
                                                                                                   gcTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCt                                                           >  1:102426/1‑130 (MQ=255)
                                                                                                         gTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTgg                                                             <  8:145434/122‑1 (MQ=255)
                                                                                                         gTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTgg                                                             >  7:145434/1‑122 (MQ=255)
                                                                                                          tttGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaa                                 >  6:151089/1‑148 (MQ=255)
                                                                                                                tCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGc                          >  4:276359/1‑149 (MQ=255)
                                                                                                                  cACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaaaaaagccc                         <  4:361909/148‑5 (MQ=255)
                                                                                                                    agatgagatgaATAAGCCCGCACCGGTCTACGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATATCACTACAAAAAAAAAAGCCCGcc                     <  4:304943/148‑1 (MQ=255)
                                                                                                                          gatgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaa                                 >  3:191046/1‑132 (MQ=255)
                                                                                                                          gatgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaa                                 <  4:191046/132‑1 (MQ=255)
                                                                                                                                ggcgacccGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTGGGGAACTGATGACAGCGTCTTGGCTTTGATTGTTATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCggg           <  2:4409/144‑1 (MQ=255)
                                                                                                                                 aaGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGTACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCt         <  5:16515/149‑1 (MQ=255)
                                                                                                                                   gCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAACCCGACTAAAGGCGGGCtt        >  4:337685/1‑148 (MQ=255)
                                                                                                                                       cGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTaaaa    <  5:151089/149‑1 (MQ=255)
                                                                                                                                        gCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAATTACGTCCTTGTTAACTGATGACAGCGTCTTGGCTTTGATTGTGATACCAATACAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTaaaa    <  4:271272/149‑1 (MQ=255)
                                                                                                                                             cGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGGTTAAAACg  >  3:217914/1‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACG  >  NC_002947/5015645‑5015929

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: