Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,491,722 T→A 9.0% intergenic (‑102/+3) yhfG ← / ← ppiA putative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,491,7220TA 9.0% 2266.2 / 26.9 756intergenic (‑102/+3)yhfG/ppiAputative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base T (236/423);  new base A (90/2);  total (331/425)
Fisher's exact test for biased strand distribution p-value = 1.36e-33
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 765 total aligned reads displayed.

CTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGAAATC  >  NC_000913/3491573‑3491871
                                                                                                                                                     |                                                                                                                                                     
   aaGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                                                    <  2:12924277/150‑1 (MQ=255)
        cccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGAATGCCCGCGCGAGAAATAAATACgg                                                                                                                                               >  1:6247164/1‑150 (MQ=255)
         ccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGAATGCCCGCGCGCGAAATCATTAAGGc                                                                                                                                              >  1:13022650/1‑150 (MQ=255)
         ccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTTCCCGCGCGAGAAAACATTAAGGc                                                                                                                                              >  1:9491674/1‑150 (MQ=255)
          cAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCa                                                                                                                                             <  2:14920365/150‑1 (MQ=255)
               aCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACt                                                                                                                                        >  2:10432455/1‑150 (MQ=255)
               aCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTGCCCGCGCGAGAAATCATTACGGCAGGACt                                                                                                                                        >  2:4885920/1‑150 (MQ=255)
                  ggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTGCCCCCGCGAGAAAACAAAAAGGCAGGACattt                                                                                                                                     >  2:13707573/1‑149 (MQ=255)
                   gACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAg                                                                                                                                    <  2:6067299/150‑1 (MQ=255)
                    aCTTTTGCTTATCGTTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGc                                                                                                                                   <  1:5691912/150‑1 (MQ=255)
                    aCTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGc                                                                                                                                   <  2:12555366/150‑1 (MQ=255)
                    aCTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGc                                                                                                                                   <  2:14756136/150‑1 (MQ=255)
                       tttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGa                                                                                                                                <  1:11238718/150‑1 (MQ=255)
                        ttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATGAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATATTGCCCGCGCGAGTAATAATTAACGCATGACATTATGGGaa                                                                                                                               >  1:2617631/1‑150 (MQ=255)
                          gCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAg                                                                                                                             >  2:2325928/1‑150 (MQ=255)
                           cTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAgg                                                                                                                            <  2:3077046/150‑1 (MQ=255)
                           cTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAgg                                                                                                                            <  1:3316862/150‑1 (MQ=255)
                            ttATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGa                                                                                                                           <  2:7245368/150‑1 (MQ=255)
                            ttATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGa                                                                                                                           <  2:7826155/150‑1 (MQ=255)
                             tATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTGCCCCCCCGCGGAAAAAATACAGCAGGACATTATCGGAAAGGAt                                                                                                                          >  1:14996037/1‑150 (MQ=255)
                              aTCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATa                                                                                                                         <  1:9877489/150‑1 (MQ=255)
                                     aGTTTCTTCACTTGATGCCTCCCTGACCGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCg                                                                                                                  <  1:6020337/150‑1 (MQ=255)
                                       tttCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTTCCCGCGCGAGAAATCAATAAAGCAGGGCATTTTCGGCAAGGATAACAAACggg                                                                                                                >  1:2749154/1‑149 (MQ=255)
                                       tttCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGAATTCCCCCGCCGCAAAACAATAACGCAGGGCGTTAGCGGCAAGGATAACAAACggg                                                                                                                >  2:14491359/1‑149 (MQ=255)
                                         tcttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCATATCGGGGTTACGCGGGGAGAGCTGAAGCGGGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGttt                                                                                                              <  1:5480298/150‑1 (MQ=255)
                                           ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGGCCGAGTGGCCGCGCGGGAACACCAAACGGCCCGACGTTAGCGGGAAGGATAAAAACCCGGTTTg                                                                                                            >  1:6757969/1‑150 (MQ=255)
                                           ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTTCCCGCGCGTGAAATCAATAAGGCAGGAGATTAGCGGGAAGGATAACAAAAGGttttt                                                                                                            >  2:5615121/1‑149 (MQ=255)
                                            tcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCAAGAAAGCCCCCCCGGGCAAAAAATAAAGCAGGACGATAGCGGAAAGGATAAAAACCCGGTTTgt                                                                                                           >  1:15391103/1‑149 (MQ=255)
                                            tcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCCCGAGAAAACATTAAGGCAGGACGTTTTCGGGAAGGATAAAAACCGGttttta                                                                                                           >  1:9186276/1‑148 (MQ=255)
                                             cACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCAGATTGCCCGCGCGCGAAATCAGTACGGCCGGGCGGTAGCGGGCAGGATAACTACCCGGGTTGAc                                                                                                          >  2:10900850/1‑150 (MQ=255)
                                             cACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGAc                                                                                                          <  2:5404375/150‑1 (MQ=255)
                                              aCTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCAGATTGGCCGCGCGGGAAATAAATAAGGCAGGACAGTTTAGGGAAGGATAACGAACGGgtttttcg                                                                                                         >  1:3817260/1‑146 (MQ=255)
                                               cTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACgg                                                                                                        <  2:11390477/150‑1 (MQ=255)
                                                ttGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGc                                                                                                       >  2:2602301/1‑150 (MQ=255)
                                                       cTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGGCAGCGTGGCCCCCCGCGAAAAAAATACGGCCGGACGTTGGGGGGAACGGTAACAAACGGGTTTTTCGGcccaatat                                                                                                >  2:9487266/1‑146 (MQ=255)
                                                         cccTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTgg                                                                                              <  1:12619343/150‑1 (MQ=255)
                                                          ccTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCACAATGCCCCCGCGCGCAAACAATACCGCCGGACTTTAGCCGAAAGGATAACAACCGGGTTTGACGGCACATTCTGGt                                                                                             >  2:5952278/1‑150 (MQ=255)
                                                          ccTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTTCCCCCCCGAGCAATCAATAACGCCCGACTTTAGGGGGAAGGATAACTACCGGGTTTTTCGGCACAATATGGt                                                                                             >  2:11096665/1‑150 (MQ=255)
                                                            tGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGAATGCCCCCCCGCGAAAACAGTACGGCCGGCCATTATCGGGAAAGATAACAACCCGGTTTGTCGGCCCAATCTGGTAc                                                                                           >  2:11965656/1‑150 (MQ=255)
                                                             gAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGGCCGCGCGAGAAAAAAATACGGCAGGACGTTAGCGGGAAGGAAAACAACCCGTTTTTAAGGCACAATATgggaag                                                                                          >  1:10361961/1‑146 (MQ=255)
                                                             gAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACg                                                                                          <  1:1002487/150‑1 (MQ=255)
                                                              aaCGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACgg                                                                                         <  1:14373854/150‑1 (MQ=255)
                                                               aCGTGTCCATTTTCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCTGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGa                                                                                        <  1:502771/150‑1 (MQ=255)
                                                               aCGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGa                                                                                        <  2:333762/150‑1 (MQ=255)
                                                                    tCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAAc                                                                                   <  2:641771/150‑1 (MQ=255)
                                                                    tCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAAc                                                                                   <  2:8149026/150‑1 (MQ=255)
                                                                           gCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCCGATTGCCCCCCCCCGAAAACAATAACGCACGCCTGTAGAGGAAACGCAAACAACCCCGTTTGAAGGCACCATAAGTTACCGACCAAccacaaga                                                                            >  2:471726/1‑143 (MQ=255)
                                                                           gCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGa                                                                            <  2:9286167/150‑1 (MQ=255)
                                                                           gCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGa                                                                            <  1:13537749/150‑1 (MQ=255)
                                                                            ccACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACGGTAGCGGAAAGGATAACTACCGGGTTTGACGGCACATTCTGGTACGGACCCACGTCATGAg                                                                           >  2:8980067/1‑150 (MQ=255)
                                                                              acaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTTACGGCCCAATATGGTACGGACCAACCTCAAGAGTc                                                                         >  2:8774774/1‑150 (MQ=255)
                                                                              acaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTTCCCCCGCGAGGAAACATTAAGGCAGGACGTTAGCGGAAAGGATAACTACCGGGTTTTACGGCACAATCTGGTACGGACCAACGGCATGAGTc                                                                         >  2:12296471/1‑150 (MQ=255)
                                                                                aaGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTTCCCGCGCGAGAAATAATTACGGCCGGACTTTAGCGGAAAGGATAACAACCGGTTTTGACGGCACATTATGGTACGGACCAACGTCATGAGTCgg                                                                       >  2:3718266/1‑150 (MQ=255)
                                                                                   tataAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCAc                                                                    <  1:7078954/150‑1 (MQ=255)
                                                                                     taAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCt                                                                  <  2:10217667/150‑1 (MQ=255)
                                                                                        gCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCAACAATGGCCGCGCGGGGAATAATTACCGCAGGACATTAGCGGCAAGGGTGACAACCCGGTTTTACGGCACAATCTGGTAAGGACCAACCGCCTGAGTCCGCCCCAggg                                                               >  2:10538236/1‑150 (MQ=255)
                                                                                         cGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGACAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGa                                                              >  2:14162879/1‑150 (MQ=255)
                                                                                         cGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGACAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGa                                                              >  2:1657476/1‑150 (MQ=255)
                                                                                          ggCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGaa                                                             <  1:13824503/150‑1 (MQ=255)
                                                                                            cAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCAGATTGGCCCCCCGGGAAAACAAAAAGGCAGGGCATTATCGGAAAAGGTAACAAACCGTTTTTACGGCACAATAATGTACGGACCCACCTCATGAGTCCGCCCCGGGGaaaa                                                           >  2:6939924/1‑149 (MQ=255)
                                                                                             aaaTCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATc                                                          <  1:4020858/150‑1 (MQ=255)
                                                                                                  cGGGTTACGCCGGGAGAGCAGAAGCGCGGGCGAGATGGCCCGCGCGGGAGAACATTACGGCAGGGCTTTAGCGGCAAGGCTAACGTACGGTTTTTCCGGCACATTCGGGTACGGACCCACGTCCAGAGTCGGCACCTGGGACATCTtctc                                                     >  2:3616257/1‑150 (MQ=255)
                                                                                                  cGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTTCCCGCCCGCGAAATAAATACGGCCGGACTTTAGCGGGAAGGATAACTAACGGGTTTTACGGCACATTCTGGTACGGACCAACGTCAAGAGTCGGCACCGGGGCAATCTTATc                                                     >  2:13510797/1‑150 (MQ=255)
                                                                                                   gggTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTTCCCGCGCGAGAAATAATTAACGCAGGACTTTAGCGGAAAGGATAACAACCCGTTTTTAAGGCACATTATGGTAAGGCCCAACCTCATGAGTCGGCACCGGGGGAATCTTATCg                                                    >  2:7677347/1‑150 (MQ=255)
                                                                                                          gCCGGGAGAGCAGAAGCGCGGGCAAGATTTCCCCCGCGAGAAATCATTAAGGCAGGACATTAGCGGGAAGGGTAACTACCGGGTTTTACGGCCCAATATGGTACCGACCAACGTCATGAGTCGGCACCAGGGGAATCTTATCGGCCAcct                                             >  2:9722869/1‑148 (MQ=255)
                                                                                                           ccGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTc                                            <  1:9643765/150‑1 (MQ=255)
                                                                                                             gggAGAGCAGATGCGCGGGCATGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGAACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCa                                          <  1:2044199/150‑1 (MQ=255)
                                                                                                              ggAGAGCAGAAGCGCGGGCCAGATGGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCAt                                         <  2:187557/150‑1 (MQ=255)
                                                                                                              ggAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCTGTCCTTTAGCGGAAAGGATAACTACCGTTTTTGACGGCACATTCTGGTACGGACCAACGTAATGAGTCGGCACCTGGTAAATCTTATCGGCAACGTCCAt                                         <  1:8125360/150‑1 (MQ=255)
                                                                                                                agagCAGATGCGCGGGAAAGTTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAATGGATAACTACCGGTTTTGACGGCACTTTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGc                                       <  1:11965656/150‑1 (MQ=255)
                                                                                                                agagCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGc                                       <  2:279851/150‑1 (MQ=255)
                                                                                                                agagCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGc                                       <  1:1097987/150‑1 (MQ=255)
                                                                                                                agagCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGc                                       <  2:13147837/150‑1 (MQ=255)
                                                                                                                  agCAGAAGCGCGGGCAAGATTTCCCGCGCGAGAAAACAATACGGCAGGACTTTAGCGGAAAGGATAAATAACGGTTTTTACGGCACAATCTGGTAAGGACCAACCACATGAGAAGGCACCTGGGAAAACATATCGGCAAACTCCAAGCCt                                     >  2:15181210/1‑150 (MQ=255)
                                                                                                                   gCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAACATTACGGCCGGACTTTAGCGGAAAGGATAACTACCGGTTTTGAAGGCACAATCTGGTACGGACCAACGTCATGAGTCCGCACCTGGGAAATCTTATCGGCAACGTCCCTGCCtt                                    >  2:8068108/1‑150 (MQ=255)
                                                                                                                    cAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCttt                                   <  2:4355605/150‑1 (MQ=255)
                                                                                                                     aGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc                                  <  1:5953079/150‑1 (MQ=255)
                                                                                                                       aaGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcac                                <  1:13652490/150‑1 (MQ=255)
                                                                                                                            cgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACtt                           >  1:11395059/1‑150 (MQ=255)
                                                                                                                              ggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa                         <  1:8774774/150‑1 (MQ=255)
                                                                                                                              ggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa                         >  2:15457529/1‑150 (MQ=255)
                                                                                                                              ggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa                         >  1:12215003/1‑150 (MQ=255)
                                                                                                                                cAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTTGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAcc                       <  1:43266/150‑1 (MQ=255)
                                                                                                                                cAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAcc                       <  1:654441/150‑1 (MQ=255)
                                                                                                                                 aaGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGACATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTATACCa                      >  2:11680316/1‑150 (MQ=255)
                                                                                                                                   gagtGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaaa                    <  1:10673728/147‑1 (MQ=255)
                                                                                                                                   gATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaaa                    <  1:7456460/150‑1 (MQ=255)
                                                                                                                                   gATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaaa                    <  1:1664182/150‑1 (MQ=255)
                                                                                                                                       gCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATAcc                <  1:8407649/150‑1 (MQ=255)
                                                                                                                                         ccGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgc              <  2:8531603/150‑1 (MQ=255)
                                                                                                                                         ccGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgc              <  2:13631793/150‑1 (MQ=255)
                                                                                                                                          cgcgcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgcg             >  1:3681103/1‑150 (MQ=255)
                                                                                                                                          cgcgcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgcg             >  1:771055/1‑150 (MQ=255)
                                                                                                                                           gcgcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGt            <  1:2298318/150‑1 (MQ=255)
                                                                                                                                            cgcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTa           >  2:12695913/1‑150 (MQ=255)
                                                                                                                                             gcgAGAAATCATTACGGCATGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTaa          <  1:471726/150‑1 (MQ=255)
                                                                                                                                             gcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTaa          <  2:214250/150‑1 (MQ=255)
                                                                                                                                              cgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAAc         >  2:5809823/1‑150 (MQ=255)
                                                                                                                                                 gaAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGa      <  1:14812489/150‑1 (MQ=255)
                                                                                                                                                  aaaTCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGaa     >  2:14023141/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGAAATC  >  NC_000913/3491573‑3491871

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑