Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,885,053 G→T 7.7% R173L (CGT→CTT) ‡ PP_2540 → short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,885,0530GT 7.7% 206.3 / 6.5 65R173L (CGT→CTT) ‡PP_2540short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base T (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

GTAGATGGCATGGCCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGTCGCGTTTTGG  >  NC_002947/2884908‑2885200
                                                                                                                                                 |                                                                                                                                                   
gtAGATGGCATGGCCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTc                                                                                                                                                   <  6:192366/148‑1 (MQ=255)
  aGATGGCATGGCCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTgg                                                                                                                                                >  7:55879/1‑149 (MQ=255)
   gATGGCATGGCCGAGAGAGGTATTGGCCGGGTGGTGAATATAACCTCAGGCGCGGTGAAGGCACCGATCGACAGGCTGGGATTATCCAAGGGTGCCCGCAGTGGCGGGCCCGGATTCATAGCCGCGCTGGCCCGGCACCGGCGGCTggg                                                                                                                                               >  2:290350/1‑148 (MQ=255)
       gcatggcCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGCCAGACGCTGCTGCCGGGc                                                                                                                                           <  4:101860/149‑1 (MQ=255)
         atggcCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTg                                                                                                                                         >  5:15343/1‑149 (MQ=255)
         atggcCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTg                                                                                                                                         >  8:2325/1‑149 (MQ=255)
                aGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAAtgtg                                                                                                                                  <  7:2325/149‑1 (MQ=255)
                 gCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGa                                                                                                                                 <  5:60092/149‑1 (MQ=255)
                       ttttGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTca                                                                                                                           <  8:54709/149‑1 (MQ=255)
                       ttttGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGCCAGACGCTGCTGCCGGGc                                                                                                                                           >  4:303043/1‑133 (MQ=255)
                       ttttGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGCCAGACGCTGCTGCCGGGc                                                                                                                                           <  3:303043/133‑1 (MQ=255)
                            gCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCaacaac                                                                                                                      >  5:106045/1‑149 (MQ=255)
                                 gtggtgAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGCCAGACGCTGCTGCCGGGcc                                                                                                                                          >  1:259977/1‑123 (MQ=255)
                                 gtggtgAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGCCAGACGCTGCTGCCGGGcc                                                                                                                                          <  2:259977/124‑2 (MQ=255)
                                      gAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACACCCTGTTGCCTg                                                                                                            >  4:66178/1‑149 (MQ=255)
                                        atatCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                          <  6:15343/149‑1 (MQ=255)
                                            cACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAAtgt                                                                                                                                   >  2:55657/1‑120 (MQ=255)
                                            cACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAAtgt                                                                                                                                   <  1:55657/120‑1 (MQ=255)
                                             aCCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGtt                                                                                                     <  1:78125/149‑1 (MQ=255)
                                               cTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTAACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAAtgt                                                                                                                                   >  3:171255/1‑117 (MQ=255)
                                               cTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTAACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAAtgt                                                                                                                                   <  4:171255/117‑1 (MQ=255)
                                                     cgcgGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACgg                                                                                             <  7:20858/149‑1 (MQ=255)
                                                          tGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGa                                                                                                  <  7:200201/139‑1 (MQ=255)
                                                              ggCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGc                                                                                    <  6:106045/149‑1 (MQ=255)
                                                               gCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGa                                                                                                  <  7:58308/134‑1 (MQ=255)
                                                               gCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGa                                                                                                  >  8:58308/1‑134 (MQ=255)
                                                                 aCCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCaa                                                                                                                                      <  6:8661/96‑1 (MQ=255)
                                                                 aCCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCaa                                                                                                                                      >  5:8661/1‑96 (MQ=255)
                                                                   cGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCaa                                                                                  >  8:292813/1‑146 (MQ=255)
                                                                   cGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAg                                                                               <  7:292813/149‑1 (MQ=255)
                                                                             gCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                          >  4:234078/1‑112 (MQ=255)
                                                                             gCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                          <  3:234078/112‑1 (MQ=255)
                                                                              cTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGAGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAgg                                                                     <  8:55879/148‑1 (MQ=255)
                                                                               tGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGCCACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCt                                                                   >  5:227249/1‑149 (MQ=255)
                                                                                gggTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTg                                                                  >  3:115659/1‑149 (MQ=255)
                                                                                   tttATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCgg                                                               >  7:229958/1‑149 (MQ=255)
                                                                                      aTCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                          <  6:95599/103‑1 (MQ=255)
                                                                                      aTCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                          >  5:95599/1‑103 (MQ=255)
                                                                                       tCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTc                                                                                                    >  8:230763/1‑108 (MQ=255)
                                                                                       tCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTc                                                                                                    <  7:230763/108‑1 (MQ=255)
                                                                                         cAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGggtggt                                                          <  8:229958/148‑1 (MQ=255)
                                                                                          aaCGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCt                                                                                      <  2:306474/119‑1 (MQ=255)
                                                                                          aaCGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGa                                                        >  4:207422/1‑149 (MQ=255)
                                                                                          aaCGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGa                                                        >  4:134660/1‑149 (MQ=255)
                                                                                           aCGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGAt                                                       <  3:134660/149‑1 (MQ=255)
                                                                                              gTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGAtgtg                                                    <  3:153746/149‑1 (MQ=255)
                                                                                                 cccGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTgg                                                   >  1:12034/1‑147 (MQ=255)
                                                                                                  ccGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAg                                                >  7:4165/1‑149 (MQ=255)
                                                                                                    ggaGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTGAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTc                                              <  1:290350/147‑1 (MQ=255)
                                                                                                        tGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                          >  5:35111/1‑85 (MQ=255)
                                                                                                        tGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTggg                                                                                                          <  6:35111/85‑1 (MQ=255)
                                                                                                            cTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGcc                                        >  4:270937/1‑147 (MQ=255)
                                                                                                                  ggATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAggcgcg                                >  3:168980/1‑149 (MQ=255)
                                                                                                                  ggATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAggcgcg                                >  6:40050/1‑149 (MQ=255)
                                                                                                                    aTTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAggcgcggc                              <  3:270937/149‑1 (MQ=255)
                                                                                                                         ttGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGTCGCGCCTGCAACCGACCTTGAAGGCTGCGGGCGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAg                         >  7:205928/1‑149 (MQ=255)
                                                                                                                           gCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAgg                        <  8:4165/148‑1 (MQ=255)
                                                                                                                            cccgCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACACAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCACCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGc                       >  6:27849/4‑148 (MQ=255)
                                                                                                                                  tGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCCTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCg                >  4:133416/1‑149 (MQ=255)
                                                                                                                                      ccGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGt            >  3:97539/1‑149 (MQ=255)
                                                                                                                                           agcagcGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTCGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTAGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGTCGCGt       >  8:266601/1‑149 (MQ=255)
                                                                                                                                            gcagcGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGTCGCGtt      <  3:207422/149‑1 (MQ=255)
                                                                                                                                               gcGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGTCGCGTTTTg   >  3:275029/1‑149 (MQ=255)
                                                                                                                                                cGTCTGGCGGGCTTCAATGTGACACTCATCATCCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCTCTGCCCAGGCGCGGCGCAAGGCAATTCCGGGGTCGCGTTTTgg  <  5:27849/149‑1 (MQ=255)
                                                                                                                                                cGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGTCCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGTCGCGTTTTgg  <  6:260239/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTAGATGGCATGGCCGAGCGTGGTTTTGGCCGGGTGGTGAATATCACCTCCGGCGCGGTGAAGGCACCGATCGACATGCTGGGTTTATCCAACGGTGCCCGCAGTGGCCTGACCGGATTCATTGCCGGCCTGGCCCGGCAGCAGCGTCTGGCGGGCTGCAATGTGACACTCAACAACCTGTTGCCTGGGGCGTTCGACACGGCGCGCCTGCAAAAGACCTTGAAGGCTGCGGGTGGTGATGTGGAAGTCACTGCCCAGGCGCGGCGCAAGGCAATTCCGGCGTCGCGTTTTGG  >  NC_002947/2884908‑2885200

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