Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,436 T→C T38T (ACT→ACC)  deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,4360TC98.6% 200.7 / ‑3.5 70coding (114/780 nt)deoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (49/20);  minor base . (1/0);  total (50/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

GCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTGATCGCC‑CTGTGTCATCAGGCCAA‑AACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT  >  NC_000913/4617190‑4617679
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   
gCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATTGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGAAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCgg                                                                                                                                                                                                                                                 >  1:79280/1‑251 (MQ=255)
           aTGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATAcc                                                                                                                                                                                                                                      <  1:155285/251‑1 (MQ=255)
                     gTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGCGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCCGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGACTGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCt                                                                                                                                                                                                                            >  1:596364/1‑251 (MQ=255)
                                 ccTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTctcctc                                                                                                                                                                                                                >  1:514669/1‑251 (MQ=255)
                                     aattgtgatgtagtcctaagggtgtTGGGGAGTAGATGGGAGAATACGAAGAACCTCGCAAG‑TGTATTTTTAGTGGCGACAAGCTAGGAAAAGGAAATAACTGATCTGAAAGCAA‑GCACCCGGCTTGCACGGAAATTGTTGGAGCTTACCTCCCTGAAGGTCGACGACACCGTC‑AGTAAGGTAATTCCT‑TTGTGTCATCGGGCGAA‑AACCCCTGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCt                                                                                                                                                                                                             <  1:409477/229‑1 (MQ=255)
                                        tgtgATGTATATCGAAGTGTGTTGCGAAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATATGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTAt                                                                                                                                                                                                         >  2:536895/1‑251 (MQ=255)
                                           gATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATccc                                                                                                                                                                                                      >  1:143696/1‑251 (MQ=255)
                                           gATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGAAAGGTG‑ACTTTTATTGGCGACAAGACAGGCGAATGAAATGCCTGATCTGAAAGCAA‑GCAGCCTACGTGCACTGAAATTGCTGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AACGTAATTGCT‑CTGTGTCATCAGGCAAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATccc                                                                                                                                                                                                      >  1:595581/1‑251 (MQ=255)
                                               tatATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGAtt                                                                                                                                                                                                  >  1:289964/3‑251 (MQ=255)
                                                     gAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACACACTGGCCCGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGAAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGc                                                                                                                                                                                            >  2:43996/1‑251 (MQ=255)
                                                           ggtgggggAGTAAATGTTAAAATACAAAAAAACTCGCAAGGTT‑AATTTTATTGGCGACAAGCCAGGAGAAAGAAATGACTGATCTGAAAGCAA‑GCAGTCTGCGTGCACTGAAATGAATGGCCCTGCCCCCCCTGAATGACGCCGACACCGACGAA‑AAAGTATTTGCT‑CTGTTTCATCAGGCCAA‑ACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCTCTTTATCCCGATTGCTCGCAAAAca                                                                                                                                                                                      <  2:327131/245‑2 (MQ=255)
                                                                  gAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGaa                                                                                                                                                                                 >  1:58940/1‑249 (MQ=255)
                                                                    gTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAag                                                                                                                                                                               >  1:444435/1‑249 (MQ=255)
                                                                       gATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGc                                                                                                                                                                            >  1:329654/1‑249 (MQ=255)
                                                                          gTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATAcc                                                                                                                                                                                                                                      <  2:190875/188‑1 (MQ=255)
                                                                           ttAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACAGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGc                                                                                                                                                                       >  1:1075/1‑250 (MQ=255)
                                                                              gAATAATAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTCTCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCAccc                                                                                                                                                                   >  2:614593/1‑251 (MQ=255)
                                                                               aaataataaaaaactgaaaagggaaaTTTTATTGGCAACAAGCCAGAAGAATGAAATGTCTGATC‑GAAAACAAACCAGCCGGCGTGCCCGGAAATTGCTGGACCTGACCACCCTGAATGACCACACCCCCGCAGAG‑AAAGTAATTGCT‑CTGTTTAAACAGCCCAA‑ACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCAccc                                                                                                                                                                   <  2:888/227‑1 (MQ=255)
                                                                               aaTACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACCAGCCAGGAGCATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAC‑ACCCCCGGTCGGCACTACCGCCGCTATCTGTCTCTATCCTCGCTTTATCCCGATTGCTCGCAACACACTGAAAGCTCAGGTCAcccc                                                                                                                                                                  >  2:55042/1‑251 (MQ=255)
                                                                               aaTACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCAcccc                                                                                                                                                                  >  1:434162/1‑251 (MQ=255)
                                                                                     aaCAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAAt                                                                                                                                                            >  1:82453/1‑251 (MQ=255)
                                                                                        aaaCTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCg                                                                                                                                                         >  1:16728/1‑251 (MQ=255)
                                                                                          aacccgaaagggaaaTTTTATTGGGAACAAGCCAGAAAAATGAAATGATGGAACGGAAAGCAA‑GCAGCCTGCGTGCACTAAAATTGAGGGCCCTGACCCCCCTAAAGGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCGGGCAAA‑ACCCCCGGGCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGt                                                                                                                                                        <  2:257466/238‑1 (MQ=255)
                                                                                          aCTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGCAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAACTCCGt                                                                                                                                                        >  2:523174/1‑250 (MQ=255)
                                                                                           cTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTAt                                                                                                                                                      >  1:439727/1‑251 (MQ=255)
                                                                                           atCGCACGGTG‑AGTTTTATTGGCGACCAGCCAGGCGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTAt                                                                                                                                                      >  2:125960/2‑251 (MQ=255)
                                                                                               cAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCATGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCt                                                                                                                                                  >  1:222799/1‑251 (MQ=255)
                                                                                                aaGGTG‑AATTTTATGGGCGACAAGCAAGGAGAATGAAATGAATGATCTGAAAGCAA‑GACGCCGGCGGGCGCTGAAATGGATGGCCCTGACCCCCCTGAAGGACAACGACACCGACGAG‑AAAGTAAT‑GCCTCGGGGTCATCAGGCCAA‑ACCCCGGGCCGGAAATACCGCCGCTATCTGTATCGATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTa                                                                                                                                                 <  1:614593/251‑1 (MQ=255)
                                                                                                  ggTG‑AATTTTATTGGCGACAAGCCAGGAGATTGAAAGGACGGATCGAAAAGAAA‑GCAGCCTGCGTGCACTGAATGTGATGGACCTGACCCCCCTGAATGACGACGACCCCGACGAG‑AAAGTAATTGCT‑CTGTGTCTTAAGGCAAA‑AACCCCGTTCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACg                                                                                                                                               <  1:332094/251‑1 (MQ=255)
                                                                                                        aTTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAc                                                                                                                                          >  1:494934/1‑251 (MQ=255)
                                                                                                                    aCAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCcac                                                                                                                               >  1:108032/1‑250 (MQ=255)
                                                                                                                     cAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCACAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCacac                                                                                                                             >  2:417511/1‑251 (MQ=255)
                                                                                                                      aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                                                                            >  1:120103/1‑251 (MQ=255)
                                                                                                                      aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                                                                            >  1:316608/1‑251 (MQ=255)
                                                                                                                         ccAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGt                                                                                                                          >  2:154579/1‑250 (MQ=255)
                                                                                                                          cAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATACTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTa                                                                                                                         >  1:530941/1‑250 (MQ=255)
                                                                                                                           aGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGACTGACGACAACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGATCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGCCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTcacaacaggtaac                                                                                                                       >  2:523596/1‑239 (MQ=255)
                                                                                                                           aGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAac                                                                                                                       >  1:227865/1‑251 (MQ=255)
                                                                                                                                  tgaaatgaCTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAACGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTACCgatgac                                                                                                                 >  1:506848/1‑250 (MQ=255)
                                                                                                                                    aaatgaCTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGCG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGAtgacatc                                                                                                              >  1:224896/1‑251 (MQ=255)
                                                                                                                                            tgatctgaacgcacaaacccTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGCAGACACCGCCGAG‑ACACTAATTGCC‑CTGTGTCATCAGGCCAACAAC‑CCGGCCGGCAATACCGCCGCTACCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCCAACCACTGAAAGCGCAGGGACCCCCGGAAATACGTATTGCTACGGTACCCAACTTCCCACACGGCAACGAagccatcgaaat                                                                                                         >  2:269142/19‑246 (MQ=255)
                                                                                                                                            aacaagAAAGCAA‑GGAGCTTGGGTCCACTGAAATTGAGGGTCCTCCCCCCCGTGAATGCCAACGCCCCCGAAGAG‑AAATTAATTTCC‑CTTTGTCATCAGTCAAA‑ACCCCCGGCCGGCATTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATgacatcgaaatc                                                                                                        <  2:565572/244‑1 (MQ=255)
                                                                                                                                                      aaa‑gcaGCCTGAGTGCACTGAAATTGATGGACCTGACCACCATGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGAAGGGCACCCCGGAAATCCGTATTGCTACGGTAAACAACTTCACACACGGTACCGATGACATAGAAATCGCGCTGGCAGa                                                                                             >  2:443218/2‑250 (MQ=255)
                                                                                                                                                       aa‑gcaGCCTGGGTGCACTGAAATTGATGGACCTGACCCCCCTGAAGGACGACGACACCGACGAG‑AAAGTATTTGCT‑CTGTGTAATCAGGCCAA‑AACCCCGGCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                                                                           <  1:109799/251‑1 (MQ=255)
                                                                                                                                                       aa‑gcaGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATAcc                                                                                                                                                                                                                                      <  2:351546/112‑1 (MQ=255)
                                                                                                                                                       aa‑gcaGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATAcc                                                                                                                                                                                                                                      >  1:351546/1‑112 (MQ=255)
                                                                                                                                                                 gCGTGCACTGAACTTGATGGACCTGACAACCCTGAATGACGACGACACCGACGAG‑AACGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCAGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTCTCCCGATTGCTCGCAAACCACTGAACGAGCAGGGCACCCCGGCAATCCGTATTGCTACGGTAACCAACTTCCCAAACGGTAACGATGACATCGAAATCGCGCTGGCAGAACCCCGTGCgg                                                                                  >  2:37992/1‑251 (MQ=255)
                                                                                                                                                                      cacgaaaaTTGAGGGCCCTGACCCCCCTAAATGACGCCCACCCCTACGAA‑AAAGTAATGGCT‑CTGTTCCACCAGGCCAA‑AACCCCGGCCTGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTGCCCACACGGTAACTATGACAGCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                             <  2:82453/246‑1 (MQ=255)
                                                                                                                                                                      cACTGAACTTGATGGACCTGACCACCCTGACTGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                             >  2:65502/1‑251 (MQ=255)
                                                                                                                                                                           aaatggagggcccTGACCCCCCTGAATGACGACGCCACCAACGAG‑AAAGTAATTGCT‑TTGTGTCATCTGGCCAA‑CACCCCGGTCGGCAGTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGAGTTCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTATCCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTa                                                                        <  1:91159/240‑1 (MQ=255)
                                                                                                                                                                           aaaTTGATGGCCCTCACCACCCTGAATCACGACGCCACCGACGAA‑AAAGTATTTGCT‑CTGTGTCATCAGGCCAA‑ACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTa                                                                        <  2:76264/251‑1 (MQ=255)
                                                                                                                                                                            aaTTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGAT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTAc                                                                       >  2:243875/1‑251 (MQ=255)
                                                                                                                                                                                   ggccctgccccccctaaagacggacgacCCGGACGAA‑AAAGTAATTGCT‑CTGTGTCATCAGACAAA‑ACCCCCGGTCGACATTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGccg                                                                <  2:181774/229‑3 (MQ=255)
                                                                                                                                                                                          accaccCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACAGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGt                                                         >  1:561909/1‑251 (MQ=255)
                                                                                                                                                                                                    aggacgccgacACCGACGAA‑AAATTAATTGCT‑TTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACgtggtg                                               <  1:545218/249‑1 (MQ=255)
                                                                                                                                                                                                      gacgacgacACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGtt                                             >  1:57641/1‑251 (MQ=255)
                                                                                                                                                                                                            gacACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGt                                        >  1:327681/1‑250 (MQ=255)
                                                                                                                                                                                                            gacACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGt                                        >  1:182200/1‑250 (MQ=255)
                                                                                                                                                                                                                  gacgaG‑AAAGTCATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGACATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATAGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTAAGGTGCCGTTGAAGTTGACGTGGTGTTCCCGTACcgcgc                                 >  2:516144/1‑251 (MQ=255)
                                                                                                                                                                                                                    cgaG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCACTCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCgcgcgc                               >  2:574632/1‑251 (MQ=255)
                                                                                                                                                                                                                     gag‑aAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTACCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCt                              >  2:463881/1‑251 (MQ=255)
                                                                                                                                                                                                                       aagaaaggaattgctcTGTGACATCAGGCCAA‑AACCCCGGTCGGCAATACCGACGCTATCTGTATCTATCCTCGCTATATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCACACACGGTAACGATGACATCGAAATCGCGATGGCAGAAACCCGCGAGACAATCGCCTACGATGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCt                              >  2:342111/16‑251 (MQ=255)
                                                                                                                                                                                                                               aattgctcTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGAACCCCGGCAATCCGTATTGCTACGGTAACCAACTTCCCACACGGAAACGATGACATCGCAATCGCGCTGGCAGACACCCGGGCGGCAATCGCCTACGGTGCCGATGAGGTCGACGTGGTGTTCCCGTACCGCGCGCTGATGGCgg                      >  2:514123/8‑251 (MQ=255)
                                                                                                                                                                                                                                 ttgctcTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGt                    >  1:581302/6‑251 (MQ=255)
                                                                                                                                                                                                                                     gctgttccatcaGCCCAA‑ACCCCCGGGCGACAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCGGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACg                <  2:17534/244‑1 (MQ=255)
                                                                                                                                                                                                                                        gtgctatcaGCCCAA‑ACCCCCGGCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTAGCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGc             <  2:143696/246‑1 (MQ=255)
                                                                                                                                                                                                                                          ttcatcaGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAACGAGCAGGGCACCCCGGAAATCCGTCTTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCTTACCGCGCGCTGATGGCggggccatcgcac           >  2:214177/2‑241 (MQ=255)
                                                                                                                                                                                                                                             atcaGGCCAA‑ACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggtt        <  1:283307/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tcaGGCCAA‑ACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGAGGAAGTGGACGTGGTGTTCCCGGACCGCGCGCTGATGGCGGGTAACGAGCAggtt        <  2:444435/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGCGCAGGTTGGtttt  >  1:440625/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   
GCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTGATCGCC‑CTGTGTCATCAGGCCAA‑AACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT  >  NC_000913/4617190‑4617679

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