Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,858 A→T 6.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,8580AT 6.5% 201.9 / 4.7 77intergenic (+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 (36/35);  new base T (2/3);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

AAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATA  >  NC_000913/3171712‑3172005
                                                                                                                                                  |                                                                                                                                                   
aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTTACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                   >  1:218531/1‑149 (MQ=255)
aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                   <  1:438310/149‑1 (MQ=255)
   ccGCTACCGGCAGCGGGCGTGGGCTATCGATTGTCCAGCGACTCGCCAAATTGCAGGGCATGAATGTTGAATTTGGTCTTGCGGAACAACGGGGAGTTGAGGCGAAGTTAAGCTGGTAACTTCTCTGTCTATGAGAGCAGTTAACACgc                                                                                                                                                >  1:676168/1‑148 (MQ=255)
    cGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAATGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAc                                                                                                                                               <  1:181242/149‑1 (MQ=255)
    cGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGAGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTTGATGTGAGGCGAAGGTAAGCACGTAACTTCTCTGTCTATGAGAGCCGTTAAATCgcc                                                                                                                                               >  2:698107/1‑147 (MQ=255)
      gtACCGGCTGCGGGGTTGGGCTTTGGATTGTCCTGCGAATCGCCATATTGGTTGGCTTGAATGTTGAATTTGGGAATGCGGAACATGGTGGATTTGTGGCGAAGGTATGCTGGTAAATTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                             <  2:437921/148‑1 (MQ=255)
       tACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGATTCGCCATATTGCATGGCATGTATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGctct                                                                                                                                            <  2:523103/149‑1 (MQ=255)
          cGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                                                                         <  1:272708/149‑1 (MQ=255)
          cGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                                                                         <  2:430697/149‑1 (MQ=255)
             cAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGATTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGa                                                                                                                                      >  1:481352/1‑149 (MQ=255)
                cGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCTGGGCAGGAATGTTGAATTTGGGACTGCGGACCAAGGTGGATTTGCGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTACAACGACTCTCATAGAttt                                                                                                                                   <  1:722683/149‑1 (MQ=255)
                  ggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATTTTGAATTTGGGAATGCTGAACAAGGTGGATTTTAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                                                                 <  2:10565/149‑1 (MQ=255)
                   gCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                                                                <  2:17317/149‑1 (MQ=255)
                   gCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGAGTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                                                                <  1:558810/149‑1 (MQ=255)
                    cTTGGGCTATCGATTGTCCAGCGATTCGCCAAATTGCATGCCATGAATGTTGAATTTGTGAATGCGGAACATGGTGTATTTGAGGCGTAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                                                                               <  2:576071/149‑1 (MQ=255)
                       gggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAt                                                                                                                            <  1:238886/149‑1 (MQ=255)
                        ggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGGGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                           <  1:203202/149‑1 (MQ=255)
                        ggCTATCGACTGTCCAGCGAATCCCCAAATTGCATGGCATGAATGTTGAACTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                           <  2:340194/149‑1 (MQ=255)
                                    tCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                                                               <  1:750739/149‑1 (MQ=255)
                                         cGAATAGCCAAATTTCATGGCATGAATGTTTAATTTGGGAATGCGGAACAAGTTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                                                                          <  1:64886/149‑1 (MQ=255)
                                            aTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcccc                                                                                                       <  2:525201/149‑1 (MQ=255)
                                             tCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcacaa                                                                                                      >  2:261880/1‑145 (MQ=255)
                                             tCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGAGTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCa                                                                                                      >  1:558351/1‑149 (MQ=255)
                                                  aactctcatggcctgAATGTTTAATTTGCGAATCCGTAACAAGGTGGCTTTGAGGCGAATGTCAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                 <  2:951360/143‑1 (MQ=255)
                                                    aTTGCATGGCTTGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTGTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                                                               <  2:418125/149‑1 (MQ=255)
                                                    aTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                                                               <  2:771914/149‑1 (MQ=255)
                                                        catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATATGAGGCGAAGGTAAGCTGGAAACATCTCTGTCTATGAGAGACGTTAAAACGACTCTCATATATTTTAATAATAGCAAAATAAAAAAATACCCCAAAAAACccacc                                                                                           >  2:585195/1‑149 (MQ=255)
                                                         atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAACCGACTCACATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAACAGacacca                                                                                          >  2:572253/1‑149 (MQ=255)
                                                         atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGacacca                                                                                          >  1:917733/1‑149 (MQ=255)
                                                           ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACACTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                                                                        >  1:362711/1‑149 (MQ=255)
                                                              atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                                                                     >  1:287343/1‑149 (MQ=255)
                                                                     ttGAATTTGGGAATGCGGAACACGGTGGATTTGAGGCGATGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCACAAAAAAGCCAccaaccaccaaccac                                                                              >  2:773531/1‑148 (MQ=255)
                                                                      tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                                                             <  1:438356/149‑1 (MQ=255)
                                                                      tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTAAAAACGACTCTCATAGATTTTACTAATATCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAgc                                                                             >  2:216856/1‑148 (MQ=255)
                                                                      tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGGCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                                                             <  1:770230/149‑1 (MQ=255)
                                                                       gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGTTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                                                                            <  1:89081/149‑1 (MQ=255)
                                                                       gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCCTAGATTTTACTAATAGCAATATATAAGCCGTCCCCAAAAAAGCCGCCAACCACAAACCAGAc                                                                            >  2:469734/1‑149 (MQ=255)
                                                                        aaTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCAAAAAAAGCCACCAACCACAAACCAGACa                                                                           >  1:575920/1‑149 (MQ=255)
                                                                         aTTTGGGAATTCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATAAAAACCGTCCCCATAAAAGCCCCCAACCACAAACCAGACaa                                                                          >  2:482221/1‑149 (MQ=255)
                                                                               gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTATCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg                                                                    <  2:833949/149‑1 (MQ=255)
                                                                                 aTGCTGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                  <  2:212092/149‑1 (MQ=255)
                                                                                   gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                                                                >  1:761801/1‑149 (MQ=255)
                                                                                    cGGAACAAGGTGGATTTGAGTCGAATGTAAGCTGGTAACTTCTCTGTCTATGATAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCTCCAACCACAAACCAGACAAATACcgagcga                                                               <  1:874527/149‑1 (MQ=255)
                                                                                      gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                                                             >  1:458923/1‑149 (MQ=255)
                                                                                            ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACat                                                       <  1:19426/149‑1 (MQ=255)
                                                                                            ggTGGATTGGAGGTGAAGGTTATCTTGTAATTTCTCGTGATTTGAGAGCCGTTAAAAAGACTTTCATAGTTTTTACTAATATACATATATAAAACGTCCCCATCAACGCCACCAACCACAAACCAGACAAATACCGCGCGAGAATAcct                                                       <  2:676168/149‑3 (MQ=255)
                                                                                             gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACACACCAGACAAATACCGAGCGAGAATACata                                                      >  2:826525/1‑149 (MQ=255)
                                                                                                gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAc                                                   <  2:761483/149‑1 (MQ=255)
                                                                                                 aTTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACTACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc                                                  <  1:504916/149‑1 (MQ=255)
                                                                                                 aTTTGAGGCGAAGGTAACCTGCAAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCACAAAGATTTTACTAAAAGCAAAATATAACCCGTCCCCAAAAAAGCCACCAACCCCAAACCCGCCAAATGCCGAGCGCGAATACAGATccc                                                  >  2:699389/1‑149 (MQ=255)
                                                                                                       ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAAc                                            >  1:350900/1‑149 (MQ=255)
                                                                                                        gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                           >  2:477877/1‑149 (MQ=255)
                                                                                                        gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCACATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCATCCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                           >  1:91777/1‑149 (MQ=255)
                                                                                                          gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTATTAGCAAAATATAAACCGTCCCCAAAAAAGCCACAAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCg                                         >  1:34132/1‑149 (MQ=255)
                                                                                                          gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCg                                         <  1:897927/149‑1 (MQ=255)
                                                                                                             ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAACCCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAAc                                      >  1:642668/1‑149 (MQ=255)
                                                                                                             ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAAc                                      <  2:860657/149‑1 (MQ=255)
                                                                                                              gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAGACCGTCCCCAAAAAAGCCACCAACCACAAACCACACAAATACCGCGCGAGAATACATATAACCAAAACTCGAACa                                     >  2:921050/1‑149 (MQ=255)
                                                                                                                aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGAGAGAATACATATACCCAAAACTCGAACAtt                                   >  1:126578/1‑149 (MQ=255)
                                                                                                                        aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCAt                           >  1:279280/1‑149 (MQ=255)
                                                                                                                           ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAATCCGTCCCCAAACAAGCCACCAACCACACACCAGCCAAATACCGCGCGAGAATACATCTACCCAAAACTCGAACAGTTCCCGCATaaa                        >  2:790230/1‑149 (MQ=255)
                                                                                                                           ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCCCCATCCACAAAGCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACAGTTCCCGCATaag                        >  1:54201/1‑148 (MQ=255)
                                                                                                                           ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                        >  2:551152/1‑149 (MQ=255)
                                                                                                                             ctctGTCTATGATATCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAACGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAaga                      <  2:143725/149‑1 (MQ=255)
                                                                                                                              tctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAagag                     <  1:976602/149‑1 (MQ=255)
                                                                                                                                  tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                 <  2:939540/149‑1 (MQ=255)
                                                                                                                                     aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt              <  2:390091/149‑1 (MQ=255)
                                                                                                                                      tGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGCACATTTCCCGCATAAAGAGTTTTCTta             >  1:469263/1‑149 (MQ=255)
                                                                                                                                      tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAACAGACAAAAACCGAGAGAGAATACATATACCCAAAACTCGAACATTAAACGAATAAAGAGTTTTCTta             >  2:790707/1‑149 (MQ=255)
                                                                                                                                       gagaGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa            >  1:138497/1‑149 (MQ=255)
                                                                                                                                        agagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAAATATACCCAAAACTCGAACATTACCCGCATAAAGCGTTTTCTtaag           >  1:169240/1‑149 (MQ=255)
                                                                                                                                        agagCCGGTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCACACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaag           <  2:558351/149‑1 (MQ=255)
                                                                                                                                         gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGGAAAATATAAACCGTCCCCAAAAAAGCCACCACCCACAAACCAGACAAATACCGAGCGAGAATACATATCACCAAAACTCGAACATTTCCCGCATAAAGCGTTTTCTtaagt          >  2:672554/1‑148 (MQ=255)
                                                                                                                                         gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAACCTCGAACATTTCCCGCATAAAGAGTTTTCTtaaga          <  2:176993/149‑1 (MQ=255)
                                                                                                                                          agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagat         <  2:75160/149‑1 (MQ=255)
                                                                                                                                             cGTTAAAACGACTCTCATAGATTTTACTAAAAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataag      <  2:215477/149‑1 (MQ=255)
                                                                                                                                              gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataaga     >  1:453812/1‑149 (MQ=255)
                                                                                                                                               ttAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTAAAGATAAGaa    >  2:457337/1‑149 (MQ=255)
                                                                                                                                               ttAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGATTTCTTAAGATAAGaa    >  2:20103/1‑149 (MQ=255)
                                                                                                                                                 aaaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaata  >  2:32628/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATA  >  NC_000913/3171712‑3172005

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: