Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,977,231 Δ1 bp 16.6% intergenic (‑92/+35) 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,2300T.16.6% 421.6 / 82.7 145intergenic (‑91/+36)motA/flhCproton conductor component of flagella motor/flagellar class II regulon transcriptional activator, with FlhD
Reads supporting (aligned to +/- strand):  ref base T (62/59);  new base . (16/8);  total (78/67)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 152 total aligned reads displayed.

TAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAA  >  NC_000913/1977133‑1977329
                                                                                                  |                                                                                                   
taagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTa                                                                                                   >  2:1354306/1‑100 (MQ=255)
    cACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCGAAAAAAAGCAGCGTT‑ACGTCGTTACCgc                                                                                               >  1:1033879/1‑100 (MQ=255)
    cACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgc                                                                                               <  2:3016533/100‑1 (MQ=255)
     aCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgct                                                                                              <  2:273193/100‑1 (MQ=255)
     aCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgct                                                                                              <  2:426543/100‑1 (MQ=255)
      cGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgctg                                                                                             >  1:2085104/1‑100 (MQ=255)
       gACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgc                                                                                            <  2:2250096/100‑1 (MQ=255)
        aCATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgct                                                                                           <  1:2841331/100‑1 (MQ=38)
        aCATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgct                                                                                           <  2:1020729/100‑1 (MQ=38)
         catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGctgctg                                                                                          <  2:1271918/100‑1 (MQ=255)
           tcatcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGa                                                                                        >  1:3535431/1‑100 (MQ=255)
           tcatcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGa                                                                                        >  2:2851920/1‑100 (MQ=255)
            catcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgcttttt                                                                                       <  2:2031751/100‑5 (MQ=37)
            catcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGaa                                                                                       >  1:2962217/1‑100 (MQ=255)
             atcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAAt                                                                                      <  2:2916714/100‑1 (MQ=255)
               ccTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGt                                                                                    >  1:1327543/1‑100 (MQ=255)
                cTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCGCTGCTAGatcgga                                                                                   >  1:1507681/1‑95 (MQ=37)
                cTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGtt                                                                                   >  2:3381453/1‑100 (MQ=255)
                cTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGtt                                                                                   >  2:758163/1‑100 (MQ=255)
                  tCCACTGTTGACAATGACAGGATCTTCAGTCGTCAGTCGTTAAAGAGCGATTGGATTAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTgc                                                                                 <  1:1592680/100‑1 (MQ=255)
                    cACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTgcgc                                                                               <  2:194396/100‑1 (MQ=255)
                     aCTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                              <  2:1149216/100‑1 (MQ=255)
                        gTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCa                                                                           >  2:3310897/1‑100 (MQ=255)
                         ttGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAc                                                                          <  1:2675729/100‑1 (MQ=255)
                         ttGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAc                                                                          <  2:766750/100‑1 (MQ=255)
                          tGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAcc                                                                         <  1:2357957/100‑1 (MQ=255)
                          tGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAcc                                                                         >  2:949362/1‑100 (MQ=255)
                           gACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg                                                                        <  2:1480404/100‑1 (MQ=255)
                           gACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg                                                                        >  2:2319779/1‑100 (MQ=255)
                           gACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg                                                                        >  2:3157309/1‑100 (MQ=255)
                            aCCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                                                       >  1:1168837/1‑100 (MQ=255)
                             ccATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTa                                                                      <  1:2561970/100‑1 (MQ=255)
                             ccATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTa                                                                      >  2:3447182/1‑100 (MQ=255)
                               aTGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATc                                                                    <  2:5039/100‑1 (MQ=255)
                                tGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa                                                                   >  1:3682238/1‑100 (MQ=255)
                                tGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa                                                                   <  1:1395056/100‑1 (MQ=255)
                                 gACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAg                                                                  >  2:2535246/1‑100 (MQ=255)
                                  aCAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGt                                                                 >  1:848846/1‑100 (MQ=255)
                                  aCAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGt                                                                 >  2:456767/1‑100 (MQ=255)
                                   cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                                >  1:1505945/1‑100 (MQ=255)
                                   cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                                >  1:1122250/1‑100 (MQ=255)
                                   cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                                >  1:3182033/1‑100 (MQ=255)
                                    aGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTa                                                               <  2:872930/100‑1 (MQ=255)
                                     ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                              >  1:762696/1‑100 (MQ=255)
                                     ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                              <  2:920467/100‑1 (MQ=255)
                                     ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                              >  2:800143/1‑100 (MQ=255)
                                      gATGTTCAGTCGTCAGGCGTTAACGCGAAAAAAAAAAAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaaa                                                             <  1:1248970/100‑1 (MQ=255)
                                       aTGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAAc                                                            <  2:1327551/100‑1 (MQ=255)
                                         gTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAg                                                          >  2:3139971/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)
                                           tCAGTCGTCAGGCGTTAACGCGCGATGGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                        <  1:1480074/100‑1 (MQ=255)
                                             aGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTg                                                      <  2:3312017/100‑1 (MQ=255)
                                               tcgtcAGGCGTTCACTCGTAAATGTGATCAAAAAAAGCAGCGTT‑ACGTAGTTAGCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                                       <  1:3445029/97‑1 (MQ=255)
                                                cgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAc                                                   >  1:3720183/1‑100 (MQ=255)
                                                 gtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAct                                                  >  2:1465262/1‑100 (MQ=255)
                                                 gtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTACCAGTTAAACAGCCTGTAct                                                  >  1:1012545/1‑100 (MQ=255)
                                                    aGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctc                                               <  2:2869101/100‑1 (MQ=255)
                                                     ggCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctct                                              >  2:1688319/1‑100 (MQ=255)
                                                      gCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTg                                             <  2:349802/100‑1 (MQ=255)
                                                       cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGt                                            <  2:1487460/100‑1 (MQ=255)
                                                       cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCGGt                                            >  2:1535177/1‑100 (MQ=255)
                                                        gTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                           <  2:1166948/100‑1 (MQ=255)
                                                          tAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCa                                         <  2:2141573/100‑1 (MQ=255)
                                                           aaCGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCAt                                        <  2:3411267/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)
                                                                  gATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagcag                                 >  2:1990293/1‑100 (MQ=255)
                                                                       gggCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtgtg                            >  1:3437286/1‑100 (MQ=255)
                                                                       gggCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtgtg                            <  1:983902/100‑1 (MQ=255)
                                                                        ggCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTgg                           >  1:527148/1‑100 (MQ=255)
                                                                           aaaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGat                        <  2:2999674/100‑1 (MQ=255)
                                                                           aaaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGat                        >  1:1948070/1‑100 (MQ=255)
                                                                            cattccagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCCCCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCACCTGTTGAGGGata                       >  1:3201738/7‑100 (MQ=21)
                                                                            cattccagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGata                       >  1:3201828/7‑100 (MQ=37)
                                                                            aaaaaaaGCAGCGGT‑ACGTCGTTACCGATGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGata                       >  2:2661740/1‑100 (MQ=255)
                                                                             attccagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGataa                      >  1:3799627/6‑100 (MQ=37)
                                                                              aaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGataat                     >  1:3577677/1‑100 (MQ=255)
                                                                                ccagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAatat                   >  1:2901011/3‑100 (MQ=38)
                                                                                 cagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc                  >  2:2332317/2‑100 (MQ=38)
                                                                                  agcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg                 >  2:433479/1‑100 (MQ=38)
                                                                                  agcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg                 >  2:2254507/1‑100 (MQ=38)
                                                                                  agcagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg                 <  2:1258496/100‑1 (MQ=255)
                                                                                   gcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg                >  1:1980310/1‑100 (MQ=255)
                                                                                   gcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg                >  2:38890/1‑100 (MQ=255)
                                                                                   gcagcGGT‑AGGTCGTTACCGCTGCTGGAATGTTTCGGTTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg                <  1:1717418/100‑1 (MQ=255)
                                                                                    cagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGc               <  1:3839715/100‑1 (MQ=255)
                                                                                    cagcGGA‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTAATCCAGCAGTTATGGGATAATATCGGc               >  1:3122082/1‑100 (MQ=255)
                                                                                     agcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCa              <  2:3694825/100‑1 (MQ=255)
                                                                                      gcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAg             >  2:1834119/1‑100 (MQ=255)
                                                                                       cGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAg             <  2:2088808/99‑1 (MQ=255)
                                                                                        ggT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGa           >  2:2776699/1‑100 (MQ=255)
                                                                                        ggT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGa           >  1:3404035/1‑100 (MQ=255)
                                                                                             cgtcgtTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCt       >  2:669578/1‑100 (MQ=255)
                                                                                             cgtcgtTACCGCTGCTGGAATGGTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCTAGCAGTTGCGGGATAATATCGGCAGGATTCt       <  1:2249021/100‑1 (MQ=255)
                                                                                              gtcgtTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTg      >  2:1420987/1‑100 (MQ=255)
                                                                                               tcgtTACCGCTGCTGGAATGTTGCGCCTCACCGTAGCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTgg     <  1:955740/100‑1 (MQ=255)
                                                                                                  ttACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGaa  <  1:3310871/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAA  >  NC_000913/1977133‑1977329

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: ‑