Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,409 G→A V29V (GTG→GTA)  deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,4090GA98.4% 193.1 / ‑6.0 64V29V (GTG→GTA) deoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (47/16);  minor base C (1/0);  total (48/16)
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.20e-01

AACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTGATCGCC‑CTGTGTCATCAGGCCAA‑AACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCT  >  NC_000913/4617181‑4617651
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aaCCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGAAATTGCT‑CTGTGTCATCAGGcc                                                                                                                                                                                                                                >  2:545218/1‑250 (MQ=255)
  ccAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCaa                                                                                                                                                                                                                              >  1:76264/1‑250 (MQ=255)
   cAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGCGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCaa‑a                                                                                                                                                                                                                            >  2:283307/1‑250 (MQ=255)
         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:316608/1‑251 (MQ=255)
                                                                                                                               aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                                                >  1:120103/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)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
AACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTG‑AATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAA‑GCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAG‑AAAGTGATCGCC‑CTGTGTCATCAGGCCAA‑AACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCT  >  NC_000913/4617181‑4617651

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