Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,529,360 C→G 6.0% intergenic (‑574/+20) 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,3600CG 6.0% 197.9 / 4.1 67intergenic (‑574/+20)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 (33/29);  new base G (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

CGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGACAGTTTCAGCGCATGGTTG  >  NC_002947/5529216‑5529506
                                                                                                                                                |                                                                                                                                                  
cGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGctgcactgca                                                                                                                                                <  3:94536/149‑1 (MQ=255)
  atatTCCACCATGTTTTCCAGAGGCGGGTCCACTCTTTCCGTGGAAAAGAGTTGGCCTGTGGAAAGTATGGGTTGCGAATGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCAGATTGAGAAGGGGGCTGCAGTGCAGc                                                                                                                                              <  4:304981/149‑1 (MQ=255)
     tGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGccc                                                                                                                                            <  6:218859/148‑1 (MQ=255)
     tGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGccc                                                                                                                                            <  4:180757/148‑1 (MQ=255)
      gCCACCATGTTTTCCCCAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGGGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCt                                                                                                                                          >  7:319358/1‑149 (MQ=255)
      gCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCt                                                                                                                                          >  1:164255/1‑149 (MQ=255)
        caccaTGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAAGGGGCTGCAGTGCAGccccctt                                                                                                                                        <  3:97564/149‑4 (MQ=255)
        caccaTGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAAGGGGCTGCAGTGCAGccccct                                                                                                                                         >  4:97564/1‑146 (MQ=255)
         accaTGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTaa                                                                                                                                       >  5:118999/1‑149 (MQ=255)
         accaTGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTaa                                                                                                                                       <  4:292278/149‑1 (MQ=255)
           caTGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAaca                                                                                                                                     >  8:62785/1‑149 (MQ=255)
           caTGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAaca                                                                                                                                     <  8:192642/149‑1 (MQ=255)
           caTGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAaca                                                                                                                                     >  4:184053/1‑149 (MQ=255)
                   ccACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAg                                                                                                                             >  2:118721/1‑149 (MQ=255)
                   ccACAGGCGGCTCCACAATATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCa                                                                                                                              <  7:192516/148‑1 (MQ=255)
                    cacaGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGTGTGTGTTTCGGGAGAAGGGAACGAAAAGTTTTCAACATTGAGAAGGGGCCTGAAAAGCAGCCCCTTAACACAACtaatt                                                                                                                            >  6:251530/1‑145 (MQ=255)
                        ggcggcTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTgccc                                                                                                                        <  8:371666/149‑1 (MQ=255)
                                 acacTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGc                                                                                                                                              >  5:83383/1‑118 (MQ=255)
                                 acacTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGc                                                                                                                                              <  6:83383/118‑1 (MQ=255)
                                   acTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAg                                                                                                             >  4:288737/1‑149 (MQ=255)
                                    cTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGc                                                                                                            <  1:118721/149‑1 (MQ=255)
                                    cTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGc                                                                                                            <  8:231189/149‑1 (MQ=255)
                                       tCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGtt                                                                                                         <  3:184053/149‑1 (MQ=255)
                                         cGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTa                                                                                                        <  6:118999/148‑1 (MQ=255)
                                         cGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTAc                                                                                                       <  5:330292/149‑1 (MQ=255)
                                                    gTTGGCCTGTGGAAACTATGGCTTGCGAAGGAATTTGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGGCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCa                                                                                            <  7:292515/149‑1 (MQ=255)
                                                      tGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGGGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACTACTCAGTGCCGGCCCAGATAAGCGTTACGCACTTCCTCAtt                                                                                          <  7:168664/149‑1 (MQ=255)
                                                       ggCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCAtt                                                                                          <  6:309804/148‑1 (MQ=255)
                                                       ggCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCAtt                                                                                          <  7:62785/148‑1 (MQ=255)
                                                        gCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGa                                                                                        >  6:107066/1‑149 (MQ=255)
                                                             tgGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGACCAGATAAGCGTTACGCACTTCCTCATTGCCCaa                                                                                    >  7:303247/1‑148 (MQ=255)
                                                               gAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCCACAg                                                                                 >  8:260311/1‑149 (MQ=255)
                                                               gAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACa                                                                                  <  5:107066/148‑1 (MQ=255)
                                                               gAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAACGGGGCTGCACTGCATCCCCTTAACACAACTCAGTGCCCGCCCAGATAACCGTTACGCACTTCATCTTTGACAAACAg                                                                                 >  5:319185/1‑149 (MQ=255)
                                                                 aaCAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGt                                                                                >  8:244634/1‑148 (MQ=255)
                                                                    aaTGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGTCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCt                                                                            <  8:110811/149‑1 (MQ=255)
                                                                         cTTGCGAAGGAATATGGGTGTGTTGCGGGAGAAGGGAACGACAAGGTTTCCCCATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGcccg                                                                       <  8:319358/149‑1 (MQ=255)
                                                                          ttGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGcccgg                                                                      >  4:385869/1‑149 (MQ=255)
                                                                             cGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTgcc                                                                                                                         <  4:61418/95‑1 (MQ=255)
                                                                             cGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTgcc                                                                                                                         >  3:61418/1‑95 (MQ=255)
                                                                              gAAGGAATATGGGTGTGTTTCGGGAGAATGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGcccggtccc                                                                  <  5:251530/149‑1 (MQ=255)
                                                                               aaGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGcccggtcccg                                                                 <  2:361852/149‑1 (MQ=255)
                                                                                     tatGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATg                                                           <  3:288737/149‑1 (MQ=255)
                                                                                           tgtgtTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATc                                                     <  3:358864/149‑1 (MQ=255)
                                                                                            gtgtTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCATCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCt                                                    >  6:12677/1‑149 (MQ=255)
                                                                                            gtgtTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCt                                                    <  2:164255/149‑1 (MQ=255)
                                                                                               tttCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGcc                                                 >  3:217831/1‑149 (MQ=255)
                                                                                               tttCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCGTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGcc                                                 >  4:341015/1‑149 (MQ=255)
                                                                                                    ggAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTg                                            >  1:175588/1‑149 (MQ=255)
                                                                                                            gAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATca                                     >  7:48578/1‑148 (MQ=255)
                                                                                                                 aCAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGc                                                          <  5:350913/122‑1 (MQ=255)
                                                                                                                 aCAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGc                                                          >  6:350913/1‑122 (MQ=255)
                                                                                                                   aaGTTTTCCTGATTGAGTAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGcc                             <  6:319185/149‑1 (MQ=255)
                                                                                                                   aaGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGcc                             >  3:394336/1‑149 (MQ=255)
                                                                                                                      ttttCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCgg                          >  6:124736/1‑149 (MQ=255)
                                                                                                                           cacaTTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGtt                                             >  7:393103/1‑125 (MQ=255)
                                                                                                                           cacaTTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGtt                                             <  8:393103/125‑1 (MQ=255)
                                                                                                                            acaTTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCTCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCAAATAGGCCCGGTCCGa                     >  7:407413/1‑148 (MQ=255)
                                                                                                                            acaTTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGa                     >  1:28396/1‑148 (MQ=255)
                                                                                                                             caTTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGACa                   >  1:241769/1‑149 (MQ=255)
                                                                                                                              aTTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGCCCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGACAg                  >  7:106552/1‑149 (MQ=255)
                                                                                                                                 gagaAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGa                                           <  3:160630/121‑1 (MQ=255)
                                                                                                                                 gagaAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGa                                           >  4:160630/1‑121 (MQ=255)
                                                                                                                                    aaGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGACAGTTTCAg            >  6:211425/1‑149 (MQ=255)
                                                                                                                                      gggggCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCACCAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCGGGCAGTTTCAgcg          >  4:315158/1‑149 (MQ=255)
                                                                                                                                       ggggCTGCAGTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGACAGTTTCAgcg          >  5:16481/1‑148 (MQ=255)
                                                                                                                                         ggCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGACAGTTTCAGCGCAt       >  5:78281/1‑149 (MQ=255)
                                                                                                                                          gCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGACAGTTTCAGCGCATg      >  1:94372/1‑149 (MQ=255)
                                                                                                                                             gcactgcaGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGACAGTTTCAGCGCATGGtt   >  5:210918/1‑149 (MQ=255)
                                                                                                                                              cactgcaGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGACAGTTTCAGCGCATGGTTg  <  8:48578/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CGATATGCCACCATGTTTTCCACAGGCGGCTCCACACTATCCGTGGAAAAGAGTTGGCCTGTGGAAACAATGGCTTGCGAAGGAATATGGGTGTGTTTCGGGAGAAGGGAACGACAAGTTTTCCACATTGAGAAGGGGGCTGCACTGCAGCCCCTTAACACAACTCAGTGCCCGCCCAGATAAGCGTTACGCACTTCCTCATTGACCAACAGTTCCTGCCCGGTCCCGGTCATGCGAATCTGCCCGTTGACCATCACATAGGCCCGGTCCGACAGTTTCAGCGCATGGTTG  >  NC_002947/5529216‑5529506

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: