Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,858 A→T 15.5% intergenic (+25/+21) 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,8580AT15.5% 572.3 / 85.5 219intergenic (+25/+21)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (85/99);  new base T (15/19);  total (101/118)
Fisher's exact test for biased strand distribution p-value = 8.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 237 total aligned reads displayed.

GCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACC  >  NC_000913/3171759‑3171957
                                                                                                   |                                                                                                   
gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTGTGTCTATGAGAGCCGTTaa                                                                                                     >  2:3900365/1‑100 (MQ=255)
gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                     >  1:3785353/1‑100 (MQ=255)
 ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                                    >  2:1670432/1‑100 (MQ=255)
    aaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACg                                                                                                 <  1:1250061/100‑1 (MQ=255)
      ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAc                                                                                               <  1:826100/100‑1 (MQ=255)
       tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAct                                                                                              <  2:1074213/100‑1 (MQ=255)
         catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGctct                                                                                            <  2:4133574/100‑1 (MQ=38)
         catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActct                                                                                            <  2:2142551/100‑1 (MQ=255)
            ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                         <  1:1189426/100‑1 (MQ=38)
            ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                         <  2:3088322/100‑1 (MQ=38)
            ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCAt                                                                                         <  1:252828/100‑1 (MQ=255)
             gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                        <  1:62327/100‑1 (MQ=38)
             gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                        <  1:3838859/100‑1 (MQ=38)
             gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                        <  1:3114103/100‑1 (MQ=38)
             gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                        <  1:230933/100‑1 (MQ=38)
             gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                        <  2:1230569/100‑1 (MQ=38)
              catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAg                                                                                       >  1:1070402/1‑100 (MQ=255)
                tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGac                                                                                     <  1:1638566/100‑2 (MQ=37)
                 gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtt                                                                                    <  1:2641162/100‑1 (MQ=255)
                   aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGTGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtttt                                                                                  <  1:3992019/100‑1 (MQ=255)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                <  1:4127943/100‑1 (MQ=255)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                >  1:3227341/1‑100 (MQ=255)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                <  1:3985892/100‑1 (MQ=255)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                >  2:1579091/1‑100 (MQ=255)
                     gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                >  2:3119621/1‑100 (MQ=255)
                      ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                               <  1:3511515/100‑1 (MQ=255)
                       tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                              >  2:2544589/1‑100 (MQ=255)
                       tGAATGTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTTTCTGTCTAGGAGAGCCGTTAAAACGACTCTCATAGAtttt                                                                                  <  2:1433074/96‑1 (MQ=255)
                           tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                          >  1:1058802/1‑100 (MQ=255)
                           tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCCCATAGATTTTACTAATAg                                                                          >  2:754949/1‑100 (MQ=255)
                           nttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTATGTCTATGAGCTCCGTTAAAAATACTCTCATAGATTTTACTAATAg                                                                          >  1:1010708/2‑100 (MQ=255)
                               ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                      >  1:2253985/1‑100 (MQ=255)
                               ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                      >  2:3132998/1‑100 (MQ=255)
                                gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                                     >  1:3705844/1‑100 (MQ=255)
                                    gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                 <  1:3119634/100‑1 (MQ=255)
                                    gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                 >  2:230831/1‑100 (MQ=255)
                                    gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                 <  1:95575/100‑1 (MQ=255)
                                    gCGGAACAAGGTAGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                 >  1:2758693/1‑100 (MQ=255)
                                     cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                                >  1:3282350/1‑100 (MQ=255)
                                       gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                              <  1:916407/100‑1 (MQ=255)
                                       gAACAAGGGGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                              <  1:1839503/100‑1 (MQ=255)
                                           aaGGTGGATTTGAGGCGAAGATAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                          >  2:2561478/1‑100 (MQ=255)
                                            aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAACAGCAAAATATAAACCGTcc                                                         <  2:2877904/100‑1 (MQ=255)
                                             ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                        <  2:1536703/100‑1 (MQ=255)
                                              gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcccc                                                       <  1:2409201/100‑1 (MQ=255)
                                                ggATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACCCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaa                                                     >  1:1160923/1‑100 (MQ=255)
                                                     tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                >  1:4149205/1‑100 (MQ=255)
                                                     tGAGGCGAAGGTAAGCTGATAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGACCCCaaaaaa                                                 >  2:65488/1‑99 (MQ=255)
                                                       aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcc                                              <  1:1456947/100‑1 (MQ=255)
                                                       aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcc                                              <  1:260870/100‑1 (MQ=255)
                                                        ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                             >  2:3066614/1‑100 (MQ=255)
                                                         gcaaagctaagctGGTACATTACATGTTTATGCGAGCCGTTAAAACCACTCTAATAGATTTTACTAATATTAAAATATAAACCTTTCCCATATAAgtc                                              <  2:2440953/95‑3 (MQ=255)
                                                         gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccac                                            <  1:71581/100‑1 (MQ=255)
                                                          cGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                           >  1:561833/1‑100 (MQ=255)
                                                          cGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                           <  1:175803/100‑1 (MQ=255)
                                                                tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                     <  2:3495385/100‑1 (MQ=255)
                                                                 aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAAAAAAAAAAAAAAAGCCACCAACcaca                                    <  2:1499211/100‑1 (MQ=255)
                                                                  aGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaa                                   <  1:1532515/100‑1 (MQ=255)
                                                                   gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                  <  2:2729163/100‑1 (MQ=255)
                                                                    cTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGCCTCTCATAGATTTTCCTAATAGCAAAATATAAACCGCCCCCAAACAAGCCACCAACCCCCAAc                                 >  2:3562403/1‑100 (MQ=255)
                                                                    cTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                 <  1:1643441/100‑1 (MQ=255)
                                                                      ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                               <  2:2108091/100‑1 (MQ=255)
                                                                       gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                              <  2:1421142/100‑1 (MQ=255)
                                                                        tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                             <  1:2210935/100‑1 (MQ=255)
                                                                         aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAGCCAGAc                            <  2:1549329/100‑1 (MQ=255)
                                                                         aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                            <  2:379891/100‑1 (MQ=255)
                                                                          aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                           <  1:2853484/100‑1 (MQ=255)
                                                                           cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                          >  1:3323290/1‑100 (MQ=255)
                                                                            ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                         >  2:710223/1‑100 (MQ=255)
                                                                            ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                         <  2:3995047/100‑1 (MQ=255)
                                                                             tctcTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAt                        >  1:2895692/1‑100 (MQ=255)
                                                                              ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATa                       >  2:1358788/1‑100 (MQ=255)
                                                                                ctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                     <  1:148728/100‑1 (MQ=255)
                                                                                ctGTCTATGAGAGCCGTTAAAACACCTCTCATAGCTTTTACTCATAGCCACATATAAAGCGACCCCAAAACCGCCACCATCCGCAGTCCAGACAAATa                       <  2:1207505/98‑1 (MQ=255)
                                                                                 tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg                    >  1:1793687/1‑100 (MQ=255)
                                                                                  gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                   >  2:4160659/1‑100 (MQ=255)
                                                                                   tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                  >  1:1897973/1‑100 (MQ=38)
                                                                                    cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                 >  1:4061058/1‑100 (MQ=38)
                                                                                    cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                 >  1:2112490/1‑100 (MQ=38)
                                                                                    cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                 >  1:2448893/1‑100 (MQ=38)
                                                                                    cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACAAATAGCAAAATATAAACCGTCCTCAAAAAAGACACTAACCACCAACCAGACAAATCCcgaac                 >  2:3390623/1‑98 (MQ=21)
                                                                                    cTATGAGAGTCGTTTTAACGGCTCTCATAGATTGTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACACACCACAAAACAGACAAATAc                      >  2:2775196/1‑95 (MQ=25)
                                                                                    cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                 <  2:161598/100‑1 (MQ=255)
                                                                                     tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                   >  2:3612535/1‑97 (MQ=38)
                                                                                     tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                >  1:593188/1‑100 (MQ=255)
                                                                                     tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                >  2:81899/1‑100 (MQ=255)
                                                                                      aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga               >  2:2238/1‑100 (MQ=38)
                                                                                      aTGAGAGTCGTTTTAACGGCTATCATAGATTTTACTAATAGGAAAATATAAAGCGTCCCCAAAAGAGCGAGCCACCACAAACCAGACAAATACctag                  >  2:19302/1‑97 (MQ=21)
                                                                                      aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga               <  2:1052144/100‑1 (MQ=255)
                                                                                        gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgata             <  1:3884799/100‑3 (MQ=255)
                                                                                             ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa        <  2:1173295/100‑1 (MQ=255)
                                                                                             ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa        <  1:532446/100‑1 (MQ=255)
                                                                                              cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACat       >  2:3283029/1‑100 (MQ=255)
                                                                                               gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata      <  2:2390491/100‑1 (MQ=255)
                                                                                                   aaaCGCCTCTCATAGATTTTCCTAATACCAAAAAATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCGGACAAATACCGAGCGAGAATACATATAcc  <  2:2977438/100‑1 (MQ=255)
                                                                                                   aaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc  <  2:3228324/100‑1 (MQ=255)
                                                                                                   aaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc  <  2:1984235/100‑1 (MQ=255)
                                                                                                   aaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc  >  1:2052677/1‑100 (MQ=255)
                                                                                                   aaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc  <  1:3635698/100‑1 (MQ=255)
                                                                                                   aaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc  <  2:886929/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACC  >  NC_000913/3171759‑3171957

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: