Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,857 A→T 14.5% intergenic (+24/+22) qseC → / ← ygiZ sensory histidine kinase QseC/conserved inner membrane protein YgiZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,8570AT14.5% 189.4 / 22.4 76intergenic (+24/+22)qseC/ygiZsensory histidine kinase QseC/conserved inner membrane protein YgiZ
Reads supporting (aligned to +/- strand):  ref base A (38/27);  new base T (7/4);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

GGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAAT  >  NC_000913/3171708‑3172004
                                                                                                                                                     |                                                                                                                                                   
ggACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTTGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                                                                     <  2:1416056/150‑1 (MQ=255)
ggACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                                                                     <  1:230059/150‑1 (MQ=255)
    aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                                                                 <  2:1375051/150‑1 (MQ=255)
     aaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACg                                                                                                                                                >  1:178988/1‑150 (MQ=255)
       ccGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAc                                                                                                                                              >  1:1004396/1‑150 (MQ=255)
         gCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGTAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                            <  1:324928/150‑1 (MQ=255)
             ccGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCAt                                                                                                                                        >  1:970139/1‑150 (MQ=255)
               ggCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTTGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                                                                      <  2:861363/150‑1 (MQ=255)
               ggCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                                                                      <  2:1014789/150‑1 (MQ=255)
               ggCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAg                                                                                                                                      >  2:485309/1‑150 (MQ=255)
                  aGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtt                                                                                                                                   >  2:643584/1‑150 (MQ=255)
                   gCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                                                                      >  2:59750/1‑146 (MQ=255)
                   gCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                                                                      <  1:59750/146‑1 (MQ=255)
                    cGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtttt                                                                                                                                 <  1:35113/150‑1 (MQ=255)
                      ggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                                                               <  2:427764/150‑1 (MQ=255)
                      ggCTTGGGCTATCGATTCTCCATCTAATCTCAAAATTCCATGGCATGAATGATTAATTTGTCAATGCGCAACAAGGGGGATTTGATGCGAAGGTAACATGCTAACTTCTCTGTCTATGAGTGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                                                               <  2:233384/150‑1 (MQ=255)
                         ttGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                                                                            <  2:290130/150‑1 (MQ=255)
                           gggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                          <  1:1115767/150‑1 (MQ=255)
                            ggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                                                                         >  2:598245/1‑150 (MQ=255)
                                   gATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                                                                  >  1:1340165/1‑150 (MQ=255)
                                       gTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                                                              >  1:114527/1‑150 (MQ=255)
                                         ccAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAcc                                                                                                            >  1:789521/1‑150 (MQ=255)
                                             cGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcc                                                                                                        <  1:658202/150‑1 (MQ=255)
                                               aaTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcccc                                                                                                      >  1:320510/1‑150 (MQ=255)
                                               aaTCGCCAAATTGCATGGCATGAATCTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcccc                                                                                                      <  2:356946/150‑1 (MQ=255)
                                                    ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                                                                                 >  1:387488/1‑150 (MQ=255)
                                                    ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                                                                                 >  1:689542/1‑150 (MQ=255)
                                                         ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                                                            <  2:178988/150‑1 (MQ=255)
                                                             atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaa                                                                                        <  2:1420258/150‑1 (MQ=255)
                                                             atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaa                                                                                        >  1:797852/1‑150 (MQ=255)
                                                                gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacca                                                                                     <  1:1253038/150‑1 (MQ=255)
                                                                 catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTATATTAGAGCAGTTAAAACGACTCTCATAGATTTTACTATTAGCAAAATATAAACCGTGCCCAAAAAAGGCACCAACcac                                                                                    <  2:536869/150‑1 (MQ=255)
                                                                   tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                                                            >  1:717530/1‑140 (MQ=255)
                                                                   tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                                                            <  2:717530/140‑1 (MQ=255)
                                                                    gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAGCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                                                                 <  2:1074050/150‑1 (MQ=255)
                                                                        gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                                                                             >  1:466769/1‑150 (MQ=255)
                                                                          tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                                                                           >  1:487438/1‑150 (MQ=255)
                                                                            aaTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                                                                         >  2:232082/1‑150 (MQ=255)
                                                                                 gggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                                                                    >  1:892163/1‑150 (MQ=255)
                                                                                  ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg                                                                   >  1:32292/1‑150 (MQ=255)
                                                                                   gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                                                                  <  1:232082/150‑1 (MQ=255)
                                                                                         ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTACCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                                                            >  1:1409828/1‑150 (MQ=255)
                                                                                            aCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa                                                         >  1:748807/1‑150 (MQ=255)
                                                                                             cAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc                                                        <  2:242407/150‑1 (MQ=255)
                                                                                                  tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGAGAGAATACAtata                                                   >  1:184883/1‑150 (MQ=255)
                                                                                                    gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc                                                 >  1:443525/1‑150 (MQ=255)
                                                                                                           ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                          >  1:137244/1‑150 (MQ=255)
                                                                                                           ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACAGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                          >  1:359653/1‑150 (MQ=255)
                                                                                                             cGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCg                                        >  1:928408/1‑150 (MQ=255)
                                                                                                               aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                      <  2:970139/150‑1 (MQ=255)
                                                                                                               aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                      <  2:487438/150‑1 (MQ=255)
                                                                                                                    aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACAttt                                 >  2:218615/1‑150 (MQ=255)
                                                                                                                     aGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTc                                >  1:632140/1‑150 (MQ=255)
                                                                                                                      gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAGCCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTcc                               >  2:51549/1‑150 (MQ=255)
                                                                                                                      gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCCCCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTcc                               <  2:1004396/150‑1 (MQ=255)
                                                                                                                      gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTcc                               <  2:1340165/150‑1 (MQ=255)
                                                                                                                        tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCg                             <  1:643584/150‑1 (MQ=255)
                                                                                                                          gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCa                           >  2:148990/1‑150 (MQ=255)
                                                                                                                          gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCa                           <  2:32292/150‑1 (MQ=255)
                                                                                                                            aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATa                         >  1:733410/1‑150 (MQ=255)
                                                                                                                             aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGCCAAATACCGAGCGAGAATACATATACCCAAAACTCGAGCAGTTCCCGCATaa                        >  2:979931/1‑150 (MQ=255)
                                                                                                                             aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaa                        <  2:744438/150‑1 (MQ=255)
                                                                                                                              cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                       >  2:217631/1‑150 (MQ=255)
                                                                                                                              cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                       >  2:1292953/1‑150 (MQ=255)
                                                                                                                                 taaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAagag                    >  1:189755/6‑150 (MQ=39)
                                                                                                                                        tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt             >  1:1089373/1‑150 (MQ=255)
                                                                                                                                        tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCCCCAACCACAAACCAGACAAATACCGAGCGAGAATACAAATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt             >  2:738251/1‑150 (MQ=255)
                                                                                                                                         aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTta            <  2:320510/150‑1 (MQ=255)
                                                                                                                                         aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTta            <  2:689542/150‑1 (MQ=255)
                                                                                                                                          tGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAACCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa           >  1:119876/1‑150 (MQ=255)
                                                                                                                                           gagaGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGTATACATATACCCaa                                              <  2:1252676/114‑1 (MQ=255)
                                                                                                                                           gagaGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGTATACATATACCCaa                                              >  1:1252676/1‑114 (MQ=255)
                                                                                                                                           gagaGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATAAACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaag          >  2:1073177/1‑150 (MQ=255)
                                                                                                                                           gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCGTAAAGAGTTTTCTtaag          >  1:271654/1‑150 (MQ=255)
                                                                                                                                            agagTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACACAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaaga         >  2:889840/1‑150 (MQ=255)
                                                                                                                                                 cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataaga    >  1:487927/1‑150 (MQ=255)
                                                                                                                                                 cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataaga    >  2:514010/1‑150 (MQ=255)
                                                                                                                                                   ttAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaat  <  2:466769/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAAT  >  NC_000913/3171708‑3172004

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