Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,977,227 T→A 16.6% intergenic (‑88/+39) 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,2270TA16.6% 358.4 / 66.0 145intergenic (‑88/+39)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 A (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 = 4.26e-02

Only 100 of 151 total aligned reads displayed.

AGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGG  >  NC_000913/1977130‑1977326
                                                                                                  |                                                                                                   
agataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgtcg                                                                                                   >  1:282388/1‑100 (MQ=255)
       cACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCGAAAAAAAGCAGCGTT‑ACGTCGTTACCgc                                                                                            >  1:1033879/1‑100 (MQ=255)
        aCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgct                                                                                           <  2:426543/100‑1 (MQ=255)
        aCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgct                                                                                           <  2:273193/100‑1 (MQ=255)
         cGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgctg                                                                                          >  1:2085104/1‑100 (MQ=255)
           tccgatctcTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgct                                                                                        <  2:1507690/92‑1 (MQ=37)
           aCATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgct                                                                                        <  2:1020729/100‑1 (MQ=38)
           aCATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgct                                                                                        <  1:2841331/100‑1 (MQ=38)
           aCATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgctgct                                                                                        <  2:2535310/100‑1 (MQ=255)
            catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgctt                                                                                       <  2:12172/100‑2 (MQ=38)
            catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGctgctg                                                                                       <  2:1271918/100‑1 (MQ=255)
            catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGctgctg                                                                                       <  1:2952879/100‑1 (MQ=255)
               catcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgcttttt                                                                                    <  2:2031751/100‑5 (MQ=37)
               catcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGaa                                                                                    >  1:2962217/1‑100 (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:2166000/100‑1 (MQ=255)
                       cACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTgcgc                                                                            <  2:194396/100‑1 (MQ=255)
                           gTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCa                                                                        >  2:2744651/1‑100 (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                                                                      >  2:949362/1‑100 (MQ=255)
                             tGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAcc                                                                      >  1:3411239/1‑100 (MQ=255)
                             tGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAcc                                                                      <  1:2357957/100‑1 (MQ=255)
                              gACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg                                                                     >  2:3157309/1‑100 (MQ=255)
                               aCCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                                                    >  1:1168837/1‑100 (MQ=255)
                               aCCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                                                    >  1:3592223/1‑100 (MQ=255)
                                ccATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTa                                                                   >  2:3447182/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)
                                    gACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAg                                                               <  2:2853645/100‑1 (MQ=255)
                                    gACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAg                                                               >  2:2535246/1‑100 (MQ=255)
                                    gACAGGATGTCCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAg                                                               <  2:2085757/100‑1 (MQ=255)
                                    gACAGGATGGTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAg                                                               <  1:911132/100‑1 (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:3133319/100‑1 (MQ=255)
                                      cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                             >  1:3182033/1‑100 (MQ=255)
                                      cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                             >  1:1122250/1‑100 (MQ=255)
                                       aGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTa                                                            <  2:872930/100‑1 (MQ=255)
                                        ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                           <  2:920467/100‑1 (MQ=255)
                                        ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                           <  2:1782085/100‑1 (MQ=255)
                                        ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                           >  2:800143/1‑100 (MQ=255)
                                         gATGTTCAGTCGTCAGGCGTTAACGCGAAAAAAAAAAAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaaa                                                          <  1:1248970/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)
                                              tCAGTCGTCAGGCGTTAACGCGCGATGGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                     <  1:1480074/100‑1 (MQ=255)
                                                aGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTg                                                   <  2:3312017/100‑1 (MQ=255)
                                                 gtcgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGt                                                  >  1:3454297/1‑100 (MQ=255)
                                                  tcgtcAGGCGTTCACTCGTAAATGTGATCAAAAAAAGCAGCGTT‑ACGTAGTTAGCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                                    <  1:3445029/97‑1 (MQ=255)
                                                   cgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAc                                                >  2:1046300/1‑100 (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)
                                                           gTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                        <  2:1166948/100‑1 (MQ=255)
                                                             tAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCa                                      <  2:2141573/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)
                                                                    cgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagca                               >  2:1982974/1‑100 (MQ=255)
                                                                    cgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagca                               >  1:959190/1‑100 (MQ=255)
                                                                        tGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtg                           <  2:1479606/100‑1 (MQ=255)
                                                                           ggCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTgg                        >  2:1431408/1‑100 (MQ=255)
                                                                             cAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGa                      <  2:2944661/100‑1 (MQ=255)
                                                                             cAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGa                      >  1:513673/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‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGata                    >  1:60559/1‑100 (MQ=255)
                                                                                 aaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGataat                  >  1:3577677/1‑100 (MQ=255)
                                                                                   ccagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAatat                >  1:965965/3‑100 (MQ=38)
                                                                                    cagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc               >  2:2332317/2‑100 (MQ=38)
                                                                                    cagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc               >  2:129232/2‑100 (MQ=38)
                                                                                     agcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg              >  2:433479/1‑100 (MQ=38)
                                                                                     agcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg              >  2:2254507/1‑100 (MQ=38)
                                                                                      gcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg             >  2:38890/1‑100 (MQ=255)
                                                                                      gcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg             >  1:1980310/1‑100 (MQ=255)
                                                                                      gcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg             >  2:2839811/1‑100 (MQ=255)
                                                                                      gcagcGGT‑AGGTCGTTACCGCTGCTGGAATGTTTCGGTTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg             <  1:1717418/100‑1 (MQ=255)
                                                                                       cagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGc            <  1:3839715/100‑1 (MQ=255)
                                                                                        agcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCa           <  2:3694825/100‑1 (MQ=255)
                                                                                           ggT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGa        >  2:2776699/1‑100 (MQ=255)
                                                                                           ggT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGa        >  1:3404035/1‑100 (MQ=255)
                                                                                             t‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGAtt      >  2:1345952/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)
                                                                                                  |                                                                                                   
AGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGG  >  NC_000913/1977130‑1977326

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