Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,485,026 T→C 12.8% E309G (GAG→GGG)  pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,485,0260TC12.8% 123.1 / 2.8 58E309G (GAG→GGG) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base T (22/27);  new base C (4/4);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

AGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCC  >  NC_000913/4484789‑4485258
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           
aGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGCCGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCCGTTTGAGCt                                                                                                                                                                                                                                             >  1:1309048/1‑238 (MQ=255)
  cgcgATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGAGAACGTTAATCGGCAGTTTTAGCTcc                                                                                                                                                                                                                                           >  1:953386/1‑238 (MQ=255)
      aTCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGCCAGTTGAAGCCCCGCGa                                                                                                                                                                                                                                       >  1:884369/1‑238 (MQ=255)
                          cGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTACGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCACGCATGTTTTCTCAGCCTGCCAACACGCCGATTACGTTAATCGTCAGTTGAACTCCCGCGGCCATCCGCATCACGCCGTaa                                                                                                                                                                                                                   >  1:1107662/1‑238 (MQ=255)
                           gTCGCCACTCCAATCCCCGCTCCCGCCTAAAAACGCTCAACGTAAGTTAACACGTCGCACAGTCCCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACCCGTCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAAc                                                                                                                                                                                                                 <  2:830860/239‑1 (MQ=255)
                            tCGCCACGTCAATCCCCGCTTCCGGCCAAAAACGCAAAACGTAAGTTAACCCGGCGCACAGTCCCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGCCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAAc                                                                                                                                                                                                                 <  1:609420/238‑1 (MQ=255)
                               ccACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTTTAGCTCCGCGACCATCCGCCTCACCCCGTAAACCGc                                                                                                                                                                                                              >  1:657399/1‑238 (MQ=255)
                               ccACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCACGCATGTTTTCACAGCCTGCCAACACGCCGTACCAGTTACTCGGCAGTTGTAGCCCCCCGCCCAACCCCCCCACCCCGCaaaacccc                                                                                                                                                                                                             >  2:679535/1‑234 (MQ=255)
                                   gTCAATCCCCGCTTCCGGCTCAAACCGCTCAACGGAAGTTAACCCGTCGCACAGTCCCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGccg                                                                                                                                                                                                        <  1:166036/239‑1 (MQ=255)
                                    ccaacacccccTTCCGGCCAAAAAGCCTCAACGGAAGTTAACCGCCACCATAATACCGGGCGGCTTTCAGGGCCGGTGTCAGCGCCTTAAAGGGTTTGACCGAACAGGGTGGTTAACACCGCCCCGGGACGATAGGCTCGTCCGCCAGGCATTTTTTCGCAGCCTGCAAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGccgc                                                                                                                                                                                                       <  2:1143684/229‑1 (MQ=255)
                                      aaTCCCCGCTCCCGTCTAAAACCCCCCAACGTAAGTTAACACGTCCACCAGTCCCAGGGGGCCTTACGGGTGGGTGTTCGGCACTTCAACGGTTTGCCCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCAAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGccgcac                                                                                                                                                                                                     <  2:413641/239‑1 (MQ=255)
                                         ccccctttCCGGCCCAAACCCCTCAAGGAAATTAAACCGGGCCGACGGTCCCAGGCGGCCTTCGGCGTCGGTGTTAGGCACTACAGCGGTTTGACCGGACGTGGGGGTTAACAAATCGCCTGGAGGATAGGCTCGTCCGCCAGGCATGTTTTCGCCCCCGTCCACCACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGccgcaccgc                                                                                                                                                                                                  <  2:1235598/233‑1 (MQ=255)
                                             gCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCCGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACat                                                                                                                                                                                              >  1:797022/1‑239 (MQ=255)
                                                                       ggatAACACGCCCCAAGGTCCCCGGCGGCGTCCGGCGCGGGTGTTCACGCCTTCAACGGTTGGCCCGGAGAGGGTGGTTAACACGTCTCCCGGACGATAGGCTCGGCCGCCAGGCTTGTTTCCGCGCCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCtt                                                                                                                                                                    <  2:98721/236‑1 (MQ=255)
                                                                                  cGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCCTCACGCCGTAAAc                                                                                                                                                                                                                 >  1:607885/1‑184 (MQ=255)
                                                                                               ggcggcCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATTTCGTACTTCATCTCCTCCATCCCTTCTGAAGGCTTGATCGAGATAcc                                                                                                                                              >  1:1321994/1‑238 (MQ=255)
                                                                                               ggcggcCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGAACGAGATAcc                                                                                                                                              >  1:552424/1‑238 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTTTAGCTCCGCGACCATCCGCATCACGCCGGCAACCGCTCACCGCCCGCACATATCCGACTTAACCTCATCCATGCCTTCTGAAGGCTTGACCGGAATAccgcccgg                                                                                                                                        >  2:189064/1‑236 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTCAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGc                                                                                                                                                                       >  1:406987/1‑208 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGGGATACCCCCGGa                                                                                                                                        >  1:1212176/1‑239 (MQ=255)
                                                                                                        cAGCGTGGGTGTTCAGCGCTTCAACGGTTTGACCGGAAGTGGTGGTAAACCACGCGCCCGGACGAAGGGCCGCGCCGCCAGGCTGGTTTCCCGACCCTGCCAAAACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCg                                                                                                                                    <  2:567934/239‑1 (MQ=255)
                                                                                                            accGGTGTTTAGCACTTTACAGGGTTTGGCCGGGCGGGGGGGTTAACCACACCCCCGGAGAGATGGTCCTGCCGCAAGGCATTTTGTCCCGGCCGCCCAACCCGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGGCCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGa                                                                                                                                   <  2:709334/235‑1 (MQ=255)
                                                                                                                            ttCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCCTCACGCCGTAAACCGCTGCCGCACCGCAAATATCGTACTTCATCTCATCCATGCCTCCTGAAGGCTTGATCGGGCTACCGCCGGAGTCGTAAGTTAAACCtccgtc                                                                                                                   >  1:709334/1‑232 (MQ=255)
                                                                                                                            ttCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATAACGCCGTAAACCGCTGCCGCCCCGAGCATATCGGACTTCATCTCCTCCCTCCCTTCTTAAGGCATGATCCGGACACCGCCCAAGTCAAAGGGTAAACATTTAcaca                                                                                                                 >  2:276634/1‑239 (MQ=255)
                                                                                                                            ttCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCACCACGCCGTAAACCGCTGCCGCACCGCACATATCGGCCGTCCTCTCATCCCTGCCTTCTGAAGGCTTGGTCGAGATACCGCCCGAGGCGAAGGTAAAACCTTTAccc                                                                                                                  >  2:1278996/1‑238 (MQ=255)
                                                                                                                            ttCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGGTGTAGCTCCGCGACCATCCGCATCACGCCGGAAACCGCTGCCCCCCCGCCCATATCGTACTTCATCTCATCCCAGCCTTCTGAAGGCCTGATCGAGATACCCCCGGAGCCCACGGGTAAACCTTTaacaa                                                                                                                 >  2:1130707/1‑235 (MQ=255)
                                                                                                                             ccAACGGTTTGACGGGACATGGTGGTTAACCCATCGCCCGGACGATGGGCCCGTCCGCCAGGAATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACccac                                                                                                                <  1:355827/238‑1 (MQ=255)
                                                                                                                                 ccGGTTTGACCGGACAGGGTGGTTAACACATCGCCCGGACGATAGGCTCGCCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAg                                                                                                             <  1:385306/238‑1 (MQ=255)
                                                                                                                                    gTTGGACGGGACAGGGTGGTTAAACCACCGCCGGGACGATAGGCTCGTCCGCCGGGCATGTTTCCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGc                                                                                                            <  1:24232/237‑1 (MQ=255)
                                                                                                                                       tGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGACGGCTTGATCGGGCTACCGCCGGCGTCGAAGGTTAAACCTTTACCCACCACCACGAt                                                                                                       >  1:612923/1‑239 (MQ=255)
                                                                                                                                       tGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTAATTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTCAACCTTTACccacacccccgat                                                                                                       >  1:1163761/1‑230 (MQ=255)
                                                                                                                                         aCCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACCTATCGTACTTCATCTCATCCATGCCTTCGTAAGGCTTGATCGAGATACCGCCGGAGTCGACGGGTCAACCTTTCCccac                                                                                                                >  1:1357784/1‑228 (MQ=255)
                                                                                                                                             gACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCCCTGCCGCACCGCACATATCGTACTTCATCTCATCCATTCCTTCTGCAGGCTTGCTCGAGAAACCGCCGGAGTCGAAAGTTAAACCTTTACCCCCCAGCACCATTGGGCg                                                                                                >  2:673088/1‑239 (MQ=255)
                                                                                                                                                            accacgccccGGACGATAGGCGTGGCCGACAGGCATTTTTACGCAGCGTGCCAACCAGCAGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTTCTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGAcg                                                                                 <  1:931764/232‑1 (MQ=255)
                                                                                                                                                               tccGCCAGGACGATGGCCTCGCCCGCCAGGCTTGTTTTCGCAGCCTCCCAACCGGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGt                                                                              <  1:708572/237‑1 (MQ=255)
                                                                                                                                                               aTCGCCCGGACGATGGGCTCGGCCGCCCGGCATGTTTGCGCCGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAACCCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGt                                                                              <  1:668850/239‑1 (MQ=255)
                                                                                                                                                               aTCCCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAAACGGCAGTTGTAGCGCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGt                                                                              <  1:1242146/239‑1 (MQ=255)
                                                                                                                                                                 ggCCCGGCCGATAGGCTGGTCCGCCAGGCATGTTTTCGCAGCCTTCCAACACGCCGATAACGTTAATCGGCCGTTGTAGCTCCGCGCCCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTg                                                                            <  1:579589/238‑1 (MQ=255)
                                                                                                                                                                 cGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTCATCCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTTATCGAGATACCGCCGGAGTGCAAGGGTAAACCTTTACCCACCACCCCCATTGGGCGGGCACTTTCCgcccggtg                                                                             >  2:1266591/1‑231 (MQ=255)
                                                                                                                                                                 cGCCCGGACGATAGGCTCGTCCACCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCACCGGAGTCGAAAGTTAAACCTTTACCCCCCAGCCCCCTTGGGCGGGGATCTTCCGAcgcggt                                                                             >  2:1002150/1‑236 (MQ=255)
                                                                                                                                                                  gCCCGGACGTAGGGCTCTGCCGCCAGGCTTGGTTCCGCGGCCTGCCAACCCGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGt                                                                              <  1:541854/236‑1 (MQ=255)
                                                                                                                                                                       gACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCCCGATTGGGCGTGCATCTTCCCACCCGTTGCCtttt                                                                      >  1:410573/1‑238 (MQ=255)
                                                                                                                                                                            tAGGCTCGTCCGCCGGGGTGTTTTCCCCCCCTGCCAACCCCGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGGAAACCGGTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTc                                                                 <  2:1321994/239‑1 (MQ=255)
                                                                                                                                                                                  cGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCCCCGCGACCATCCGCCTCACGCCGGAAACCGCCGCCGCACCGCCCAAATCGTACTTCCTCTCATCCCTGCCTTCTGAAGGCCTGATCGCGATACCGCCGGAGCCGAAGGTTGAACCTTTACCCACCAGCCCGATGGGTCTGCACCTTTCCGCCCCTTTGCCTTTTTACTCAAacac                                                           >  2:853866/1‑239 (MQ=255)
                                                                                                                                                                                         ccAGGCATTTTTCCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACAt                                                     <  1:1159401/239‑1 (MQ=255)
                                                                                                                                                                                           gcgcttttttcCGCCGCCCCCCAACCCACCGATAACGTTAATCGGCAGTTGTAGCCCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCCCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCa                                                   <  2:579313/232‑1 (MQ=255)
                                                                                                                                                                                                        accccTGCCCCCACCCCGATTCGGTTAATCGGCATTTGTATCTCCGCGACCATCCCCATCCCGCCGTAAACCGCTGCCGCACCGCACATACCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTg                                      <  2:884369/236‑1 (MQ=255)
                                                                                                                                                                                                              gCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGCGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCCCCCGCACGAATGGGCGGTGCTCTTCCCACGCGTGGCCTTTTTACCCAACCACCCACACCAGCGGTTCGTTTTGCCaaa                                 >  2:1071804/1‑237 (MQ=255)
                                                                                                                                                                                                                     ccGCCGAAAACGTTAAGCGGGAGTTGTAGCTCCGCGACCTTCCGCATCACGCCGTAAACCGCTGCCGCACCTCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCGgacc                          <  1:1010038/238‑1 (MQ=255)
                                                                                                                                                                                                                     acgcgtaTAACGTTAATCGGCCGTTGAAGCTCCACGACCATCCGCAACACGCCGTAAACCGCTGCCGCCACGCACATATCGTACTTCATCGCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgacc                          <  1:799495/233‑1 (MQ=255)
                                                                                                                                                                                                                        tcgttAACGTTATCCGGCGGTTGTAGCTCCGCGCCCATCCGCATCACGCCGTAACCCCCTGCCGCCCCCCCCATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgac                       <  2:72624/235‑1 (MQ=255)
                                                                                                                                                                                                                         ccgtaacagTTAATCGGCCGTTGACGCCCCCCCACCACCCCCACCCCGCCGTAACCCGCTGCCGCACCGCACATTTCGTACTTCAGCCCGTCCATGCCTTCGGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCGTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgac                       <  2:1033229/231‑1 (MQ=255)
                                                                                                                                                                                                                         cGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGAACTCAATCACCGACATCAGCGATTCGTTTTCCGAACCCTgaccgac                       >  1:623475/1‑238 (MQ=255)
                                                                                                                                                                                                                         cGATAAAGTTAATCGGCAGTTGTAGTCCCGCGACCCTCCGCATCACGCCGTAAACCGGTGCCGCACCGCACATATCGTAATTCATCCCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgac                       <  1:1247278/238‑1 (MQ=255)
                                                                                                                                                                                                                          ggtAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGGTGCCGCACCGAACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGaccgaccg                     <  1:1071804/237‑1 (MQ=255)
                                                                                                                                                                                                                          gATAACGTTAATCGGCAGTTTTAGCTCCGCGACCATCCGCATCCCGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGaccgaccg                     <  1:926985/239‑1 (MQ=255)
                                                                                                                                                                                                                                   tAACCGGCAGTTTTACCTCCGGGACCATCCGCATCACGCCGTAAACCGCCGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCCGACCGGCCGCCAGATa              <  1:1002150/238‑1 (MQ=255)
                                                                                                                                                                                                                                             ttgccgCTCCGCAGCCATCCGCCTCACGCCGTAAACCGCTGCCGCCCCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc   <  2:1212176/234‑1 (MQ=255)
                                                                                                                                                                                                                                              tggaGCCCCCCGACCATCCGCATCACGCCTTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATccc  <  1:801001/236‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           
AGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCC  >  NC_000913/4484789‑4485258

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