Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,268,089 C→T 6.9% A187V (GCG→GTG)  azoR1 → FMN‑dependent NADH‑azoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,268,0890CT 6.9% 171.6 / 2.1 58A187V (GCG→GTG) azoR1FMN‑dependent NADH‑azoreductase
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base T (2/2);  total (26/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 = 9.12e-01

CACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGTTCAGGCCACCGGGATGAGCGGAT  >  NC_002947/3267943‑3268232
                                                                                                                                                  |                                                                                                                                               
cACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGcga                                                                                                                                               <  3:100780/149‑1 (MQ=255)
cACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGcga                                                                                                                                               <  3:95821/149‑1 (MQ=255)
        gTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGTGATCGAtgct                                                                                                                                       >  6:302208/1‑149 (MQ=255)
        gTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGTGATCGAtgct                                                                                                                                       >  6:212240/1‑149 (MQ=255)
                  cgatGGCCAGCCGACTGGAGCCGGGCATGAAGACTTCCTGAAACTGTTCATGCGCTTCATCGGGATTACTGACCTGGAAATATTCCGGGCAAACGGCCTGGCCTATGGAACGGAGCAGCGCAGCCAGGCTATCGATGCTGCACAGgcg                                                                                                                              >  8:68123/4‑148 (MQ=255)
                  cGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGgcgc                                                                                                                             <  7:295581/149‑1 (MQ=255)
                  cGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGgcgc                                                                                                                             <  5:64467/149‑1 (MQ=255)
                   gCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATTGGGATTACTGACCTGGAAATCGTCCGGGCCCATGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGcac                                                                                                                                   <  8:127094/142‑1 (MQ=255)
                   gCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATTGGGATTACTGACCTGGAAATCGTCCGGGCCCATGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGcac                                                                                                                                   >  7:127094/1‑142 (MQ=255)
                    cTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAg                                                                                                                           <  5:358432/149‑1 (MQ=255)
                       gCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGTGATCGATGCTGCACAGGCGCAGAt                                                                                                                         <  5:212240/148‑1 (MQ=255)
                       gCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGTGATCGATGCTGCACAGGCGCAGAt                                                                                                                         <  5:302208/148‑1 (MQ=255)
                       gCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCGTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGAtt                                                                                                                        <  4:70209/149‑1 (MQ=255)
                         cAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGc                                                                                                                      <  5:47111/149‑1 (MQ=255)
                             cGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGc                                                                                                                  <  7:68123/149‑1 (MQ=255)
                                           cATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGc                                                                                                    >  3:184212/1‑149 (MQ=255)
                                                 gACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGt                                                                                                           <  2:14739/136‑1 (MQ=255)
                                                 gACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGt                                                                                                           >  1:14739/1‑136 (MQ=255)
                                                    ttCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGAt                                                                                            >  5:286488/1‑148 (MQ=255)
                                                               gtTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGtt                                                                                >  6:111101/1‑149 (MQ=255)
                                                                    tGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGtt                                                                                                          <  1:193377/118‑1 (MQ=255)
                                                                    tGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGtt                                                                                                          >  2:193377/1‑118 (MQ=255)
                                                                    tGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTgg                                                                             <  1:215990/147‑1 (MQ=255)
                                                                    tGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTgg                                                                             >  2:215990/1‑147 (MQ=255)
                                                                      ggCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGc                                                                         <  5:137188/149‑1 (MQ=255)
                                                                             tCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTg                                                                  >  4:78108/1‑149 (MQ=255)
                                                                                ggATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGcc                                                               <  6:337157/149‑1 (MQ=255)
                                                                                ggATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGcc                                                               <  6:286488/149‑1 (MQ=255)
                                                                                 gATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGcgatcga                                                                                                                                           >  8:91115/1‑72 (MQ=255)
                                                                                 gATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGcgatcga                                                                                                                                           <  7:91115/72‑1 (MQ=255)
                                                                                     aCTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGt                                                           <  3:18275/148‑1 (MQ=255)
                                                                                      cggACCTGGAAAGCCGGTGGGGCCACGGCCTGGCCTATGGCGCGGAGCAGCGCAGGCTGTCGATCGATGCTGCAGAGGCTCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGAGCTTGGGAGCGTTTGTTTGCAAGCGTGGCCATGTc                                                          <  8:238138/146‑1 (MQ=255)
                                                                                       tGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGgcgc                                                                                                                             >  6:95320/1‑80 (MQ=255)
                                                                                       tGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGgcgc                                                                                                                             <  5:95320/80‑1 (MQ=255)
                                                                                       tGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGcc                                                               <  8:312006/142‑1 (MQ=255)
                                                                                       tGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGcc                                                               >  7:312006/1‑142 (MQ=255)
                                                                                           cTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTgg                                                                 >  5:198836/1‑136 (MQ=255)
                                                                                           cTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTgg                                                                 <  6:198836/136‑1 (MQ=255)
                                                                                             ggAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGa                                                        >  7:94176/1‑143 (MQ=255)
                                                                                             ggAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGa                                                        <  8:94176/143‑1 (MQ=255)
                                                                                               aaaTCGTCCGGGCCCAAGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAg                                                <  4:184212/149‑1 (MQ=255)
                                                                                                  tgctCCGGGCCATAGGCCAGGCCGTGGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGc                                              >  4:195671/4‑148 (MQ=255)
                                                                                                    gTCCGTGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGaaa                                                      <  3:127995/138‑1 (MQ=255)
                                                                                                    gTCCGGGCCCAAGCCCTGGCCTATGGCCCGTAGCAGAGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATAGCCAGCGAGCTGTTTGCGGATGCGTGATCAAGCGACCGGTTGGTTGCAAGCGTGGCCCTGTCGaaa                                                      >  4:127995/1‑138 (MQ=255)
                                                                                                            ccACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGc                                                                     >  6:262831/1‑115 (MQ=255)
                                                                                                            ccACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGc                                                                     <  7:265788/115‑1 (MQ=255)
                                                                                                            ccACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGc                                                                     >  8:265788/1‑115 (MQ=255)
                                                                                                            ccACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGc                                                                     <  5:262831/115‑1 (MQ=255)
                                                                                                            ccACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCt                                   >  5:252047/1‑149 (MQ=255)
                                                                                                              aCGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTg                                  >  2:209712/1‑148 (MQ=255)
                                                                                                                ggccttgcctATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAg                               <  6:252047/149‑1 (MQ=255)
                                                                                                                  cctggcctATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGaa                             <  4:300473/149‑1 (MQ=255)
                                                                                                                          aTGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTg                                                            <  5:256136/110‑1 (MQ=255)
                                                                                                                          aTGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTg                                                            >  6:256136/1‑110 (MQ=255)
                                                                                                                                     cagcgcagcCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGTTCAGGCCACCGGGAt          <  7:178712/149‑1 (MQ=255)
                                                                                                                                        cgcagcCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGTTCAg                     >  5:258091/1‑135 (MQ=255)
                                                                                                                                        cgcagcCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGTTCAg                     <  6:258091/135‑1 (MQ=255)
                                                                                                                                        cgcagcCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGTTCAGGCCACCGGGATGAg       >  7:89305/1‑149 (MQ=255)
                                                                                                                                         gcagcCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGTTCAGGCCACCGGGATGAGc      <  5:111101/149‑1 (MQ=255)
                                                                                                                                             ccAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGt                                 <  4:273524/118‑1 (MQ=255)
                                                                                                                                             ccAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGt                                 >  3:273524/1‑118 (MQ=255)
                                                                                                                                               aGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGTTCAGGCCACCGGGATGAGCGGAt  >  7:306637/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGTTCAGGCCACCGGGATGAGCGGAT  >  NC_002947/3267943‑3268232

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: