Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,529,190 C→A 7.3% intergenic (‑404/+190) 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,1900CA 7.3% 199.1 / 4.1 63intergenic (‑404/+190)ydfG/livF‑IIbifunctional NADP‑dependent 3‑hydroxy acid dehydrogenase/malonic semialdehyde reductase/ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (32/26);  new base A (3/2);  total (35/28)
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.00e+00

TGCCGCCTGCAGCGGTCTGTCCAAGGTTTTTACACAGAGTTATCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTC  >  NC_002947/5529047‑5529338
                                                                                                                                               |                                                                                                                                                         
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                                                               >  4:151321/1‑146 (MQ=255)
                                                                                     gAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGACGGCTCCACACTATCCGTGGAAAAGAGTTGGCCtgt                                                               <  3:151321/146‑1 (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                                 <  4:177420/149‑1 (MQ=255)
                                                                                                                gggATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGgt                                 <  6:32523/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                                   <  1:49356/138‑1 (MQ=255)
                                                                                                                         aCAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATgg                                   >  2:49356/1‑138 (MQ=255)
                                                                                                                         aCAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCg                         <  4:138546/148‑1 (MQ=255)
                                                                                                                         aCAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCg                         <  7:247289/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                                                                                     >  6:327272/1‑80 (MQ=255)
                                                                                                                                 tGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTAt                                                                                     <  5:327272/80‑1 (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                                                             >  1:248038/1‑98 (MQ=255)
                                                                                                                                       aaTCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTgg                                                             <  2:248038/98‑1 (MQ=255)
                                                                                                                                        aTCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAggg                <  4:171860/142‑1 (MQ=255)
                                                                                                                                        aTCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAggg                >  3:171860/1‑142 (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)
                                                                                                                                               |                                                                                                                                                         
TGCCGCCTGCAGCGGTCTGTCCAAGGTTTTTACACAGAGTTATCCACAGGATTGGCACTGCACTTGAAGCACACGCAGAGATGTTGAAAGCTTGGGAAAATCATCCACAACGGGGATATGAACAAAAAATGATCAAATCGCTACAGGGCTTGATACTCAAGGCCTTG‑‑‑‑AGCGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTC  >  NC_002947/5529047‑5529338

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: