Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 850,991 C→G 5.3% intergenic (‑574/‑231) pcs ← / → hemA phosphatidylcholine synthase/glutamyl‑tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947850,9910CG 5.3% 242.0 / 3.2 75intergenic (‑574/‑231)pcs/hemAphosphatidylcholine synthase/glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base G (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

GCTCTGGCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGC  >  NC_002947/850845‑851136
                                                                                                                                                  |                                                                                                                                                  
gCTCTGGCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATcc                                                                                                                                                   <  5:237432/148‑1 (MQ=255)
gCTCTGGCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATcc                                                                                                                                                   <  4:111867/148‑1 (MQ=255)
gCTCTGGCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATcc                                                                                                                                                   <  2:69072/148‑1 (MQ=255)
  tctGGCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCAGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATg                                                                                                                                                >  7:19570/1‑149 (MQ=255)
      gCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATGGATGCCCg                                                                                                                                            <  4:34548/149‑1 (MQ=255)
       cTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACATCGCTTGACGGGCATCCATGCCCGt                                                                                                                                           <  7:44528/149‑1 (MQ=255)
       cTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATGGATGCCCGt                                                                                                                                           <  2:168593/149‑1 (MQ=255)
             ctctTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCAt                                                                                                                                                 >  6:319930/1‑137 (MQ=255)
             ctctTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCAt                                                                                                                                                 <  5:319930/137‑1 (MQ=255)
             ctctTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTc                                                                                                                                     >  4:215633/1‑149 (MQ=255)
              tctTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTcc                                                                                                                                    >  7:147989/1‑149 (MQ=255)
                ttACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCGCCATCCTCTGGCGCAGCGCTTGACGGGCATCCATGCCCGTCACCTccca                                                                                                                                  >  7:66798/1‑148 (MQ=255)
                ttACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                                          >  4:140709/1‑141 (MQ=255)
                ttACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                                          <  3:140709/141‑1 (MQ=255)
                ttACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGAGCACAGCGCTTGACGGGCATCCATGCCCGTCacctacct                                                                                                                                  >  1:144906/1‑149 (MQ=255)
                  aCTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTcc                                                                                                                                    <  8:66798/145‑1 (MQ=255)
                    tGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTgcg                                                                                                                               >  5:236490/1‑148 (MQ=255)
                      tCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTgcgc                                                                                                                              >  5:304403/1‑147 (MQ=255)
                      tCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCaa                                                                                                                            <  6:236490/149‑1 (MQ=255)
                        gCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAg                                                                                                                           <  6:163128/148‑1 (MQ=255)
                         cTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc                                                                                                                         <  6:304403/149‑1 (MQ=255)
                             tAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATg                                                                                                                                                <  8:35529/122‑1 (MQ=255)
                              aaGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCg                                                                                                                    >  3:254793/1‑149 (MQ=255)
                                 cgcgcgGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTATGGCTTGCTCCAGGGAGGACCTCGCTGCAGGCCCCATCCTAGGGCACAGCGCATGACGGTCATCCATGCCCGTCACCTACCTGCTCAAGCACTGcctac                                                                                                                 >  8:227997/1‑145 (MQ=21)
                                              aGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGTGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAg                                                                                                    <  4:254793/148‑1 (MQ=255)
                                              aGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGTCATCCATGACCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGt                                                                                                   >  6:232275/1‑149 (MQ=255)
                                                 cgcTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcgc                                                                                                <  7:227997/149‑1 (MQ=255)
                                                 cgcTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcgc                                                                                                <  8:19570/149‑1 (MQ=255)
                                                      cAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCt                                                                                                                                  <  5:304471/111‑1 (MQ=255)
                                                      cAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCt                                                                                                                                  >  6:304471/1‑111 (MQ=255)
                                                        gATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTGGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatc                                                                                         <  2:168641/149‑1 (MQ=255)
                                                          ttGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaa                                                                                       <  2:265781/149‑1 (MQ=255)
                                                           tGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagata                                                                                         >  5:217317/1‑145 (MQ=255)
                                                               aCTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagtggt                                                                                  >  7:287253/1‑145 (MQ=255)
                                                                 ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc                                                                                                                         <  5:213690/109‑1 (MQ=255)
                                                                 ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc                                                                                                                         >  6:213690/1‑109 (MQ=255)
                                                                 ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTTCTGAAGTCGCGAAGATCAAGATCAg                                                                                 >  4:310120/1‑148 (MQ=255)
                                                                   cAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCCGGCCTCCTGAAGTCGCGAAGATCAAGATCAgg                                                                                <  1:201767/148‑1 (MQ=255)
                                                                   cAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTc                                                                                                                <  2:176869/115‑1 (MQ=255)
                                                                   cAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTc                                                                                                                >  1:176869/1‑115 (MQ=255)
                                                                      gaggCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCa                                                                            >  6:326574/1‑148 (MQ=255)
                                                                          ccGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGc                                                                       >  8:63310/1‑149 (MQ=255)
                                                                            gAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGc                                                                       >  6:167918/1‑147 (MQ=255)
                                                                             aGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAg                                                                    <  3:69358/149‑1 (MQ=255)
                                                                                    ggCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCaaaa                                                             <  2:40782/149‑1 (MQ=255)
                                                                                    ggCGTGCTCCCGGGAGGACCGCGGTGCTGGCCCCTTCCTCGGGCACAGCGCTTGAGGGGCATCCATGCCCGTCACCTCCGTGCGCATGGCGTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCaaaa                                                             <  4:101755/149‑1 (MQ=18)
                                                                                     gCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGcc                                                                                                                        >  8:84442/1‑90 (MQ=255)
                                                                                     gCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGcc                                                                                                                        <  7:84442/90‑1 (MQ=255)
                                                                                     gCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAg                                                            <  3:215633/149‑1 (MQ=255)
                                                                                           tCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGc                                                       >  7:305030/1‑148 (MQ=255)
                                                                                           tCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGc                                                       >  4:228920/1‑148 (MQ=255)
                                                                                             ccGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAg                                                                        <  2:216821/129‑1 (MQ=255)
                                                                                             ccGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAg                                                                        >  1:216821/1‑129 (MQ=255)
                                                                                             ccGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCaa                                                     <  7:63310/148‑1 (MQ=255)
                                                                                                    gaCCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGa                                              >  3:180137/1‑148 (MQ=255)
                                                                                                    gaCCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGa                                              <  4:180137/148‑1 (MQ=255)
                                                                                                       cgcgCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaaa                                           >  7:259427/1‑148 (MQ=255)
                                                                                                       cgcgCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaaa                                           >  4:52754/1‑148 (MQ=255)
                                                                                                         cgcTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAacac                                        >  4:217394/1‑149 (MQ=255)
                                                                                                         cgcTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAacac                                        <  5:167918/149‑1 (MQ=255)
                                                                                                          gcTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACAcc                                       >  4:152401/1‑149 (MQ=255)
                                                                                                            tGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGc                                     >  7:71037/1‑149 (MQ=255)
                                                                                                                 gCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAggg                                >  8:8661/1‑149 (MQ=255)
                                                                                                                  ccccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAggg                                >  6:56685/1‑148 (MQ=255)
                                                                                                                      aTCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGc                           <  5:326574/149‑1 (MQ=255)
                                                                                                                      aTCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGc                           <  8:305030/149‑1 (MQ=255)
                                                                                                                       tCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCt                          <  3:228920/149‑1 (MQ=255)
                                                                                                                                acaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAgg                 <  6:295510/149‑1 (MQ=255)
                                                                                                                                acaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAgg                 <  5:56685/149‑1 (MQ=255)
                                                                                                                                 caGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGc                <  8:147989/149‑1 (MQ=255)
                                                                                                                                  aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCg               <  8:259427/149‑1 (MQ=255)
                                                                                                                                       ttGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAgg                 <  4:201124/142‑1 (MQ=255)
                                                                                                                                       ttGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAgg                 >  3:201124/1‑142 (MQ=255)
                                                                                                                                        tGACGGGCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGcc         >  8:270503/1‑149 (MQ=255)
                                                                                                                                         gACGGGCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCg        >  1:161846/1‑149 (MQ=255)
                                                                                                                                            gggCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTg     >  1:158765/1‑149 (MQ=255)
                                                                                                                                              gCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgcg   >  3:259636/1‑149 (MQ=255)
                                                                                                                                               cATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgcgc  <  3:152401/149‑1 (MQ=255)
                                                                                                                                                  ccATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGAc                                                  >  3:92656/1‑98 (MQ=255)
                                                                                                                                                  ccATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGAc                                                  <  4:92656/98‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCTCTGGCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGC  >  NC_002947/850845‑851136

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: