Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,857 A→T 17.4% intergenic (+24/+22) 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,8570AT17.4% 380.1 / 68.4 154intergenic (+24/+22)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (67/60);  new base T (15/12);  total (82/72)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

Only 100 of 164 total aligned reads displayed.

CGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAC  >  NC_000913/3171758‑3171956
                                                                                                   |                                                                                                   
cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                     <  2:3379960/100‑1 (MQ=255)
cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                     >  2:1800973/1‑100 (MQ=255)
 gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                    <  1:534038/100‑1 (MQ=255)
 gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                    >  2:1449113/1‑100 (MQ=255)
    aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                 <  2:443809/100‑1 (MQ=255)
    aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                 <  1:1972436/100‑1 (MQ=255)
     aaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                 <  1:498102/99‑1 (MQ=255)
          catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGctct                                                                                           <  1:1878276/100‑1 (MQ=38)
           atggcatgAATGTTGAATTTGGGAATGCGGAACAAGCTGGATTTGAGGCGAAGGTAAGCTGGTATCTTCTCTGTCTATGAGAGCCGTTAAAACGACtctc                                                                                          <  1:1385795/100‑1 (MQ=255)
            tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                         >  2:3357296/1‑100 (MQ=255)
            tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                         <  1:2018128/100‑1 (MQ=255)
             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)
                 tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGac                                                                                    <  2:2294402/100‑2 (MQ=37)
                 tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGac                                                                                    <  2:2889255/100‑2 (MQ=37)
                       ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                              <  2:3530484/100‑1 (MQ=255)
                       ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                              <  2:1220015/100‑1 (MQ=255)
                       ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                              >  1:2749455/1‑100 (MQ=255)
                        tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                             >  1:866753/1‑100 (MQ=255)
                         gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                            >  1:88861/1‑100 (MQ=255)
                         gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                            >  1:1259122/1‑100 (MQ=255)
                         gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGCAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                            >  1:806306/1‑100 (MQ=255)
                           aTTTGGGAATTCGGAACAAGGGGGTTTGGAGGCGGAGTTAAGCTGGTAACTTCTCTGTCTATGAGATCCGTTAAAACTACTCTCATAGATTTTACTAATa                                                                          >  1:3183933/1‑100 (MQ=255)
                           aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                          >  2:2718794/1‑100 (MQ=255)
                           aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                          <  2:524706/100‑1 (MQ=255)
                           aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTGAAAACGACTCTCATAGATTTTACTAATa                                                                          >  2:1677168/1‑100 (MQ=255)
                            tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTATATGAGAGCCGTTAAAAAGACTCTCATAGATTTTACTAATa                                                                          >  1:802150/1‑99 (MQ=255)
                             ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                        <  1:336346/100‑1 (MQ=255)
                               gggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCa                                                                       >  1:1093152/1‑99 (MQ=255)
                                  aaTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                                   <  2:602704/100‑1 (MQ=255)
                                   aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                  >  2:783544/1‑100 (MQ=255)
                                   aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                  <  1:2246134/100‑1 (MQ=255)
                                   aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                  >  2:1402778/1‑100 (MQ=255)
                                     gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                <  2:3658949/100‑1 (MQ=255)
                                      cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                               >  1:1453044/1‑100 (MQ=255)
                                        gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                             <  1:3133966/100‑1 (MQ=255)
                                              ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                       >  1:564830/1‑100 (MQ=255)
                                                tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATGTAAACCGTCCCCa                                                     <  2:2590553/100‑1 (MQ=255)
                                                  gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaa                                                   >  2:23555/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)
                                                      tGAGGCGAAGGTAAACTGGTAACTTCTCTGTCTATGAGAGCAGTAAAAACGACTATCATAGCTCTTACTACTAGCAAAAAAGAAACAGCACCCaaaaaac                                               >  1:3539741/1‑99 (MQ=255)
                                                         ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                            >  1:2259522/1‑100 (MQ=255)
                                                         ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                            >  2:3717921/1‑100 (MQ=255)
                                                             aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaa                                        <  2:2852551/100‑1 (MQ=255)
                                                               ggTAAGCTGGTAACTTCTCTGTCTATGCGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                      <  2:1510251/100‑1 (MQ=255)
                                                               ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                      >  2:661306/1‑100 (MQ=255)
                                                               ggTAAGCTGGTAACTAGTCTGTCTATGAGAGCCGTTAAAACGACTCACATAGATTTTACTAATAGCAAAATATAAAAAAAAAAAAAAAAAGCCaccaacc                                      <  1:782290/100‑1 (MQ=255)
                                                                 tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                    <  2:826228/100‑1 (MQ=255)
                                                                  aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                   <  1:1115252/100‑1 (MQ=255)
                                                                  aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                   >  1:3187785/1‑100 (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:2729946/1‑100 (MQ=255)
                                                                         tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                            >  2:3786860/1‑100 (MQ=255)
                                                                         tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                            >  2:1422031/1‑100 (MQ=255)
                                                                          aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCTTCCCCAAAAAAGCCACCAACCACAAACCAGAc                           >  1:371069/1‑100 (MQ=255)
                                                                          aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                           <  2:2126094/100‑1 (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)
                                                                                   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)
                                                                                     cTATGAGAGTCGTATTAACGGATCGCATAGAATTTACTAATGGCAAAATACAAAACATCCCCAAAAAAACAAACAACCACAAACCAGACCAGTACcga                  >  1:706315/1‑98 (MQ=255)
                                                                                       aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga              >  2:516698/1‑100 (MQ=38)
                                                                                       aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga              >  2:2916841/1‑100 (MQ=38)
                                                                                       aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACAAGACAAATACcgaacga              >  1:246703/1‑100 (MQ=37)
                                                                                       aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga              >  2:410984/1‑100 (MQ=255)
                                                                                       aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga              >  1:2531626/1‑100 (MQ=255)
                                                                                       aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga              >  2:693625/1‑100 (MQ=255)
                                                                                        tGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag             >  2:1966230/1‑100 (MQ=38)
                                                                                        tGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTACCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag             >  1:1723494/1‑100 (MQ=38)
                                                                                        tGAGAGCCGTTAAAACGACTCTCATAGATTTTACGAACAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag             <  1:3494794/100‑1 (MQ=255)
                                                                                         gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCACCCACAAACCAGACAAATACCGAGCgaga            >  2:343037/1‑100 (MQ=255)
                                                                                         gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga            >  1:48370/1‑100 (MQ=255)
                                                                                         gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga            >  2:3680602/1‑100 (MQ=255)
                                                                                         gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga            >  2:3205796/1‑100 (MQ=255)
                                                                                         gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga            <  2:1749451/100‑1 (MQ=255)
                                                                                         cttccgatctAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga            <  1:1914228/91‑1 (MQ=255)
                                                                                           gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt          <  2:499519/100‑1 (MQ=255)
                                                                                           gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt          >  1:2633678/1‑100 (MQ=255)
                                                                                           gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCAACAACCACAAACCAGACAAATACCGAGCGAGAAt          >  1:2421826/1‑100 (MQ=255)
                                                                                             gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc        >  2:1450483/1‑100 (MQ=255)
                                                                                             gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc        <  2:333397/100‑1 (MQ=255)
                                                                                              ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa       <  2:3296432/100‑1 (MQ=255)
                                                                                               cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACat      <  1:1205408/100‑1 (MQ=255)
                                                                                                gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata     >  1:113547/1‑100 (MQ=255)
                                                                                                   aaaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaagatcggaag  >  2:1914237/1‑91 (MQ=255)
                                                                                                   |                                                                                                   
CGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAC  >  NC_000913/3171758‑3171956

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: