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

AGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTC‑GCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCC  >  NC_000913/4484789‑4485258
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           
aGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTC‑GCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGCCGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCCGTTTGAGCt                                                                                                                                                                                                                                             >  1:1309048/1‑238 (MQ=255)
  cgcgATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTC‑GCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGAGAAC‑GTTAATCGGCAGTTTTAGCTcc                                                                                                                                                                                                                                           >  1:953386/1‑238 (MQ=255)
      aTCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTC‑GCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGCCAGTTGAAGCCCCGCGa                                                                                                                                                                                                                                       >  1:884369/1‑238 (MQ=255)
                          cGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTACGTTAACACGTC‑GCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCACGCATGTTTTCTCAGCCTGCCAA‑CACGCCGATTAC‑GTTAATCGTCAGTTGAACTCCCGCGGCCATCCGCATCACGCCGTaa                                                                                                                                                                                                                   >  1:1107662/1‑238 (MQ=255)
                           gTCGCCACTCCAATCCCCGCTCCCGCCTAAAAACGCTCAACGTAAGTTAACACGTC‑GCACAGTCCCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CCCGTCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAAc                                                                                                                                                                                                                 <  2:830860/239‑1 (MQ=255)
                            tCGCCACGTCAATCCCCGCTTCCGGCCAAAAACGCAAAACGTAAGTTAACCCGGC‑GCACAGTCCCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGC‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAAc                                                                                                                                                                                                                 <  1:609420/238‑1 (MQ=255)
                               ccACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTC‑GCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTTTAGCTCCGCGACCATCCGCCTCACCCCGTAAACCGc                                                                                                                                                                                                              >  1:657399/1‑238 (MQ=255)
                               ccACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTC‑GCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCACGCATGTTTTCACAGCCTGCCAA‑CACGCCG‑TACCAGTTACTCGGCAGTTGTAGCCCCCCGCCCAACCCCCCCACCCCGCaaaacccc                                                                                                                                                                                                             >  2:679535/1‑234 (MQ=255)
                                   gTCAATCCCCGCTTCCGGCTCAAACCGCTCAACGGAAGTTAACCCGTC‑GCACAGTCCCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑Gccg                                                                                                                                                                                                        <  1:166036/239‑1 (MQ=255)
                                    ccaacacccccTTCCGGCCAAAAAGCCTCAACGGAAGTTAAC‑CGCCACCATAATACCGGGCGGCTTTCAGGGCCGGTG‑TCAGCGCCTTAA‑AGGGTTTGACCGAACAGGGTGGTTAA‑CACCGCCCCGGGACGATAGGCTCGT‑CCGCCAGGCATTTTTTCGCAGCCTGCAAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑Gccgc                                                                                                                                                                                                       <  2:1143684/229‑1 (MQ=255)
                                      aaTCCCCGCTCCCGTCTAAAACCCCCCAACGTAAGTTAACACGTC‑CACCAGTCCCAGGGGGCCTTACGGGTGGGTGTTCGGC‑ACTTCA‑ACGGTTTGCCCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCAAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑Gccgcac                                                                                                                                                                                                     <  2:413641/239‑1 (MQ=255)
                                         ccccctttCCGGCCCAAACCCCTCAAGGAAATTAAACCGGGC‑CGACGGTCCCAGGCGGCCTTCGGCGTCGGTGTTAGGC‑ACTACA‑GCGGTTTGACCGGACGTGGGGGTTAA‑CAAATCGCCTGGAGGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCCCCCGTCCAC‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑Gccgcaccgc                                                                                                                                                                                                  <  2:1235598/233‑1 (MQ=255)
                                             gCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTC‑GCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCCGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACat                                                                                                                                                                                              >  1:797022/1‑239 (MQ=255)
                                                                       ggatAACACGCC‑CCAAGGTCCCCGGCGGCGTCCGGCGCGGGTGTTCACG‑CCTTCA‑ACGGTTGGCCCGGAGAGGGTGGTTAA‑CACGTCTCCCGGACGATAGGCTCGG‑CCGCCAGGCTTGTTTCCGCGCCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCtt                                                                                                                                                                    <  2:98721/236‑1 (MQ=255)
                                                                                  c‑GCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCCTCACGCCGTAAAc                                                                                                                                                                                                                 >  1:607885/1‑184 (MQ=255)
                                                                                               ggcggcCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATTTCGTACTTCATCTCCTCCATCCCTTCTGAAGGCTTGATCGAGATAcc                                                                                                                                              >  1:1321994/1‑238 (MQ=255)
                                                                                               ggcggcCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGAACGAGATAcc                                                                                                                                              >  1:552424/1‑238 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTTTAGCTCCGCGACCATCCGCATCACGCCGGCAACCGCTCACCGC‑CCGCACATATCCGACTTAACCTCATCCATGCCTTCTGAAGGCTTGACCGGAATAccgcccgg                                                                                                                                        >  2:189064/1‑236 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTCAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGc                                                                                                                                                                       >  1:406987/1‑208 (MQ=255)
                                                                                                    ccTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGGGATACCCCCGGa                                                                                                                                        >  1:1212176/1‑239 (MQ=255)
                                                                                                        cAGCGTGGGTGTTCAGC‑GCTTCA‑ACGGTTTGACCGGAAGTGGTGGTAAA‑CCACGCGCCCGGACGAAGGGCCGCG‑CCGCCAGGCTGGTTTCCCGACCCTGCCAA‑AACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCg                                                                                                                                    <  2:567934/239‑1 (MQ=255)
                                                                                                            accGGTGTTTAGC‑ACTTTACAGGGTTTGGCCGGGCGGGGGGGTTAACCACA‑CCCCCGGA‑GAGATGGTCCTGCCGCAAGGCATTTTGTCCCGGCCGCCCAA‑CCCGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGGCCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGa                                                                                                                                   <  2:709334/235‑1 (MQ=255)
                                                                                                                            ttCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCCTCACGCCGTAAACCGCT‑GCCGCACCGCAAATATCGTACTTCATCTCATCCATGCCTCCTGAAGGCTTGATCGGGCTACCGCCGGAGTCGTAAGTTAAACCtccgtc                                                                                                                   >  1:709334/1‑232 (MQ=255)
                                                                                                                            ttCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATAACGCCGTAAACCGCT‑GCCGCCCCGAGCATATCGGACTTCATCTCCTCCCTCCCTTCTTAAGGCATGATCCGGACACCGCCCAAGTCAAAGGGTAAACATTTAcaca                                                                                                                 >  2:276634/1‑239 (MQ=255)
                                                                                                                            ttCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCACCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGGCCGTCCTCTCATCCCTGCCTTCTGAAGGCTTGGTCGAGATACCGCCCGAGGCGAAGGTAAAACCTTTAccc                                                                                                                  >  2:1278996/1‑238 (MQ=255)
                                                                                                                            ttCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGGTGTAGCTCCGCGACCATCCGCATCACGCCGGAAACCGCT‑GCCCCCCCGCCCATATCGTACTTCATCTCATCCCAGCCTTCTGAAGGCCTGATCGAGATACCCCCGGAGCCCACGGGTAAACCTTTaacaa                                                                                                                 >  2:1130707/1‑235 (MQ=255)
                                                                                                                             ccA‑ACGGTTTGACGGGACATGGTGGTTAA‑CCCATCGCCCGGACGATGGGCCCGT‑CCGCCAGGAATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACccac                                                                                                                <  1:355827/238‑1 (MQ=255)
                                                                                                                                 ccGGTTTGACCGGACAGGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGC‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAg                                                                                                             <  1:385306/238‑1 (MQ=255)
                                                                                                                                    gTTGGACGGGACAGGGTGGTTAA‑ACCACCGCCGGGACGATAGGCTCGT‑CCGCCGGGCATGTTTCCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGc                                                                                                            <  1:24232/237‑1 (MQ=255)
                                                                                                                                       tGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGACGGCTTGATCGGGCTACCGCCGGCGTCGAAGGTTAAACCTTTACCCACCACCACGAt                                                                                                       >  1:612923/1‑239 (MQ=255)
                                                                                                                                       tGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTAATTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTCAACCTTTACccacacccccgat                                                                                                       >  1:1163761/1‑230 (MQ=255)
                                                                                                                                         aCCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACCTATCGTACTTCATCTCATCCATGCCTTCGTAAGGCTTGATCGAGATACCGCCGGAGTCGACGGGTCAACCTTTCCccac                                                                                                                >  1:1357784/1‑228 (MQ=255)
                                                                                                                                             gACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCCCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATTCCTTCTGCAGGCTTGCTCGAGAAACCGCCGGAGTCGAAAGTTAAACCTTTACCCCCCAGCACCATTGGGC‑g                                                                                                >  2:673088/1‑239 (MQ=255)
                                                                                                                                                            accacgccccGGACGATAGGCGTGG‑CCGACAGGCATTTTTACGCAGCGTGCCAA‑CCAGCAGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTTCTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGAcg                                                                                 <  1:931764/232‑1 (MQ=255)
                                                                                                                                                               tccGCCAGGACGATGGCCTCGC‑CCGCCAGGCTTGTTTTCGCAGCCTCCCAA‑CCGGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGt                                                                              <  1:708572/237‑1 (MQ=255)
                                                                                                                                                               aTCGCCCGGACGATGGGCTCGG‑CCGCCCGGCATGTTTGCGCCGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAACCCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGt                                                                              <  1:668850/239‑1 (MQ=255)
                                                                                                                                                               aTCCCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAAACGGCAGTTGTAGCGCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGt                                                                              <  1:1242146/239‑1 (MQ=255)
                                                                                                                                                                 ggCCCGGCCGATAGGCTGGT‑CCGCCAGGCATGTTTTCGCAGCCTTCCAA‑CACGCCGATAAC‑GTTAATCGGCCGTTGTAGCTCCGCGCCCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTg                                                                            <  1:579589/238‑1 (MQ=255)
                                                                                                                                                                 cGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTCATCCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTTATCGAGATACCGCCGGAGTGCAAGGGTAAACCTTTACCCACCACCCCCATTGGGC‑GGGCACT‑TTCCgcccggtg                                                                             >  2:1266591/1‑231 (MQ=255)
                                                                                                                                                                 cGCCCGGACGATAGGCTCGT‑CCACCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCACCGGAGTCGAAAGTTAAACCTTTACCCCCCAGCCCCCTTGGGC‑GGGGATC‑TTCCGAcgcggt                                                                             >  2:1002150/1‑236 (MQ=255)
                                                                                                                                                                  gCCCGGACGTAGGGCTCTG‑CCGCCAGGCTTGGTTCCGCGGCCTGCCAA‑CCCGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGt                                                                              <  1:541854/236‑1 (MQ=255)
                                                                                                                                                                       gACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCCCGATTGGGC‑GTGCATC‑TTCCCACCCGTTGCCtttt                                                                      >  1:410573/1‑238 (MQ=255)
                                                                                                                                                                            tAGGCTCGT‑CCGCC‑GGGGTGTTTTCCCCCCCTGCCAACCCCGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGGAAACCGGT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTc                                                                 <  2:1321994/239‑1 (MQ=255)
                                                                                                                                                                                  cGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCCCCGCGACCATCCGCCTCACGCCGGAAACCGCC‑GCCGCACCGCCCAAATCGTACTTCCTCTCATCCCTGCCTTCTGAAGGCCTGATCGCGATACCGCCGGAGCCGAAGGTTGAACCTTTACCCACCAGCCCGA‑TGGGT‑CTGCACCTTTCCGCCCCTTTGCCTTTTTACTCAAacac                                                           >  2:853866/1‑239 (MQ=255)
                                                                                                                                                                                         ccAGGCATTTTTCCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACAt                                                     <  1:1159401/239‑1 (MQ=255)
                                                                                                                                                                                           gcgcttttttcCGCCGCCCCCCAA‑CCCACCGATAAC‑GTTAATCGGCAGTTGTAGCCCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCCCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCa                                                   <  2:579313/232‑1 (MQ=255)
                                                                                                                                                                                                        accccTGCCCC‑CACCCCGATTCG‑GTTAATCGGCATTTGTATCTCCGCGACCATCCCCATCCCGCCGTAAACCGCT‑GCCGCACCGCACATACCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTg                                      <  2:884369/236‑1 (MQ=255)
                                                                                                                                                                                                              gCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGCGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCCCCCGCACGAATGGGCGGTGC‑TC‑TTCCCACGCGTGGCCTTTTTACCCAACCACCCACACCAGCGGTTCGTTTTGCCaaa                                 >  2:1071804/1‑237 (MQ=255)
                                                                                                                                                                                                                     ccGCCGAAAAC‑GTTAAGCGGGAGTTGTAGCTCCGCGACCTTCCGCATCACGCCGTAAACCGCT‑GCCGCACCTCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCGgacc                          <  1:1010038/238‑1 (MQ=255)
                                                                                                                                                                                                                     acgcgtaTAAC‑GTTAATCGGCCGTTGAAGCTCCACGACCATCCGCAACACGCCGTAAACCGCT‑GCCGCCACGCACATATCGTACTTCATCGCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgacc                          <  1:799495/233‑1 (MQ=255)
                                                                                                                                                                                                                        tcgttAAC‑GTTATCCGGCGGTTGTAGCTCCGCGCCCATCCGCATCACGCCGTAACCCCCT‑GCCGCCCCCCCCATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgac                       <  2:72624/235‑1 (MQ=255)
                                                                                                                                                                                                                         ccgtaacagTTAATCGGCCGTTGACGCCCCCCCACCACCCCCACCCCGCCGTAACCCGCT‑GCCGCACCGCACATTTCGTACTTCAGCCCGTCCATGCCTTCGGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCGTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgac                       <  2:1033229/231‑1 (MQ=255)
                                                                                                                                                                                                                         cGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGAACTCAATCACCGACATCAGCGATTCGTTTTCCGAACCCTgaccgac                       >  1:623475/1‑238 (MQ=255)
                                                                                                                                                                                                                         cGATAAA‑GTTAATCGGCAGTTGTAGTCCCGCGACCCTCCGCATCACGCCGTAAACCGGT‑GCCGCACCGCACATATCGTAATTCATCCCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgac                       <  1:1247278/238‑1 (MQ=255)
                                                                                                                                                                                                                          ggtAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGGT‑GCCGCACCGAACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGaccgaccg                     <  1:1071804/237‑1 (MQ=255)
                                                                                                                                                                                                                          gATAAC‑GTTAATCGGCAGTTTTAGCTCCGCGACCATCCGCATCCCGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGaccgaccg                     <  1:926985/239‑1 (MQ=255)
                                                                                                                                                                                                                                   tAACCGGCAGTTTTACCTCCGGGACCATCCGCATCACGCCGTAAACCGCC‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCCGACCGGCCGCCAGATa              <  1:1002150/238‑1 (MQ=255)
                                                                                                                                                                                                                                             ttgccgCTCCGCAGCCATCCGCCTCACGCCGTAAACCGCT‑GCCGCCCCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc   <  2:1212176/234‑1 (MQ=255)
                                                                                                                                                                                                                                              tggaGCCCCCCGACCATCCGCATCACGCCTTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATccc  <  1:801001/236‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           
AGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGTAAGTTAACACGTC‑GCACAGTACCAGGCGGCCTTCAGCGTCGGTGTTCAGC‑ACTTCA‑ACGGTTTGACCGGACATGGTGGTTAA‑CACATCGCCCGGACGATAGGCTCGT‑CCGCCAGGCATGTTTTCGCAGCCTGCCAA‑CACGCCGATAAC‑GTTAATCGGCAGTTGTAGCTCCGCGACCATCCGCATCACGCCGTAAACCGCT‑GCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGC‑GTGCATC‑TTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCC  >  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: ‑