Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,540,155 G→C 12.7% A277A (GCC→GCG) ‡ PP_4872 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,540,1550GC12.7% 158.2 / 10.0 60A277A (GCC→GCG) ‡PP_4872hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/23);  new base C (4/3);  total (33/27)
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.99e-01

GGGCATAGAGCACCGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGGCGCACGAAGGCGGTGCCGCGCACGTCGGTGATG  >  NC_002947/5540010‑5540302
                                                                                                                                                  |                                                                                                                                                   
gggCATAGAACACCGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGgcg                                                                                                                                                  >  3:292217/1‑149 (MQ=255)
     tAGAGCACCGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGc                                                                                                                                             >  1:112927/1‑149 (MQ=255)
        agCACCGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCGCCATCGAGCAGGccaa                                                                                                                                          <  5:194853/149‑3 (MQ=255)
         gCACCGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGcctgg                                                                                                                                         <  6:203730/149‑1 (MQ=255)
         acACCGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGcctgg                                                                                                                                         <  4:292217/148‑1 (MQ=255)
           aCCGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGcctgggt                                                                                                                                       <  8:266949/149‑1 (MQ=255)
           aCCGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGcctgggt                                                                                                                                       >  6:196179/1‑149 (MQ=255)
           aCCGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCAGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGcctgggt                                                                                                                                       <  8:10195/149‑1 (MQ=255)
             cGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTATAAGTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGcctgggtg                                                                                                                                      >  1:265790/1‑149 (MQ=255)
                    ttACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCGCCATCGAGCAGGccaagg                                                                                                                                        <  1:197635/139‑5 (MQ=255)
                    ttACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCGCCATCGAGCAGGccaagg                                                                                                                                        >  2:197635/1‑135 (MQ=255)
                    ttACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCGCCATCGAGCAGGccaa                                                                                                                                          >  4:310670/1‑135 (MQ=255)
                    ttACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCGCCATCGAGCAGGccaa                                                                                                                                          <  3:310670/137‑3 (MQ=255)
                      aCGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGcctgggtgtcctgggtgt                                                                                                                            >  5:284347/1‑149 (MQ=255)
                      aCGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGcctgggtgtcctgggtgt                                                                                                                            >  2:296377/1‑149 (MQ=255)
                              tttCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCAGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCATGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTAc                                                                                                                    >  4:32027/1‑149 (MQ=255)
                                 cGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTc                                                                                                                          <  3:191691/140‑1 (MQ=255)
                                 cGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTc                                                                                                                          >  4:191691/1‑140 (MQ=255)
                                  gAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGt                                                                                                                <  8:83792/149‑1 (MQ=255)
                                         ggCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTgatgat                                                                                                         >  6:225771/1‑149 (MQ=255)
                                            aGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCtt                                                                                                      <  6:284347/149‑1 (MQ=255)
                                            aGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATAGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCtt                                                                                                      <  5:225567/149‑1 (MQ=255)
                                                     aCCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTgatgat                                                                                                         <  3:249698/137‑1 (MQ=255)
                                                     aCCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTgatgat                                                                                                         >  4:249698/1‑137 (MQ=255)
                                                      ccGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACgcgc                                                                                            >  1:26707/1‑149 (MQ=255)
                                                                ccGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGcctgggt                                                                                                                                       >  6:92892/1‑96 (MQ=255)
                                                                ccGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGcctgggt                                                                                                                                       <  5:92892/96‑1 (MQ=255)
                                                                ccGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATc                                                                                                        <  3:217691/127‑1 (MQ=255)
                                                                ccGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATc                                                                                                        >  4:217691/1‑127 (MQ=255)
                                                                ccGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCg                                                                                    >  3:4186/1‑147 (MQ=255)
                                                                ccGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCg                                                                                    <  4:4186/147‑1 (MQ=255)
                                                                         ggCCTGTCTTTCGAGAATCTGCAGGACCTTTGCTTCTTTCGGGTCTTCCTGGTCGCCCTTGTCCTGCTCTATGGCGCAGGCCTGTGTGTCCTGTGTGTTAAAGTACGGGTTTATTATCTTTTCAACGCGCCCAGGGCGGTCTATAc                                                                            <  8:131575/146‑1 (MQ=255)
                                                                         ggCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTgatgat                                                                                                         <  1:75148/117‑1 (MQ=255)
                                                                         ggCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTgatgat                                                                                                         >  2:75148/1‑117 (MQ=255)
                                                                         ggCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACCTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACga                                                                          >  7:131575/1‑147 (MQ=255)
                                                                           ccGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTg                                                                       <  5:225771/149‑1 (MQ=255)
                                                                              ggCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAg                                                                    <  7:41469/149‑1 (MQ=255)
                                                                                   tCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCATCGCGCCCCGGGCGGTCGATACGGGTGGAGTCGa                                                                <  1:296377/148‑1 (MQ=255)
                                                                                      agaATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGAGCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGAtctc                                                            >  6:328786/1‑149 (MQ=255)
                                                                                          tCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCAGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCAAGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGCCGTCGATCTCGTAc                                                        >  2:56140/1‑149 (MQ=255)
                                                                                             gCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAg                                                     >  2:24748/1‑149 (MQ=255)
                                                                                               accaccTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGGTCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGcc                                                   >  5:60833/1‑149 (MQ=255)
                                                                                                   ccTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGc                                               >  1:81854/1‑149 (MQ=255)
                                                                                                          ttcgttcgGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATAc                                                                            >  8:122234/1‑113 (MQ=255)
                                                                                                          ttcgttcgGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATAc                                                                            <  7:122234/113‑1 (MQ=255)
                                                                                                          ttcgttcgGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTAt                                        >  6:14142/1‑149 (MQ=255)
                                                                                                             cgttcgGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTAt                                        <  2:26707/147‑1 (MQ=255)
                                                                                                                tcgGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGAtc                                                              >  2:23908/1‑122 (MQ=255)
                                                                                                                tcgGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGAtc                                                              <  1:23908/122‑1 (MQ=255)
                                                                                                                 cgGGTCTTCCTGGTCGCCCTCGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCgcggc                                  >  7:134959/1‑149 (MQ=255)
                                                                                                                         ccTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGGCGCACGAAg                          >  4:123504/1‑149 (MQ=255)
                                                                                                                          cTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGGCGCACGAAg                          <  5:14142/148‑1 (MQ=255)
                                                                                                                              tCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGGCGCACGAAGGCGGt                     <  2:81854/149‑1 (MQ=255)
                                                                                                                                   ccaggCCTGCGCCATCGAGCAGGCCTGGGTGTACTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGCCGCACGAAGGCGGTGCcgc                >  3:190786/4‑149 (MQ=255)
                                                                                                                                     aggCCTGCGCCATCGAGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGTAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCTGCGCACGAAGGCGGTGCcgcg               >  4:276679/2‑148 (MQ=255)
                                                                                                                                         cTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCACGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATAGCGGCTCACGAAGGCGGTGCCGCGCACGt          >  3:311367/1‑149 (MQ=255)
                                                                                                                                          tGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGGCGCACGAAGGCGGTGCCGCGCACGTc         >  8:76707/1‑149 (MQ=255)
                                                                                                                                            cTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGGTCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGGCGCACGAAGGCGGTGCCGCGCACGTCg        <  6:60833/148‑1 (MQ=255)
                                                                                                                                            cTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGGCGCACGAAGGCGGTGCCGCGCACGTCgg       >  5:47220/1‑149 (MQ=255)
                                                                                                                                                 tGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGGCGCACGAAGGCGGTGCCGCGCACGTCGGTGATg  <  1:24748/149‑1 (MQ=255)
                                                                                                                                                  ggCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGGCGCACGAAGGCGGTGCCGCGCACGTCGGTGATg  <  8:134959/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGGCATAGAGCACCGGGTCGTTACGCACCTTTTCGAACGATGGCAGGCGGATCACCGACTTTTCCCGGGTCACGGCCGGGCTGTCGAGAATCTGCACCACCTTGGCTTCGTTCGGGTCTTCCTGGTCGCCCTTGGCCTGCTCGATGGCGCAGGCCTGGGTGTCCTGGGTGTTCACGTACGGGTTGATGATCTTGTCAACGCGCCCCGGGCGGTCGATACGGGTGGAGTCGATCTCGTACCAGCCCTGCGGGGTATCGCGGCGCACGAAGGCGGTGCCGCGCACGTCGGTGATG  >  NC_002947/5540010‑5540302

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: