Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,977,224:1 +A 16.4% intergenic (‑85/+42) motA ← / ← flhC proton conductor component of flagella motor/flagellar class II regulon transcriptional activator, with FlhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,977,2231.A16.4% 469.3 / 103.6 146intergenic (‑84/+43)motA/flhCproton conductor component of flagella motor/flagellar class II regulon transcriptional activator, with FlhD
Reads supporting (aligned to +/- strand):  ref base . (63/59);  new base A (16/8);  total (79/67)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

Only 100 of 150 total aligned reads displayed.

CTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATT  >  NC_000913/1977124‑1977322
                                                                                                    |                                                                                                   
cTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGt                                                                                                      >  1:2518832/1‑100 (MQ=255)
 taataaGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑a                                                                                                    >  1:899311/1‑100 (MQ=255)
      agataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgtcg                                                                                               >  1:282388/1‑100 (MQ=255)
         taagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTa                                                                                            >  2:1354306/1‑100 (MQ=255)
             cACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCGAAAAAAAGCAGCGTT‑ACGTCGTTACCgc                                                                                        >  1:1033879/1‑100 (MQ=255)
             cACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgc                                                                                        <  2:3016533/100‑1 (MQ=255)
               cGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgctg                                                                                      >  1:2085104/1‑100 (MQ=255)
                gACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgc                                                                                     <  2:2250096/100‑1 (MQ=255)
                 tccgatctcTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgct                                                                                    <  2:1507690/92‑1 (MQ=37)
                 aCATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgct                                                                                    <  1:2841331/100‑1 (MQ=38)
                 aCATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgct                                                                                    <  2:1020729/100‑1 (MQ=38)
                 aCATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgctgct                                                                                    <  2:2535310/100‑1 (MQ=255)
                  catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgctt                                                                                   <  2:12172/100‑2 (MQ=38)
                  catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGctgctg                                                                                   <  1:2952879/100‑1 (MQ=255)
                  catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGctgctg                                                                                   <  2:1271918/100‑1 (MQ=255)
                  catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑AAGTCGTTACCGctgctg                                                                                   >  1:2261400/1‑100 (MQ=255)
                    tcatcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGa                                                                                 >  2:2851920/1‑100 (MQ=255)
                     catcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgcttttt                                                                                <  2:651327/100‑5 (MQ=38)
                     catcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGaa                                                                                >  1:2962217/1‑100 (MQ=255)
                      atcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAAt                                                                               <  2:2916714/100‑1 (MQ=255)
                          ttCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgcttttttttgc                                                                           <  2:923578/100‑10 (MQ=25)
                           tCCACTGTTGACAATGACAGGATCTTCAGTCGTCAGTCGTTAAAGAGCGATTGGATTAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTgc                                                                          <  1:1592680/100‑1 (MQ=255)
                            ccACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTgcg                                                                         >  2:437820/1‑100 (MQ=255)
                            ccACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTgcg                                                                         <  2:76361/100‑1 (MQ=255)
                             cACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTgcgc                                                                        <  2:2166000/100‑1 (MQ=255)
                                 gTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCa                                                                    >  2:2744651/1‑100 (MQ=255)
                                 gTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCa                                                                    >  2:3310897/1‑100 (MQ=255)
                                  ttGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAc                                                                   <  2:766750/100‑1 (MQ=255)
                                  ttGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAc                                                                   <  1:2675729/100‑1 (MQ=255)
                                   tGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAcc                                                                  >  2:949362/1‑100 (MQ=255)
                                   tGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAcc                                                                  >  1:3411239/1‑100 (MQ=255)
                                    gACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg                                                                 <  2:1480404/100‑1 (MQ=255)
                                    gACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg                                                                 >  2:3157309/1‑100 (MQ=255)
                                    gACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg                                                                 >  2:2319779/1‑100 (MQ=255)
                                     aCCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                                                >  1:3592223/1‑100 (MQ=255)
                                     aCCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                                                >  1:1168837/1‑100 (MQ=255)
                                      ccATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTa                                                               <  1:2561970/100‑1 (MQ=255)
                                        aTGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATc                                                             <  2:5039/100‑1 (MQ=255)
                                         tGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa                                                            <  1:1395056/100‑1 (MQ=255)
                                          ggcAGGATGTTCAGTAGTCAGGCGTTAACGCCCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAg                                                           <  2:2261415/98‑1 (MQ=255)
                                          gACAGGATGTCCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAg                                                           <  2:2085757/100‑1 (MQ=255)
                                           aCAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGt                                                          >  1:848846/1‑100 (MQ=255)
                                            cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                         >  1:1122250/1‑100 (MQ=255)
                                            cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                         <  1:3133319/100‑1 (MQ=255)
                                            cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                         >  1:3182033/1‑100 (MQ=255)
                                            cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                         >  1:1505945/1‑100 (MQ=255)
                                             aGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTa                                                        <  2:872930/100‑1 (MQ=255)
                                              ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                       <  2:1782085/100‑1 (MQ=255)
                                               gATGTTCAGTCGTCAGGCGTTAACGCGAAAAAAAAAAAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaaa                                                      <  1:1248970/100‑1 (MQ=255)
                                                aTGTTCAGTCGTCAGGCGTTAACGCGAAATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAAc                                                     <  1:1771723/100‑1 (MQ=255)
                                                  gTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAg                                                   >  1:1487451/1‑100 (MQ=255)
                                                  gTTCAGTCGGCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAg                                                   <  1:2714063/100‑1 (MQ=255)
                                                   ttCAGTCGTCAGGCGTTAACGCGCGATTGAGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGc                                                  <  1:541359/100‑1 (MQ=255)
                                                    tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                 >  2:1435377/1‑100 (MQ=255)
                                                     cAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                                <  2:3325046/100‑1 (MQ=255)
                                                       gtcgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGt                                              >  1:3454297/1‑100 (MQ=255)
                                                        tcgtcAGGCGTTCACTCGTAAATGTGATCAAAAAAAGCAGCGTT‑ACGTAGTTAGCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                                <  1:3445029/97‑1 (MQ=255)
                                                         cgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAc                                            >  1:3720183/1‑100 (MQ=255)
                                                         cgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAc                                            >  2:1046300/1‑100 (MQ=255)
                                                          gtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAct                                           <  2:3597021/100‑1 (MQ=255)
                                                          gtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTACCAGTTAAACAGCCTGTAct                                           >  1:1012545/1‑100 (MQ=255)
                                                             aGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctc                                        <  2:2869101/100‑1 (MQ=255)
                                                               gCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTg                                      <  2:349802/100‑1 (MQ=255)
                                                                cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCGGt                                     >  2:1535177/1‑100 (MQ=255)
                                                                 gTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                    <  2:1166948/100‑1 (MQ=255)
                                                                    aaCGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCAt                                 <  2:3411267/100‑1 (MQ=255)
                                                                      cgcgcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCTCTGCTGGAATGTTTCGCCTCACCGTATCAGTTAACCCGCCTGTACTCTCTGTTCATcc                               <  1:3728132/100‑1 (MQ=255)
                                                                      cgcgcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATcc                               >  2:2486467/1‑100 (MQ=255)
                                                                      cgcgcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATcc                               <  2:836432/100‑1 (MQ=255)
                                                                        cgcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcag                             <  2:2521117/100‑1 (MQ=255)
                                                                        cgcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcag                             >  2:2764712/1‑100 (MQ=255)
                                                                          cgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagca                           >  1:959190/1‑100 (MQ=255)
                                                                          cgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagca                           >  2:1982974/1‑100 (MQ=255)
                                                                           gATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagcag                          >  2:1990293/1‑100 (MQ=255)
                                                                              tGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtg                       <  2:1479606/100‑1 (MQ=255)
                                                                               ggggCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtgt                      >  1:1524597/1‑100 (MQ=255)
                                                                                gggCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtgtg                     <  1:13888/100‑1 (MQ=255)
                                                                                 ggCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTgg                    >  2:1431408/1‑100 (MQ=255)
                                                                                   cAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGa                  >  1:513673/1‑100 (MQ=255)
                                                                                    aaaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGat                 >  1:1948070/1‑100 (MQ=255)
                                                                                    aaaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGat                 <  2:2999674/100‑1 (MQ=255)
                                                                                     cattccagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCCCCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCACCTGTTGAGGGata                >  1:3201738/7‑100 (MQ=21)
                                                                                     cattccagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGata                >  1:3201828/7‑100 (MQ=37)
                                                                                     aaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGata                >  1:60559/1‑100 (MQ=255)
                                                                                     aaaaaaaGCAGCGGT‑ACGTCGTTACCGATGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGata                >  2:2661740/1‑100 (MQ=255)
                                                                                       aaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGataat              >  1:3577677/1‑100 (MQ=255)
                                                                                         ccagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAatat            >  1:965965/3‑100 (MQ=38)
                                                                                         ccagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAatat            >  2:3431459/3‑100 (MQ=38)
                                                                                          cagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc           >  2:129232/2‑100 (MQ=38)
                                                                                           agcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg          >  2:3301775/1‑100 (MQ=38)
                                                                                           agcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg          >  2:433479/1‑100 (MQ=38)
                                                                                            gcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg         >  2:2839811/1‑100 (MQ=255)
                                                                                            gcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg         >  2:38890/1‑100 (MQ=255)
                                                                                            gcagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg         <  2:419721/100‑1 (MQ=255)
                                                                                             cagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGc        <  1:3839715/100‑1 (MQ=255)
                                                                                             cagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGc        >  1:1233811/1‑100 (MQ=255)
                                                                                               gcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAg      >  2:1834119/1‑100 (MQ=255)
                                                                                                cGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAg      <  2:2088808/99‑1 (MQ=255)
                                                                                                 ggT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGa    >  1:3404035/1‑100 (MQ=255)
                                                                                                   t‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGAtt  >  2:1345952/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATT  >  NC_000913/1977124‑1977322

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