Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,679,142 G→T 6.0% E256D (GAG→GAT)  bioA → adenosylmethionine‑‑8‑amino‑7‑oxononanoate aminotransferase BioA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,679,1420GT 6.0% 215.1 / 3.8 67E256D (GAG→GAT) bioAadenosylmethionine‑‑8‑amino‑7‑oxononanoate aminotransferase BioA
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base T (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

CCGAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTACCA  >  NC_002947/5679000‑5679285
                                                                                                                                              |                                                                                                                                               
ccGAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGcc                                                                                                                                           >  4:150404/1‑149 (MQ=255)
   aaCACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTg                                                                                                                                        <  4:128765/149‑1 (MQ=255)
    acacCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTgg                                                                                                                                       <  2:98651/149‑1 (MQ=255)
    acacCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTgg                                                                                                                                       >  5:267029/1‑149 (MQ=255)
     caccacGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGGCGATCTCGTCGTgga                                                                                                                                      <  2:129254/149‑2 (MQ=255)
        cacGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCtt                                                                                                                                   >  2:93122/1‑149 (MQ=255)
             cTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGGCGATCTCGTCGTggataag                                                                                                                                  <  2:57485/145‑6 (MQ=255)
             cTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGGCGATCTCGTCGTggataag                                                                                                                                  >  1:57485/1‑140 (MQ=255)
             cTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGcc                                                                                                                              >  1:217016/1‑149 (MQ=255)
                 aTCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCAc                                                                                                                          <  3:150404/149‑1 (MQ=255)
                       gcCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCAc                                                                                                                    >  3:203579/1‑149 (MQ=255)
                              tcgtcgAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTc                                                                                                             <  4:153098/149‑1 (MQ=255)
                               cgtcgAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCg                                                                                                            >  4:256782/1‑149 (MQ=255)
                                    aGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGc                                                                                                       <  7:187305/149‑1 (MQ=255)
                                      ccGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGa                                                                                                     <  5:345384/149‑1 (MQ=255)
                                       cGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGc                                                                                                       >  2:356966/1‑146 (MQ=255)
                                       cGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGc                                                                                                       <  1:356966/146‑1 (MQ=255)
                                               cAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCgg                                                                                            >  8:31430/1‑149 (MQ=255)
                                                      ccGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGc                                                                                     >  5:4171/1‑149 (MQ=255)
                                                       cGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGcc                                                                                    <  2:217016/149‑1 (MQ=255)
                                                         ggggCATGCGCATGTACCACCCGGTGTACCTCAATTTGTTGGGCGCTGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCg                                                                                  <  7:53201/147‑1 (MQ=255)
                                                         gTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCg                                                                                  <  3:256782/149‑1 (MQ=255)
                                                              aTGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATAGCCGTGGGCTTCGGCCGCACCGGCACGATGTTAGCCTGCGAGCAAGCAGGCATCCTCCAGGAgtt                                                                             >  7:215629/1‑149 (MQ=255)
                                                                   cATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCTCCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCtgt                                                                        <  7:137017/149‑1 (MQ=255)
                                                                   cATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCTCCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCtgt                                                                        <  5:56277/149‑1 (MQ=255)
                                                                      gTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGAAGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGAGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGcc                                                                     >  5:226016/1‑149 (MQ=255)
                                                                        accaccCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGc                                                                                                       >  3:323384/1‑113 (MQ=255)
                                                                        accaccCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGc                                                                                                       <  4:323384/113‑1 (MQ=255)
                                                                           accCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTc                                                                 >  8:305338/1‑148 (MQ=255)
                                                                                      cTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTc                                                                            <  6:269504/126‑1 (MQ=255)
                                                                                      cTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTc                                                                            >  5:269504/1‑126 (MQ=255)
                                                                                      cTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGAc                                                     >  5:39109/1‑149 (MQ=255)
                                                                                       tCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGAcc                                                    >  5:124321/1‑149 (MQ=255)
                                                                                        cAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACcg                                                   >  2:238965/1‑149 (MQ=255)
                                                                                               ttgCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGAt                                                                               <  3:67050/114‑1 (MQ=255)
                                                                                               ttgCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGAt                                                                               >  4:67050/1‑114 (MQ=255)
                                                                                               ttgCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCt                                                                          >  4:168863/1‑119 (MQ=255)
                                                                                                   gcgAGGCCTGCGACCGCTATTACGTGCCCCTTATCCACGACGAGATCGCCGTGTGCTTCGGCAGCAACGGCACGATGTTCGCCTGCTAGCAAGCCTGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCACTGACCGGCGGCTACCTg                                        <  8:215629/149‑1 (MQ=255)
                                                                                                   gcgAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGcc                                                          >  7:236483/1‑131 (MQ=255)
                                                                                                   gcgAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGcc                                                          <  8:236483/131‑1 (MQ=255)
                                                                                                   gcgAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTg                                        >  1:300118/1‑149 (MQ=255)
                                                                                                    cgAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGc                                       <  7:305338/149‑1 (MQ=255)
                                                                                                       ggCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGc                                    >  8:164122/1‑149 (MQ=255)
                                                                                                         ccTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTg                                  <  3:57111/149‑1 (MQ=255)
                                                                                                             cGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGgccg                              >  3:373366/1‑149 (MQ=255)
                                                                                                              gACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGgccgc                             <  6:39109/149‑1 (MQ=255)
                                                                                                                ccGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCaa                                                              <  2:319918/114‑1 (MQ=255)
                                                                                                                   cTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCgcctgcc                        >  8:353341/1‑149 (MQ=255)
                                                                                                                   cTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCgcctgcc                        >  5:209886/1‑149 (MQ=255)
                                                                                                                    tATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCgcctgcct                       >  8:21417/1‑149 (MQ=255)
                                                                                                                      tGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCgcctgcct                       <  6:4171/147‑1 (MQ=255)
                                                                                                                      tGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCgcctgcct                       <  6:209886/147‑1 (MQ=255)
                                                                                                                      tGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGa                     >  2:6164/1‑149 (MQ=255)
                                                                                                                      tGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGa                     <  7:164122/149‑1 (MQ=255)
                                                                                                                            gCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGacacctt               >  4:152065/1‑144 (MQ=255)
                                                                                                                                gaagcccACGGCGATCTCGTCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGAc             >  7:307525/6‑147 (MQ=255)
                                                                                                                                cTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACaa           <  7:21417/149‑1 (MQ=255)
                                                                                                                                   aTCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAg          <  1:93122/147‑1 (MQ=255)
                                                                                                                                     ccACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGTACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGgtg       >  5:341559/1‑148 (MQ=255)
                                                                                                                                     ccACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGgtg       >  7:322319/1‑148 (MQ=255)
                                                                                                                                      cACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTa     >  3:303227/1‑149 (MQ=255)
                                                                                                                                      cACGACGAGATCGCCGGGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTa     <  6:267029/149‑1 (MQ=255)
                                                                                                                                       ccgacgaGCTCGCCGTGCGCTTCGCCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGCATACCAGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCAGACCCTGGCCGCCTCCCTGACCCCCGACAAGGTGTAc    >  6:50098/2‑149 (MQ=255)
                                                                                                                                       acgacgAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTAc    <  4:303227/149‑1 (MQ=255)
                                                                                                                                       acgacgAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCAGTGCCTGTCCAAGGCCCTGCCCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTAc    >  6:55693/1‑149 (MQ=255)
                                                                                                                                         gacgaGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGccc                                                         <  6:263524/94‑1 (MQ=255)
                                                                                                                                         gacgaGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTACCa  <  4:373366/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CCGAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTACCA  >  NC_002947/5679000‑5679285

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