Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,485,027 T→C 13.3% 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,0270TC13.3% 148.7 / 3.7 73E309G (GAG→GGG) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base T (28/33);  new base C (5/6);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCC  >  NC_000913/4484790‑4485259
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
aGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGCCGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCCGTTTGAGCt                                                                                                                                                                                                                                              >  1:1309048/1‑238 (MQ=255)
  cgcgATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGAGAACGTTAATCGGCAGTTTTAGCTcc                                                                                                                                                                                                                                            >  1:953386/1‑238 (MQ=255)
      aTCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGCCAGTTGAAGCCCCGCGa                                                                                                                                                                                                                                        >  1:884369/1‑238 (MQ=255)
                 cccggtcaccgcccccgctcaactCACCTCTCCCGCCTCAAAACGCCCCACGTAAGTTAACGCGCCCCACAGTACCAGGCGTCTTTCAGCTTGGGTCTTCCGCCCTTAAACGTTTTGCCCGGACAGGTTGTTTAACACCTCCCCCGGCCGATAGGCTCGTCCGCCAGGCATGTTTCCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATc                                                                                                                                                                                                                             <  2:953386/215‑1 (MQ=255)
                  cgggtccgtgccccaaccacAATCACGGGTGCGGGCTAAAACCTCTCAGCGAAATTTAACGCGTCGCACGGTACGAGGCGGCCTTAAGTGTCGGTTTTCAGCAATTCAACTGTGTGGCCGGACATGGTGGTTAACACATCGCCCGGAGGATAGGCTCGTCCGCCAGGCTTGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCAc                                                                                                                                                                                                                           <  2:829607/220‑1 (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)
                                         ccccccttccgcgcaaaaaCCCTCAACGTAAGTTACCACGCCGCACAGCCCCAGGCGGCCTTCAGCGTGGGTTTTCAGCCCTTCACCGGTTTGACGGGACATGGTGTTTACCACACCGCCCGGACGATAGGCTCGTCCGCCCGGCATGTTTTCGCAGCCTGCCACCACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGccgcaccgc                                                                                                                                                                                                   <  2:914243/224‑1 (MQ=255)
                                             gCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCCGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACat                                                                                                                                                                                               >  1:797022/1‑239 (MQ=255)
                                                                     gtagttaaaaccgtcaccacgTACCGGGGGGCCTTCGACGTCGGGTGTCGGCCTTTCAACGGTTTGCCCGGACATGGTGTTTAACCACTCGCCCGGGCGATAGGTCCGTCCGCCAGGCATGTTTTCGCAGCCTACCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGc                                                                                                                                                                        <  2:1021193/219‑1 (MQ=255)
                                                                       ggatAACACGCCCCAAGGTCCCCGGCGGCGTCCGGCGCGGGTGTTCACGCCTTCAACGGTTGGCCCGGAGAGGGTGGTTAACACGTCTCCCGGACGATAGGCTCGGCCGCCAGGCTTGTTTCCGCGCCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCtt                                                                                                                                                                     <  2:98721/236‑1 (MQ=255)
                                                                                cgtcccccaggtccggggcgcCTTCCGGGTGGGTGTTCAGCCATTCAACGGTTGGACCGGAACTGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCGGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGc                                                                                                                                                             <  2:657399/220‑1 (MQ=255)
                                                                                  cGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCCTCACGCCGTAAAc                                                                                                                                                                                                                  >  1:607885/1‑184 (MQ=255)
                                                                                               ggcggcCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATTTCGTACTTCATCTCCTCCATCCCTTCTGAAGGCTTGATCGAGATAcc                                                                                                                                               >  1:1321994/1‑238 (MQ=255)
                                                                                               ggcggcCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGAACGAGATAcc                                                                                                                                               >  1:552424/1‑238 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCACCACGCCGATAACGTTAATCGGCAGTTGTAACTCCGCGACCATCCGCACCCACCCCTAAACCGCTGGCCCCCCCCACACACCGTACTTCCCTCCATCCCTACCTTCTGACGGccttgacgaataacccccccg                                                                                                                                         >  2:1242146/1‑219 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTTTAGCTCCGCGACCATCCGCATCACGCCGGCAACCGCTCACCGCCCGCACATATCCGACTTAACCTCATCCATGCCTTCTGAAGGCTTGACCGGAATAccgcccgg                                                                                                                                        >  2:189064/1‑236 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTCAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGc                                                                                                                                                                        >  1:406987/1‑208 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGGGATACCCCCGGa                                                                                                                                        >  1:1212176/1‑239 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCCCTGCCGCACCGCCCATATCGTACTTCATCCCATCCATGCCTTCTGAAGGCTTgtccggaaaaaccccggg                                                                                                                                         >  2:1159401/1‑222 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCCCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCCCCGCCCAACACCGACTTCATCTCATCCCTGCCTTCCAAAAGCCTGAtacaaaaaaccccgga                                                                                                                                         >  2:579589/1‑224 (MQ=255)
                                                                                                        cAGCGTGGGTGTTCAGCGCTTCAACGGTTTGACCGGAAGTGGTGGTAAACCACGCGCCCGGACGAAGGGCCGCGCCGCCAGGCTGGTTTCCCGACCCTGCCAAAACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCg                                                                                                                                    <  2:567934/239‑1 (MQ=255)
                                                                                                            accGGTGTTTAGCACTTTACAGGGTTTGGCCGGGCGGGGGGGTTAACCACACCCCCGGAGAGATGGTCCTGCCGCAAGGCATTTTGTCCCGGCCGCCCAACCCGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGGCCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGa                                                                                                                                   <  2:709334/235‑1 (MQ=255)
                                                                                                                    aaaggagcatcccaCGGTTTGCCGGGTCTGGGGGCTAACCCATCGCCCCGACCGTTGGCTCTGCCCCTCAGCCTGGTTTACGAACCCGGCCACCCCGCCGATAACGTTAATCGGCAGTTGTAGCCCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAAc                                                                                                                          <  2:1357784/226‑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)
                                                                                                                                                               aTCTCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCTCGACCATCCGCATCCCGCCGTAAACCGCTGCCGCACCCGCCATATCGCACCTCCTCTCCTCCATGCCTTCTGAACGCCTGACAGAGATACCGCCCGGCGTCGAGGTTAAACCTTTACCCCCCCACACATTTGGGCGTggaccctcccgccggt                                                                               >  2:1121769/1‑224 (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)
                                                                                                                                                                 cGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCCCCCATATCGTACTTCATCTCATCCCTGCCTTCTGAAGGCTTGATCGGGATACCGCCGGCGTCGAAGGGTAAACCCTTAACCCCCACCACGAtttgtggggacccttccgccccgtt                                                                               >  2:637901/1‑215 (MQ=255)
                                                                                                                                                                 cGCCCGGACGATAGGCTCGTCCACCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCACCGGAGTCGAAAGTTAAACCTTTACCCCCCAGCCCCCTTGGGCGGGGATCTTCCGAcgcggt                                                                             >  2:1002150/1‑236 (MQ=255)
                                                                                                                                                                  gCCCGGACGTAGGGCTCTGCCGCCAGGCTTGGTTCCGCGGCCTGCCAACCCGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGt                                                                              <  1:541854/236‑1 (MQ=255)
                                                                                                                                                                  gCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCCCCGCGCCCCTCCGCATCACGCCGTAAACCGCAGCCGCACCGCACATATCGTACTTCATCTCATCCCTGCCCTCTGAAGGCTTGATCGGGATACCTCCGGGCTCGAAGGTTTAAACTTTACCCCCCACCACCATTGGGCGTGCAtttctccccgcgtc                                                                              >  2:1295726/1‑224 (MQ=255)
                                                                                                                                                                       gACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCCCGATTGGGCGTGCATCTTCCCACCCGTTGCCtttt                                                                      >  1:410573/1‑238 (MQ=255)
                                                                                                                                                                            tAGGCTCGTCCGCCGGGGTGTTTTCCCCCCCTGCCAACCCCGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGGAAACCGGTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTc                                                                 <  2:1321994/239‑1 (MQ=255)
                                                                                                                                                                              gcggctgcgcgcacAGGCATTTTTGCGACACCCGCCAACACGCGGATACCGTTAGTCGGCAGTGGTAGCCCCTGTCCCATCCGAAACACGCCGTAAACCGCTGCCGGACCGCACATATCGGGTGTCATCCCATCCAGGCCTTCTGAAGGCTTTATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACAAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACt                                                                  <  2:1107662/225‑1 (MQ=255)
                                                                                                                                                                                  cGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCCCCGCGACCATCCGCCTCACGCCGGAAACCGCCGCCGCACCGCCCAAATCGTACTTCCTCTCATCCCTGCCTTCTGAAGGCCTGATCGCGATACCGCCGGAGCCGAAGGTTGAACCTTTACCCACCAGCCCGATGGGTCTGCACCTTTCCGCCCCTTTGCCTTTTTACTCAAacac                                                           >  2:853866/1‑239 (MQ=255)
                                                                                                                                                                                     tCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGCTAACGTTAATCGGCAGTTGTAGCTCCGCGTCCATCCGCATCCCGCCGTAAACCGCTGCCGCACCGCCCAAATCGCACTTCATCTCCTGCATGCCTTCTGAAGGCCTGAACGAGAGACCGCCCGCATCGAAGGTTAAACCCTTACCCACCCCCCCACTTGGGCGTGCCCCTTTCAAACCGTTGCAttttaacccaaatcccg                                                         >  2:299767/1‑225 (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)
                                                                                                                                                                                                                                     aatcggccagtgatgCTCCGCGACCACCCGCAGCACGCCGTAAACCGCTGCCGCACCGCACATAGCGGACTTCCTCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAgg            <  1:276634/225‑1 (MQ=255)
                                                                                                                                                                                                                                              ttgccgCTCCGCAGCCATCCGCCTCACGCCGTAAACCGCTGCCGCCCCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc   <  2:1212176/234‑1 (MQ=255)
                                                                                                                                                                                                                                               tggaGCCCCCCGACCATCCGCATCACGCCTTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATccc  <  1:801001/236‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
AGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTCGCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGCACTTCAACGGTTTGACCGGACATGGTGGTTAACACATCGCCCGGACGATAGGCTCGTCCGCCAGGCATGTTTTCGCAGCCTGCCAACACGCCGATAACGTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCC  >  NC_000913/4484790‑4485259

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