Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,529,194 G→C 6.0% intergenic (‑408/+186) ydfG ← / ← livF‑II bifunctional NADP‑dependent 3‑hydroxy acid dehydrogenase/malonic semialdehyde reductase/ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,529,1940GC 6.0% 204.1 / 4.3 67intergenic (‑408/+186)ydfG/livF‑IIbifunctional NADP‑dependent 3‑hydroxy acid dehydrogenase/malonic semialdehyde reductase/ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base C (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 = 5.68e-01

TGCCGCCTGCAGCGGTCTGTCCAAGGTTTTTACACAGAGTTATCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACA  >  NC_002947/5529047‑5529342
                                                                                                                                                   |                                                                                                                                                         
tGCCGCATGCAGCGGTATGTCCAATGTTTTTACAGAGAGTTATCCACATGATTGTCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGc                                                                                                                                                          <  1:110639/149‑1 (MQ=255)
                   tCCAAGGTTTTTACACAGAGTTATCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCtt                                                                                                                                         >  6:62697/1‑147 (MQ=255)
                     cAAGGTTTTTACACAGAGTTATCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGc                                                                                                                                >  8:212311/1‑149 (MQ=255)
                      aaGGTTTTTACACAGAGTTATCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGc                                                                                                                                >  6:85621/1‑148 (MQ=255)
                       aGGTTTTTACACAGACTTATCCACAGGATTGGCACTGCACTTGCAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGa                                                                                                                              <  4:114986/149‑1 (MQ=255)
                             ttAAACAGAGTTATCGAGAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTTGGATAATCATCCACAACGGTGATATGAACAAAAATTGATCAAATCGCTCAAGGCCTTGAGTATCAAGCCCTGT‑‑‑‑AGCGATTTga                                                                                                                         <  4:352134/148‑2 (MQ=21)
                               acacAGAGTTATCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTCAAGGCCTTGAGTATCAAGCCCTGT‑‑‑‑AGCGATTTgatc                                                                                                                       <  6:297480/148‑4 (MQ=37)
                                       ttATCCACAGGATTGGCACTGCACATGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCAACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGtttt                                                                                                              >  8:55768/1‑149 (MQ=255)
                                         aTCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTc                                                                                                             <  7:212311/148‑1 (MQ=255)
                                         aTCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTc                                                                                                             >  5:9747/1‑148 (MQ=255)
                                          tCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTcc                                                                                                            <  5:85621/148‑1 (MQ=255)
                                                ggATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGTAAAATCATCCACAACGGGTATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAggc                                                                                                      <  2:179520/148‑1 (MQ=255)
                                                   ttGGAACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCt                                                                                                  <  5:62697/149‑1 (MQ=255)
                                                     ggCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACGATGTTTTCCACAGGCGGCTc                                                                                                 <  3:36970/148‑1 (MQ=255)
                                                        actgcactTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCcaca                                                                                             >  2:300515/1‑149 (MQ=255)
                                                         ctgcactTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCacac                                                                                            >  5:32523/1‑149 (MQ=255)
                                                         ctgcactTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCacac                                                                                            >  3:177420/1‑149 (MQ=255)
                                                          tgcactTGAAGCACACGTAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCccac                                                                                            >  8:253724/1‑148 (MQ=255)
                                                            cacaTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACt                                                                                           <  7:55768/147‑1 (MQ=255)
                                                                        acGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATAGCTAAAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCAACATGTTTTCCACAGGCGGCTCCACACTATCCGTGGaa                                                                                 >  1:266575/1‑145 (MQ=255)
                                                                          gCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGt                                                                           >  6:118082/1‑149 (MQ=255)
                                                                          gCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGt                                                                           >  4:74612/1‑149 (MQ=255)
                                                                          gCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGt                                                                           >  1:184475/1‑149 (MQ=255)
                                                                            agagATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTg                                                                         >  3:175025/1‑149 (MQ=255)
                                                                             gagaTGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTgg                                                                        <  4:339090/149‑1 (MQ=255)
                                                                                aTGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCt                                                                     <  5:118082/149‑1 (MQ=255)
                                                                                    tGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCtgtg                                                                  >  3:73016/1‑148 (MQ=255)
                                                                                     gAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGACGGCTCCACACTATCCGTGGAAAAGAGTTGGCCtgt                                                                   <  3:151321/146‑1 (MQ=255)
                                                                                     gAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGACGGCTCCACACTATCCGTGGAAAAGAGTTGGCCtgt                                                                   >  4:151321/1‑146 (MQ=255)
                                                                                              ggAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATgg                                                        <  4:73016/148‑1 (MQ=255)
                                                                                                aaaaTCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCtt                                                     >  8:247289/1‑149 (MQ=255)
                                                                                                 aaaTCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTg                                                    <  8:186800/149‑1 (MQ=255)
                                                                                                      atcCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGaa                                                <  1:181611/148‑1 (MQ=255)
                                                                                                      atcCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAg                                               <  8:34010/149‑1 (MQ=255)
                                                                                                         cacaACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGaa                                            >  3:138546/1‑149 (MQ=255)
                                                                                                          acaACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCtgtg                                                                  >  8:102128/1‑126 (MQ=255)
                                                                                                          acaACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCtgtg                                                                  <  7:102128/126‑1 (MQ=255)
                                                                                                          acaACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGaa                                            >  7:42017/1‑148 (MQ=255)
                                                                                                             aCGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATg                                        <  8:42017/149‑1 (MQ=255)
                                                                                                              cGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATgg                                       >  7:126076/1‑149 (MQ=255)
                                                                                                                gggATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGgt                                     <  6:32523/149‑1 (MQ=255)
                                                                                                                gggATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGgt                                     <  4:177420/149‑1 (MQ=255)
                                                                                                                       gAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTc                              >  3:225470/1‑149 (MQ=255)
                                                                                                                        aaCAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCg                             >  2:292815/1‑149 (MQ=255)
                                                                                                                        aaCAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCg                             >  8:181986/1‑149 (MQ=255)
                                                                                                                         aCAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATgg                                       >  2:49356/1‑138 (MQ=255)
                                                                                                                         aCAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATgg                                       <  1:49356/138‑1 (MQ=255)
                                                                                                                         aCAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCg                             <  7:247289/148‑1 (MQ=255)
                                                                                                                         aCAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCg                             <  4:138546/148‑1 (MQ=255)
                                                                                                                          cAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCgg                            >  2:325375/1‑148 (MQ=255)
                                                                                                                          cAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCgg                            >  7:339083/1‑148 (MQ=255)
                                                                                                                            aaaaaTTATCAAATCGCTTCATGTCTTGTTACTAAAGGCCTTGTCTAACCCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTc                              <  1:300515/149‑1 (MQ=255)
                                                                                                                                 tGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTAt                                                                                         <  5:327272/80‑1 (MQ=255)
                                                                                                                                 tGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTAt                                                                                         >  6:327272/1‑80 (MQ=255)
                                                                                                                                 tGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAgg                     >  6:6000/1‑148 (MQ=255)
                                                                                                                                  gATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAggg                    >  2:42169/1‑148 (MQ=255)
                                                                                                                                      ataTCGCTCAAGGCCTTGAGTATCAAGCCCTGT‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGTGTGTGTTTCGGGAGAAGGGAAc                 >  6:250027/3‑147 (MQ=37)
                                                                                                                                      ataTCGCTCAAGGCCTTGAGTATCAAGCCCTGT‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGa               >  3:2622/3‑149 (MQ=38)
                                                                                                                                       aaTCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTgg                                                                 <  2:248038/98‑1 (MQ=255)
                                                                                                                                       aaTCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTgg                                                                 >  1:248038/1‑98 (MQ=255)
                                                                                                                                        aTCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAggg                    >  3:171860/1‑142 (MQ=255)
                                                                                                                                        aTCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAggg                    <  4:171860/142‑1 (MQ=255)
                                                                                                                                             tACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGttt        >  7:109121/1‑149 (MQ=255)
                                                                                                                                               cAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGtttt       <  7:253724/148‑1 (MQ=255)
                                                                                                                                               cAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTc      <  1:42169/149‑1 (MQ=255)
                                                                                                                                               cAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTc      <  3:74612/149‑1 (MQ=255)
                                                                                                                                                aGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTcc     <  5:6000/149‑1 (MQ=255)
                                                                                                                                                 gggCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCca    <  4:175025/149‑1 (MQ=255)
                                                                                                                                                 gggCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCca    <  1:325375/149‑1 (MQ=255)
                                                                                                                                                   gCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCcaca  >  7:319627/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                         
TGCCGCCTGCAGCGGTCTGTCCAAGGTTTTTACACAGAGTTATCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACA  >  NC_002947/5529047‑5529342

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: