Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,852 C→T 17.2% intergenic (+19/+27) qseC → / ← ygiZ quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,8520CT17.2% 387.2 / 60.0 163intergenic (+19/+27)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (71/64);  new base T (15/13);  total (86/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 165 total aligned reads displayed.

AATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACA  >  NC_000913/3171755‑3171951
                                                                                                 |                                                                                                   
aaTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCg                                                                                                   >  1:1055697/1‑100 (MQ=255)
aaTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCg                                                                                                   <  2:2475790/100‑1 (MQ=255)
  tCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGtt                                                                                                 <  2:2527886/100‑1 (MQ=255)
   cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                <  2:3379960/100‑1 (MQ=255)
   cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                >  2:1800973/1‑100 (MQ=255)
    gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                               >  2:1115259/1‑100 (MQ=255)
    gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                               >  1:495016/1‑100 (MQ=255)
    gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                               >  2:1449113/1‑100 (MQ=255)
    gCCAAATTGCATGGCATGAATGTGGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                <  1:3748303/99‑1 (MQ=255)
      cAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                             >  2:3446703/1‑100 (MQ=255)
       aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                            <  1:1212727/100‑1 (MQ=255)
       aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                            <  2:443809/100‑1 (MQ=255)
       aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                            <  2:1983550/100‑1 (MQ=255)
          ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAc                                                                                         <  1:2166243/100‑1 (MQ=255)
             catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGctct                                                                                      <  1:1878276/100‑1 (MQ=38)
                ggcatgAATGTTGTATTGGGAAATGCGGAACAAGGTGGTTTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                   <  1:3026647/100‑1 (MQ=37)
                ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                   <  2:1348397/100‑1 (MQ=38)
                ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                   <  1:1511015/100‑1 (MQ=38)
                ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                   <  1:2460118/100‑1 (MQ=38)
                 gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                  <  2:479444/100‑1 (MQ=38)
                  catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                 <  2:772190/100‑1 (MQ=38)
                  catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                 <  2:30878/100‑1 (MQ=38)
                  catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAg                                                                                 >  1:772798/1‑100 (MQ=255)
                    tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGac                                                                               <  2:2889255/100‑2 (MQ=37)
                         gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                          <  1:198017/100‑1 (MQ=255)
                          ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                         <  2:3530484/100‑1 (MQ=255)
                          ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                         >  1:2749455/1‑100 (MQ=255)
                           tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                        >  1:866753/1‑100 (MQ=255)
                            aaaTTTGAGAATGCGGAACAAGATGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                        >  1:2476698/2‑99 (MQ=255)
                              aTTTGGGAATTCGGAACAAGGGGGTTTGGAGGCGGAGTTAAGCTGGTAACTTCTCTGTCTATGAGATCCGTTAAAACTACTCTCATAGATTTTACTAATa                                                                     >  1:3183933/1‑100 (MQ=255)
                              aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                     <  2:524706/100‑1 (MQ=255)
                              aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACGCTCATAGATTTTACTAATa                                                                     >  1:1538786/1‑100 (MQ=255)
                              aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTGAAAACGACTCTCATAGATTTTACTAATa                                                                     >  2:1677168/1‑100 (MQ=255)
                                ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                   <  1:336346/100‑1 (MQ=255)
                                   ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                >  2:2526673/1‑100 (MQ=255)
                                     aaTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                              <  2:602704/100‑1 (MQ=255)
                                     aaTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                              >  2:2758829/1‑100 (MQ=255)
                                     aaTGCGGAACAACGTGGATTTGAGTCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                               <  1:3682378/99‑1 (MQ=255)
                                      aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                             >  2:783544/1‑100 (MQ=255)
                                      aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                             >  2:1402778/1‑100 (MQ=255)
                                      aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                             <  1:2246134/100‑1 (MQ=255)
                                       tGCGGAACAATGTGGATTTGAGGCGAAGGTAAGCTTGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCTTAGATTTTACTAATAGCAAAatat                                                            <  1:2425276/100‑1 (MQ=255)
                                        gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                           <  2:3658949/100‑1 (MQ=255)
                                           gcaccaGGTGGATTTGAGGCGAAGGTAAGCTGGTCACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCGCATACATTTTACTAATATCAAAATATAAAc                                                        <  1:702004/95‑1 (MQ=255)
                                           gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                        <  1:3133966/100‑1 (MQ=255)
                                             aCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCg                                                      >  2:1040023/1‑100 (MQ=255)
                                                 ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                  <  2:2981177/100‑1 (MQ=255)
                                                  gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcccc                                                 >  2:2164479/1‑100 (MQ=255)
                                                   tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATGTAAACCGTCCCCa                                                <  2:2590553/100‑1 (MQ=255)
                                                     gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaa                                              >  2:2351273/1‑100 (MQ=255)
                                                      aTTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaa                                             >  1:1761081/1‑100 (MQ=255)
                                                         tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCGGTTAAAACGACACACATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                          >  1:528613/1‑100 (MQ=255)
                                                         tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                          <  1:270491/100‑1 (MQ=255)
                                                            ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                       >  2:3717921/1‑100 (MQ=255)
                                                             gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccac                                      <  2:3688748/100‑1 (MQ=255)
                                                                aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaa                                   <  2:2852551/100‑1 (MQ=255)
                                                                aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaa                                   >  2:3354339/1‑100 (MQ=255)
                                                                  ggTAAGCTGGTAACTTCTCTGTCTATGCGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                 <  2:1510251/100‑1 (MQ=255)
                                                                  ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                 >  2:661306/1‑100 (MQ=255)
                                                                  ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                 <  1:1096467/100‑1 (MQ=255)
                                                                  ggTAAGCTGGTAACTAGTCTGTCTATGAGAGCCGTTAAAACGACTCACATAGATTTTACTAATAGCAAAATATAAAAAAAAAAAAAAAAAGCCaccaacc                                 <  1:782290/100‑1 (MQ=255)
                                                                    tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                               <  2:826228/100‑1 (MQ=255)
                                                                     aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                              >  1:3187785/1‑100 (MQ=255)
                                                                     aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                              <  1:1115252/100‑1 (MQ=255)
                                                                     aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACCCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                 <  1:3427823/97‑1 (MQ=255)
                                                                       gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                            >  1:556773/1‑100 (MQ=255)
                                                                        cTGGTAACTTCTCTGTCTATAAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                           <  1:1449105/100‑1 (MQ=255)
                                                                         tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                            >  1:688646/1‑98 (MQ=255)
                                                                         tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                          <  2:1087515/100‑1 (MQ=255)
                                                                          ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                         <  2:1259129/100‑1 (MQ=255)
                                                                          ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                         >  1:2396225/1‑100 (MQ=255)
                                                                          ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                         >  1:2723420/1‑100 (MQ=255)
                                                                            tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                       >  2:1422031/1‑100 (MQ=255)
                                                                             aaCTTCCCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                      <  2:733116/100‑1 (MQ=255)
                                                                              aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                     <  1:2758815/100‑1 (MQ=255)
                                                                               cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATCGATTTTACTAATAGCAAAATATAAACAGTCCCCAAAAAAGCAACGAACCACAAACCAGACaa                    >  1:1061096/1‑100 (MQ=255)
                                                                               cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                    <  1:3788438/100‑1 (MQ=255)
                                                                                   aaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                >  2:1883692/5‑100 (MQ=37)
                                                                                   aaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGCCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                >  1:803443/5‑100 (MQ=25)
                                                                                      gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga             <  2:1849575/100‑1 (MQ=255)
                                                                                      gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCAAAAAAAGCAACAAACCACAAACCAGACAAACACaga             >  1:129950/1‑100 (MQ=255)
                                                                                      atCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga             >  2:2297272/2‑100 (MQ=37)
                                                                                       tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATGAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga             >  1:2579758/1‑99 (MQ=37)
                                                                                        cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc           >  2:3184877/1‑100 (MQ=38)
                                                                                        cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAACACACCAACCACAAACCAGACAAGTAAcga             >  1:826599/1‑98 (MQ=25)
                                                                                        cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc           <  2:337389/100‑1 (MQ=255)
                                                                                        cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc           >  2:3217476/1‑100 (MQ=255)
                                                                                         tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg          >  1:1374741/1‑100 (MQ=255)
                                                                                          aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga         >  2:516698/1‑100 (MQ=38)
                                                                                          aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga         >  2:693625/1‑100 (MQ=255)
                                                                                          aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga         >  2:410984/1‑100 (MQ=255)
                                                                                           tGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag        >  2:1966230/1‑100 (MQ=38)
                                                                                           tGAGAGCCGTTAAAACGACTCTCATAGATTTTACGAACAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag        <  1:3494794/100‑1 (MQ=255)
                                                                                            gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga       <  2:1749451/100‑1 (MQ=255)
                                                                                             agagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag        <  1:1878527/98‑1 (MQ=255)
                                                                                              gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt     <  2:499519/100‑1 (MQ=255)
                                                                                              gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt     >  1:2633678/1‑100 (MQ=255)
                                                                                               agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAACGTCCCCAAAAAAGCCACCAACCACAAACAAGACAAATACCGAGCGAGAATa    >  1:5288/1‑100 (MQ=255)
                                                                                                gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc   >  2:1450483/1‑100 (MQ=255)
                                                                                                 ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa  <  2:3296432/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
AATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACA  >  NC_000913/3171755‑3171951

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