Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,838 T→A 7.8% intergenic (‑32/‑255) ygbA ← / → mutS protein YgbA/DNA mismatch repair protein MutS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,8380TA 7.8% 258.7 / 9.4 89intergenic (‑32/‑255)ygbA/mutSprotein YgbA/DNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base T (40/42);  new base A (2/5);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACAAACAACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTC  >  NC_000913/2856689‑2856987
                                                                                                                                                     |                                                                                                                                                     
aaacaaacAACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTatat                                                                                                                                                       >  2:589432/1‑150 (MQ=255)
 aacaaacaACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATAtt                                                                                                                                                      >  1:818563/1‑150 (MQ=255)
     aacaacGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACat                                                                                                                                                  >  2:1513401/1‑150 (MQ=255)
        aacGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatctt                                                                                                                                               <  1:48355/150‑1 (MQ=255)
            ccTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACTCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatc                                                                                                                                           <  2:56609/150‑1 (MQ=255)
            ccTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatc                                                                                                                                           >  2:134357/1‑150 (MQ=255)
            ccTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatc                                                                                                                                           <  2:1239430/150‑1 (MQ=255)
            ccTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTAATATAAATACatcttaatg                                                                                                                                           <  1:569121/150‑2 (MQ=255)
             cTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatct                                                                                                                                          >  1:696102/1‑150 (MQ=255)
              tCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTAATATAAATACatcttaatgag                                                                                                                                         <  2:1457841/150‑4 (MQ=255)
              tCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTAATATAAATACatcttaatgag                                                                                                                                         <  1:501536/150‑4 (MQ=255)
               cGTAATGCTCCGGCTCCGGTGAGGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTATTCGTCAGTTTTTCACGAGAGTTACGCTTGGCGGACATGCTGGCTTCACCTCATTAAGATGTAATATAATTACatcttaatgagg                                                                                                                                        <  2:92841/150‑5 (MQ=255)
               cGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTAATATAAATACatcttaatgagg                                                                                                                                        <  2:215564/150‑5 (MQ=255)
               cGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTAATATAAATACatcttaatgagg                                                                                                                                        <  2:1251437/150‑5 (MQ=255)
                gTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACAtcttaatcttaa                                                                                                                                       >  2:143342/1‑150 (MQ=255)
                gTAATGCTCCGGCTCCGCTGACGCCTCCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACAtcttaatcttaa                                                                                                                                       <  1:1403904/150‑1 (MQ=255)
                  aaTGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAg                                                                                                                                     <  2:1317465/150‑1 (MQ=255)
                    tGCTCCGGCTCCGCTGCCGCCTGCGGGCATTTCGCGGGATTAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGGGAGATGCGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATTTTCCATCTTAATCTTAAAggg                                                                                                                                   >  1:1265139/1‑150 (MQ=255)
                    tGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatctta                                                                                                                                        >  1:1342774/1‑145 (MQ=255)
                    tGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatctta                                                                                                                                        <  2:1342774/145‑1 (MQ=255)
                        ccggctacgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTAATATTACATCTTAATCTTAAAGGGCACt                                                                                                                               >  2:714468/1‑150 (MQ=255)
                         cggctccgCTGACGCCTGCGGGCATTTCGCTTGATACAGATCGATCATTTTTTTAATCGTCCGTTTTTCACGCGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTa                                                                                                                              <  2:949482/150‑1 (MQ=255)
                         cggctccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTa                                                                                                                              <  1:1298736/150‑1 (MQ=255)
                           gctccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATg                                                                                                                            <  2:296576/150‑1 (MQ=255)
                           gctccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATg                                                                                                                            <  1:1513401/150‑1 (MQ=255)
                           gctccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATg                                                                                                                            <  2:702266/150‑1 (MQ=255)
                             tccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTTACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGAc                                                                                                                          <  1:1114454/150‑1 (MQ=255)
                             tccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGAc                                                                                                                          <  1:524980/150‑1 (MQ=255)
                              ccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACt                                                                                                                         <  1:457520/150‑1 (MQ=255)
                                 cTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCa                                                                                                                      <  1:1170535/150‑1 (MQ=255)
                                 cTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCa                                                                                                                      >  2:1369267/1‑150 (MQ=255)
                                       ccTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaagaag                                                                                                                <  2:1176274/150‑1 (MQ=255)
                                       ccTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaagaag                                                                                                                >  1:332438/1‑150 (MQ=255)
                                           cGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAgg                                                                                                            <  2:757400/150‑1 (MQ=255)
                                              gCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGtt                                                                                                         <  2:818563/150‑1 (MQ=255)
                                               cATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaaga                                                                                                                  >  1:228820/1‑140 (MQ=255)
                                               cATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaaga                                                                                                                  <  2:228820/140‑1 (MQ=255)
                                               cATTTCGCTTGATAAAGATCGATCATTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAg                                                                                                       <  2:640701/150‑1 (MQ=255)
                                                    cGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCaa                                                                                                   >  2:930551/1‑150 (MQ=255)
                                                     gCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaaga                                                                                                                  >  2:381628/1‑134 (MQ=255)
                                                     gCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaaga                                                                                                                  <  1:381628/134‑1 (MQ=255)
                                                     gCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAAc                                                                                                  <  1:47322/150‑1 (MQ=255)
                                                     gCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAAc                                                                                                  <  2:814993/150‑1 (MQ=255)
                                                        tGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGa                                                                                               >  1:326224/1‑150 (MQ=255)
                                                         gATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGAt                                                                                              >  2:800861/1‑150 (MQ=255)
                                                           tAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGCCTCCAAAGAAGAAGGGTTCGCCAACCGATAc                                                                                            >  2:1336101/1‑150 (MQ=255)
                                                           tAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATAc                                                                                            <  1:589432/150‑1 (MQ=255)
                                                               catcgatcATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAAtt                                                                                        <  2:740446/149‑1 (MQ=255)
                                                                  cgatcATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTg                                                                                     <  1:143342/150‑1 (MQ=255)
                                                                  cgatcATTTTCTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTg                                                                                     <  1:247969/150‑1 (MQ=255)
                                                                       aTTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTAc                                                                                >  1:896449/1‑150 (MQ=255)
                                                                       aTTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTAc                                                                                >  1:830394/1‑150 (MQ=255)
                                                                        tttttttAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACt                                                                               >  1:1106282/1‑150 (MQ=255)
                                                                           ttttAATCGTCAGTTTTTCACGAGAGATACGCTTGCCTGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGc                                                                            <  2:818006/150‑1 (MQ=255)
                                                                               aaTCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCa                                                                        <  1:930551/150‑1 (MQ=255)
                                                                                  cGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                                     >  2:1405349/1‑150 (MQ=255)
                                                                                   gTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                                    >  1:268405/1‑150 (MQ=255)
                                                                                    tCAGTTTTTCACGATAGATACTCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGTGTTAGACAACCGATACAATTTTGCGTACTTGCTTCATAAGc                                                                   <  2:1319827/150‑1 (MQ=255)
                                                                                    tCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGc                                                                   >  1:1011052/1‑150 (MQ=255)
                                                                                    tCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGc                                                                   >  2:904059/1‑150 (MQ=255)
                                                                                                gagagaTACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAg                                                       <  2:765910/150‑1 (MQ=255)
                                                                                                  gagaTACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCt                                                     >  2:206112/1‑150 (MQ=255)
                                                                                                            tGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGc                                           >  2:22621/1‑150 (MQ=255)
                                                                                                              ccGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCAt                                         <  1:1162245/150‑1 (MQ=255)
                                                                                                                     tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                                  >  1:1423507/1‑150 (MQ=255)
                                                                                                                     tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                                  <  2:326224/150‑1 (MQ=255)
                                                                                                                      gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                                 >  2:310192/1‑150 (MQ=255)
                                                                                                                      gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                                 >  2:713439/1‑150 (MQ=255)
                                                                                                                      gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                                 <  2:576004/150‑1 (MQ=255)
                                                                                                                       ctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCgg                                <  2:1386023/150‑1 (MQ=255)
                                                                                                                       ctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCgg                                <  2:332438/150‑1 (MQ=255)
                                                                                                                         gcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGaa                              >  2:104293/1‑150 (MQ=255)
                                                                                                                            tCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCa                           >  2:935567/1‑150 (MQ=255)
                                                                                                                            tCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCa                           >  1:521992/1‑150 (MQ=255)
                                                                                                                             ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAg                          >  2:1317746/1‑150 (MQ=255)
                                                                                                                             ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAg                          >  2:655702/1‑150 (MQ=255)
                                                                                                                                ccTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCAt                       >  2:438876/1‑150 (MQ=255)
                                                                                                                                 cTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATc                      <  1:1317746/150‑1 (MQ=255)
                                                                                                                                  tCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                     <  1:66257/150‑1 (MQ=255)
                                                                                                                                   cATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCaa                    >  2:1507649/1‑150 (MQ=255)
                                                                                                                                   cATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAg                   >  1:1114537/1‑150 (MQ=255)
                                                                                                                                    aTTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAg                   >  2:1038949/1‑150 (MQ=255)
                                                                                                                                    aTTAAGATGTAATATAAATACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAg                   >  2:1283818/1‑150 (MQ=255)
                                                                                                                                     ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa                  <  1:655702/150‑1 (MQ=255)
                                                                                                                                     ttAAGATGTATTTATATTACATCTTAATATTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa                  <  1:1434573/150‑1 (MQ=255)
                                                                                                                                     ttAAGATGTAATATAAATACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa                  >  2:1261204/1‑150 (MQ=255)
                                                                                                                                      tAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGaa                 >  2:364098/1‑150 (MQ=255)
                                                                                                                                         gATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTc              >  1:813821/1‑150 (MQ=255)
                                                                                                                                         gATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTc              >  2:1195890/1‑150 (MQ=255)
                                                                                                                                            gTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGcc           <  1:77814/150‑1 (MQ=255)
                                                                                                                                             tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCg          >  2:1495374/1‑150 (MQ=255)
                                                                                                                                                  atatTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGc     <  1:206112/150‑1 (MQ=255)
                                                                                                                                                    atTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGctt   >  1:865358/1‑150 (MQ=255)
                                                                                                                                                    atTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGctt   <  1:1405349/150‑1 (MQ=255)
                                                                                                                                                    atTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGctt   >  2:8685/1‑150 (MQ=255)
                                                                                                                                                     ttACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATATTTCCCGGAACCAGCATCAAGAACTCCCCGTTCGctta  >  2:102206/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACAAACAACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTC  >  NC_000913/2856689‑2856987

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: