Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,663,901 T→C T240A (ACC→GCC)  gadW ← transcriptional activator of gadA and gadBC; repressor of gadX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,663,9010TC99.5% 749.7 / ‑4.8 215T240A (ACC→GCC) gadWtranscriptional activator of gadA and gadBC; repressor of gadX
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (108/106);  total (109/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 215 total aligned reads displayed.

AGCCGGTAGCCCGGCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGTACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCCACATTTTTC  >  NC_000913/3663804‑3663997
                                                                                                 |                                                                                                
aGCCGGTAGCCCGGCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCAc                                                                                                <  2:1124844/100‑1 (MQ=255)
 gCCGGTAGCCCGGCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCAcc                                                                                               >  1:4559122/1‑100 (MQ=255)
  ccGGTAGCCCGGCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCt                                                                                              >  1:1908980/1‑100 (MQ=255)
    ggTAGCCCGGCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTgg                                                                                            <  2:2798434/100‑1 (MQ=255)
    ggTAGCCCGGCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTgg                                                                                            >  1:3927195/1‑100 (MQ=255)
        gCCCGGCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGaa                                                                                        >  2:4679620/1‑100 (MQ=255)
        gCCCGGCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGaa                                                                                        <  1:4928643/100‑1 (MQ=255)
         cccGGCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAAt                                                                                       <  1:4014685/100‑1 (MQ=255)
             gCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTg                                                                                   >  2:2450458/1‑100 (MQ=255)
              cAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTgc                                                                                  <  2:4920749/100‑1 (MQ=255)
              cAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTgc                                                                                  <  2:1228132/100‑1 (MQ=255)
                 aaaTCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCa                                                                               <  1:1610255/100‑1 (MQ=255)
                    tcatcaGGGACGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCCGTGTCCACATGCACCCTATGTCAGGGAAAGGCCCCTGGCGTATGCTgcgcgacc                                                                            >  2:2724818/1‑96 (MQ=255)
                    tcatcaGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAAc                                                                            >  1:4936654/1‑100 (MQ=255)
                        caGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGAt                                                                        >  2:894346/1‑100 (MQ=255)
                         aGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGAtg                                                                       <  2:4776557/100‑1 (MQ=255)
                            gAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTgg                                                                    >  2:1990332/1‑100 (MQ=255)
                             aaGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGc                                                                   <  2:3277264/100‑1 (MQ=255)
                               gagTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGt                                                                 <  2:319996/100‑1 (MQ=255)
                                 gTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTa                                                               >  2:819656/1‑100 (MQ=255)
                                 gTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTa                                                               >  1:2796930/1‑100 (MQ=255)
                                   ttCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTaca                                                             >  1:1420704/1‑100 (MQ=255)
                                    tCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTacac                                                            >  2:552911/1‑100 (MQ=255)
                                     cacaTGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTCCAcc                                                           >  2:4831778/1‑100 (MQ=255)
                                     cacaTGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACAcc                                                           >  2:1451376/1‑100 (MQ=255)
                                     cacaTGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACAcc                                                           >  1:2353107/1‑100 (MQ=255)
                                     cacaTGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTAAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACAcc                                                           >  1:2870005/1‑100 (MQ=255)
                                       caTGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCat                                                         <  1:754645/100‑1 (MQ=255)
                                       caTGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCat                                                         >  1:3414245/1‑100 (MQ=255)
                                       caTGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCat                                                         >  2:741413/1‑100 (MQ=255)
                                        aTGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTGAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCata                                                        >  2:4375804/1‑100 (MQ=255)
                                         tGAAGCAGGTGTGAGATCCTGGCCAAGATTCACATGCGAAAGATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAAATGATGTGGCGTTACACCataa                                                       <  1:2062589/100‑1 (MQ=255)
                                         tGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCataa                                                       <  2:3655038/100‑1 (MQ=255)
                                          gAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCataat                                                      <  2:4468512/100‑1 (MQ=255)
                                          gAAGCAGGTGAGAGATCCTGACCAATATACAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCataat                                                      <  2:3004984/100‑1 (MQ=255)
                                           aaGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCataata                                                     <  1:1029058/100‑1 (MQ=255)
                                                ggTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTc                                                >  2:773510/1‑100 (MQ=255)
                                                ggTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTc                                                >  1:795224/1‑100 (MQ=255)
                                                cgtgtGAGATCCTGACCAATATTCAAATGCGAAAAAAAAAAAAAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTc                                                <  2:81462/99‑1 (MQ=255)
                                                 gtgtGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACTTGCCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTATCg                                               >  2:4002994/1‑100 (MQ=255)
                                                 gtgtGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCg                                               >  2:2952114/1‑100 (MQ=255)
                                                 gtgtGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCg                                               <  2:2765585/100‑1 (MQ=255)
                                                 gtgtGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCg                                               <  1:399318/100‑1 (MQ=255)
                                                 gtgtGAGAGCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCg                                               <  2:2035906/100‑1 (MQ=255)
                                                  tgtgAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGa                                              <  2:1341254/100‑1 (MQ=255)
                                                   gtgAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGaa                                             >  2:3050052/1‑100 (MQ=255)
                                                   gtgAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGaa                                             >  2:1622746/1‑100 (MQ=255)
                                                    tgAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGaaa                                            <  2:2330100/100‑1 (MQ=255)
                                                     gagaTCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAAt                                           >  1:3552124/1‑100 (MQ=255)
                                                      agaTCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAAtg                                          <  1:2985622/100‑1 (MQ=255)
                                                      agaTCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGGGGCGTTACACCATAATATTGTCGAAAtg                                          <  2:2202977/100‑1 (MQ=255)
                                                         tCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATgtgt                                       >  2:3385906/1‑100 (MQ=255)
                                                          ccTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGtt                                      <  2:3454732/100‑1 (MQ=255)
                                                             gACCAATATTCAAATGCGAAATATGTCAGGAAAAGGTACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTtat                                   >  1:99572/1‑100 (MQ=255)
                                                              aCCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACCCCATAATATTGTCGAAATGTGTTtata                                  <  2:2393745/100‑1 (MQ=255)
                                                              aCCAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTtata                                  <  1:1127207/100‑1 (MQ=255)
                                                               tcAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATaa                                 <  2:2660826/99‑1 (MQ=255)
                                                               ccAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATaa                                 >  2:141503/1‑100 (MQ=255)
                                                               ccAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATaa                                 <  1:502322/100‑1 (MQ=255)
                                                                cAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATaaa                                >  1:3439199/1‑100 (MQ=255)
                                                                cAATATTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATaaa                                >  1:2691361/1‑100 (MQ=255)
                                                                  atatTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGt                              >  1:3289193/1‑100 (MQ=255)
                                                                  atatTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGt                              <  2:1819263/100‑1 (MQ=255)
                                                                  atatTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGt                              <  2:3713186/100‑1 (MQ=255)
                                                                   tatTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTa                             <  2:1582877/100‑1 (MQ=255)
                                                                    nttCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTAc                            >  1:4831374/2‑100 (MQ=255)
                                                                    atTCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTAc                            >  1:2226124/1‑100 (MQ=255)
                                                                     ttCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACg                           >  1:4702327/1‑100 (MQ=255)
                                                                     ttCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACg                           <  1:143265/100‑1 (MQ=255)
                                                                     ttCAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACg                           >  1:1344038/1‑100 (MQ=255)
                                                                       naaaTGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGAt                         >  1:516427/2‑100 (MQ=255)
                                                                       cAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGAt                         >  2:1940418/1‑100 (MQ=255)
                                                                       cAAATGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGAt                         <  2:3876928/100‑1 (MQ=255)
                                                                        aaaTGCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATg                        >  2:3149977/1‑100 (MQ=255)
                                                                            gCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACt                    <  1:1676134/100‑1 (MQ=255)
                                                                            gCGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACt                    <  1:963929/100‑1 (MQ=255)
                                                                             cGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTg                   <  1:3193555/100‑1 (MQ=255)
                                                                             cGAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTAAACCATACTATTGTCGAAATGTGTTTATAAAGTACGATGTACTg                   <  1:263579/100‑1 (MQ=255)
                                                                              gAAATATGTCAGGAAAAGGCGCCTGGCGCATGTTGCGCAAAGTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGc                  <  2:1228966/100‑1 (MQ=255)
                                                                              gAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGc                  <  2:3513154/100‑1 (MQ=255)
                                                                              gAAATATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGc                  <  2:1719765/100‑1 (MQ=255)
                                                                               aaaTATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCt                 >  1:215209/1‑100 (MQ=255)
                                                                               aaaTATGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCt                 >  1:598975/1‑100 (MQ=255)
                                                                               aaaTATGACATCAAAATGCACCCGGCAGATGGTGCGCAAGCGGATGTGTCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCt                 <  1:1274455/100‑1 (MQ=255)
                                                                                  tatGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCtata              >  1:3826545/1‑100 (MQ=255)
                                                                                   atGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGcagatc             >  1:725836/1‑95 (MQ=255)
                                                                                   atGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAg             >  1:4603457/1‑100 (MQ=255)
                                                                                    tGTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGc            <  1:282677/100‑1 (MQ=255)
                                                                                     ttCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGcc           <  1:3563963/99‑1 (MQ=255)
                                                                                     gTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGcc           <  1:4046740/100‑1 (MQ=255)
                                                                                     gTCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGcc           >  1:10745/1‑100 (MQ=255)
                                                                                      tCAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCca          >  2:3524895/1‑100 (MQ=255)
                                                                                       cAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCcac         >  2:3395698/1‑100 (MQ=255)
                                                                                       cAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCcac         <  1:1770515/100‑1 (MQ=255)
                                                                                       cAGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCcac         >  1:4276698/1‑100 (MQ=255)
                                                                                        aGGAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCcaca        <  1:4428210/100‑1 (MQ=255)
                                                                                         ggAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCCACAt       >  2:4584767/1‑100 (MQ=255)
                                                                                         ggAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCCACAt       >  2:3342302/1‑100 (MQ=255)
                                                                                          gAAAAGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCCACAtt      <  1:1967006/100‑1 (MQ=255)
                                                                                           aaaaGGCACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCCACAttt     <  1:3147347/100‑1 (MQ=255)
                                                                                                 |                                                                                                
AGCCGGTAGCCCGGCAGAAATCATCAGGGAAGAGTTTCACATGAAGCAGGTGTGAGATCCTGACCAATATTCAAATGCGAAATATGTCAGGAAAAGGTACCTGGCGAATGTTGCGCAAACTGATGTGGCGTTACACCATAATATTGTCGAAATGTGTTTATAAAGTACGATGTACTGCTATAGCCACATTTTTC  >  NC_000913/3663804‑3663997

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: