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

GCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT  >  NC_000913/4617190‑4617679
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   
gCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATTGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGAAATTGCTCTGTGTCATCAGGCCAAAACCCCgg                                                                                                                                                                                                                                                 >  1:79280/1‑251 (MQ=255)
           aTGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATAcc                                                                                                                                                                                                                                      <  1:155285/251‑1 (MQ=255)
                     gTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGCGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCCGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGACTGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCt                                                                                                                                                                                                                            >  1:596364/1‑251 (MQ=255)
                                 ccTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTctcctc                                                                                                                                                                                                                >  1:514669/1‑251 (MQ=255)
                                     aattgtgatgtagtcctaagggtgtTGGGGAGTAGATGGGAGAATACGAAGAACCTCGCAAGTGTATTTTTAGTGGCGACAAGCTAGGAAAAGGAAATAACTGATCTGAAAGCAAGCACCCGGCTTGCACGGAAATTGTTGGAGCTTACCTCCCTGAAGGTCGACGACACCGTCAGTAAGGTAATTCCTTTGTGTCATCGGGCGAAAACCCCTGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCt                                                                                                                                                                                                             <  1:409477/229‑1 (MQ=255)
                                        tgtgATGTATATCGAAGTGTGTTGCGAAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATATGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTAt                                                                                                                                                                                                         >  2:536895/1‑251 (MQ=255)
                                           gATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATccc                                                                                                                                                                                                      >  1:143696/1‑251 (MQ=255)
                                           gATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGAAAGGTGACTTTTATTGGCGACAAGACAGGCGAATGAAATGCCTGATCTGAAAGCAAGCAGCCTACGTGCACTGAAATTGCTGGACCTGACCACCCTGAATGACGACGACACCGACGAGAACGTAATTGCTCTGTGTCATCAGGCAAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATccc                                                                                                                                                                                                      >  1:595581/1‑251 (MQ=255)
                                               tatATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGAtt                                                                                                                                                                                                  >  1:289964/3‑251 (MQ=255)
                                                     gAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACACACTGGCCCGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGAAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGc                                                                                                                                                                                            >  2:43996/1‑251 (MQ=255)
                                                           ggtgggggAGTAAATGTTAAAATACAAAAAAACTCGCAAGGTTAATTTTATTGGCGACAAGCCAGGAGAAAGAAATGACTGATCTGAAAGCAAGCAGTCTGCGTGCACTGAAATGAATGGCCCTGCCCCCCCTGAATGACGCCGACACCGACGAAAAAGTATTTGCTCTGTTTCATCAGGCCAAACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCTCTTTATCCCGATTGCTCGCAAAAca                                                                                                                                                                                      <  2:327131/245‑2 (MQ=255)
                                                                  gAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGaa                                                                                                                                                                                 >  1:58940/1‑249 (MQ=255)
                                                                    gTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAag                                                                                                                                                                               >  1:444435/1‑249 (MQ=255)
                                                                       gATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGc                                                                                                                                                                            >  1:329654/1‑249 (MQ=255)
                                                                          gTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATAcc                                                                                                                                                                                                                                      <  2:190875/188‑1 (MQ=255)
                                                                           ttAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACAGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGc                                                                                                                                                                       >  1:1075/1‑250 (MQ=255)
                                                                              gAATAATAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTCTCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCAccc                                                                                                                                                                   >  2:614593/1‑251 (MQ=255)
                                                                               aaataataaaaaactgaaaagggaaaTTTTATTGGCAACAAGCCAGAAGAATGAAATGTCTGATCGAAAACAAACCAGCCGGCGTGCCCGGAAATTGCTGGACCTGACCACCCTGAATGACCACACCCCCGCAGAGAAAGTAATTGCTCTGTTTAAACAGCCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCAccc                                                                                                                                                                   <  2:888/227‑1 (MQ=255)
                                                                               aaTACTAACAAACTCGCAAGGTGAATTTTATTGGCGACCAGCCAGGAGCATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCACACCCCCGGTCGGCACTACCGCCGCTATCTGTCTCTATCCTCGCTTTATCCCGATTGCTCGCAACACACTGAAAGCTCAGGTCAcccc                                                                                                                                                                  >  2:55042/1‑251 (MQ=255)
                                                                               aaTACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCAcccc                                                                                                                                                                  >  1:434162/1‑251 (MQ=255)
                                                                                     aaCAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAAt                                                                                                                                                            >  1:82453/1‑251 (MQ=255)
                                                                                        aaaCTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCg                                                                                                                                                         >  1:16728/1‑251 (MQ=255)
                                                                                          aacccgaaagggaaaTTTTATTGGGAACAAGCCAGAAAAATGAAATGATGGAACGGAAAGCAAGCAGCCTGCGTGCACTAAAATTGAGGGCCCTGACCCCCCTAAAGGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCGGGCAAAACCCCCGGGCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGt                                                                                                                                                        <  2:257466/238‑1 (MQ=255)
                                                                                          aCTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGCAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAACTCCGt                                                                                                                                                        >  2:523174/1‑250 (MQ=255)
                                                                                           cTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTAt                                                                                                                                                      >  1:439727/1‑251 (MQ=255)
                                                                                           atCGCACGGTGAGTTTTATTGGCGACCAGCCAGGCGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTAt                                                                                                                                                      >  2:125960/2‑251 (MQ=255)
                                                                                               cAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCATGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCt                                                                                                                                                  >  1:222799/1‑251 (MQ=255)
                                                                                                aaGGTGAATTTTATGGGCGACAAGCAAGGAGAATGAAATGAATGATCTGAAAGCAAGACGCCGGCGGGCGCTGAAATGGATGGCCCTGACCCCCCTGAAGGACAACGACACCGACGAGAAAGTAATGCCTCGGGGTCATCAGGCCAAACCCCGGGCCGGAAATACCGCCGCTATCTGTATCGATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTa                                                                                                                                                 <  1:614593/251‑1 (MQ=255)
                                                                                                  ggTGAATTTTATTGGCGACAAGCCAGGAGATTGAAAGGACGGATCGAAAAGAAAGCAGCCTGCGTGCACTGAATGTGATGGACCTGACCCCCCTGAATGACGACGACCCCGACGAGAAAGTAATTGCTCTGTGTCTTAAGGCAAAAACCCCGTTCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACg                                                                                                                                               <  1:332094/251‑1 (MQ=255)
                                                                                                        aTTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAc                                                                                                                                          >  1:494934/1‑251 (MQ=255)
                                                                                                                    aCAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCcac                                                                                                                               >  1:108032/1‑250 (MQ=255)
                                                                                                                     cAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCACAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCacac                                                                                                                             >  2:417511/1‑251 (MQ=255)
                                                                                                                      aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                                                                            >  1:120103/1‑251 (MQ=255)
                                                                                                                      aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                                                                            >  1:316608/1‑251 (MQ=255)
                                                                                                                         ccAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGt                                                                                                                          >  2:154579/1‑250 (MQ=255)
                                                                                                                          cAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATACTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTa                                                                                                                         >  1:530941/1‑250 (MQ=255)
                                                                                                                           aGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGACTGACGACAACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGATCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGCCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTcacaacaggtaac                                                                                                                       >  2:523596/1‑239 (MQ=255)
                                                                                                                           aGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAac                                                                                                                       >  1:227865/1‑251 (MQ=255)
                                                                                                                                  tgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAACGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTACCgatgac                                                                                                                 >  1:506848/1‑250 (MQ=255)
                                                                                                                                    aaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGCGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGAtgacatc                                                                                                              >  1:224896/1‑251 (MQ=255)
                                                                                                                                            tgatctgaacgcacaaacccTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGCAGACACCGCCGAGACACTAATTGCCCTGTGTCATCAGGCCAACAACCCGGCCGGCAATACCGCCGCTACCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCCAACCACTGAAAGCGCAGGGACCCCCGGAAATACGTATTGCTACGGTACCCAACTTCCCACACGGCAACGAagccatcgaaat                                                                                                         >  2:269142/19‑246 (MQ=255)
                                                                                                                                            aacaagAAAGCAAGGAGCTTGGGTCCACTGAAATTGAGGGTCCTCCCCCCCGTGAATGCCAACGCCCCCGAAGAGAAATTAATTTCCCTTTGTCATCAGTCAAAACCCCCGGCCGGCATTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATgacatcgaaatc                                                                                                        <  2:565572/244‑1 (MQ=255)
                                                                                                                                                      aaagcaGCCTGAGTGCACTGAAATTGATGGACCTGACCACCATGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGAAGGGCACCCCGGAAATCCGTATTGCTACGGTAAACAACTTCACACACGGTACCGATGACATAGAAATCGCGCTGGCAGa                                                                                             >  2:443218/2‑250 (MQ=255)
                                                                                                                                                       aagcaGCCTGGGTGCACTGAAATTGATGGACCTGACCCCCCTGAAGGACGACGACACCGACGAGAAAGTATTTGCTCTGTGTAATCAGGCCAAAACCCCGGCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                                                                           <  1:109799/251‑1 (MQ=255)
                                                                                                                                                       aagcaGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATAcc                                                                                                                                                                                                                                      <  2:351546/112‑1 (MQ=255)
                                                                                                                                                       aagcaGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATAcc                                                                                                                                                                                                                                      >  1:351546/1‑112 (MQ=255)
                                                                                                                                                                 gCGTGCACTGAACTTGATGGACCTGACAACCCTGAATGACGACGACACCGACGAGAACGTAATTGCTCTGTGTCATCAGGCCAAAACCCAGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTCTCCCGATTGCTCGCAAACCACTGAACGAGCAGGGCACCCCGGCAATCCGTATTGCTACGGTAACCAACTTCCCAAACGGTAACGATGACATCGAAATCGCGCTGGCAGAACCCCGTGCgg                                                                                  >  2:37992/1‑251 (MQ=255)
                                                                                                                                                                      cacgaaaaTTGAGGGCCCTGACCCCCCTAAATGACGCCCACCCCTACGAAAAAGTAATGGCTCTGTTCCACCAGGCCAAAACCCCGGCCTGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTGCCCACACGGTAACTATGACAGCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                             <  2:82453/246‑1 (MQ=255)
                                                                                                                                                                      cACTGAACTTGATGGACCTGACCACCCTGACTGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                             >  2:65502/1‑251 (MQ=255)
                                                                                                                                                                           aaatggagggcccTGACCCCCCTGAATGACGACGCCACCAACGAGAAAGTAATTGCTTTGTGTCATCTGGCCAACACCCCGGTCGGCAGTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGAGTTCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTATCCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTa                                                                        <  1:91159/240‑1 (MQ=255)
                                                                                                                                                                           aaaTTGATGGCCCTCACCACCCTGAATCACGACGCCACCGACGAAAAAGTATTTGCTCTGTGTCATCAGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTa                                                                        <  2:76264/251‑1 (MQ=255)
                                                                                                                                                                            aaTTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGATCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTAc                                                                       >  2:243875/1‑251 (MQ=255)
                                                                                                                                                                                   ggccctgccccccctaaagacggacgacCCGGACGAAAAAGTAATTGCTCTGTGTCATCAGACAAAACCCCCGGTCGACATTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGccg                                                                <  2:181774/229‑3 (MQ=255)
                                                                                                                                                                                          accaccCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACAGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGt                                                         >  1:561909/1‑251 (MQ=255)
                                                                                                                                                                                                    aggacgccgacACCGACGAAAAATTAATTGCTTTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACgtggtg                                               <  1:545218/249‑1 (MQ=255)
                                                                                                                                                                                                      gacgacgacACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGtt                                             >  1:57641/1‑251 (MQ=255)
                                                                                                                                                                                                            gacACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGt                                        >  1:327681/1‑250 (MQ=255)
                                                                                                                                                                                                            gacACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGt                                        >  1:182200/1‑250 (MQ=255)
                                                                                                                                                                                                                  gacgaGAAAGTCATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGACATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATAGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTAAGGTGCCGTTGAAGTTGACGTGGTGTTCCCGTACcgcgc                                 >  2:516144/1‑251 (MQ=255)
                                                                                                                                                                                                                    cgaGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCACTCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCgcgcgc                               >  2:574632/1‑251 (MQ=255)
                                                                                                                                                                                                                     gagaAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTACCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCt                              >  2:463881/1‑251 (MQ=255)
                                                                                                                                                                                                                       aagaaaggaattgctcTGTGACATCAGGCCAAAACCCCGGTCGGCAATACCGACGCTATCTGTATCTATCCTCGCTATATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCACACACGGTAACGATGACATCGAAATCGCGATGGCAGAAACCCGCGAGACAATCGCCTACGATGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCt                              >  2:342111/16‑251 (MQ=255)
                                                                                                                                                                                                                               aattgctcTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGAACCCCGGCAATCCGTATTGCTACGGTAACCAACTTCCCACACGGAAACGATGACATCGCAATCGCGCTGGCAGACACCCGGGCGGCAATCGCCTACGGTGCCGATGAGGTCGACGTGGTGTTCCCGTACCGCGCGCTGATGGCgg                      >  2:514123/8‑251 (MQ=255)
                                                                                                                                                                                                                                 ttgctcTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGt                    >  1:581302/6‑251 (MQ=255)
                                                                                                                                                                                                                                     gctgttccatcaGCCCAAACCCCCGGGCGACAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCGGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACg                <  2:17534/244‑1 (MQ=255)
                                                                                                                                                                                                                                        gtgctatcaGCCCAAACCCCCGGCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTAGCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGc             <  2:143696/246‑1 (MQ=255)
                                                                                                                                                                                                                                          ttcatcaGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAACGAGCAGGGCACCCCGGAAATCCGTCTTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCTTACCGCGCGCTGATGGCggggccatcgcac           >  2:214177/2‑241 (MQ=255)
                                                                                                                                                                                                                                             atcaGGCCAAACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggtt        <  1:283307/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tcaGGCCAAACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGAGGAAGTGGACGTGGTGTTCCCGGACCGCGCGCTGATGGCGGGTAACGAGCAggtt        <  2:444435/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGCGCAGGTTGGtttt  >  1:440625/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   
GCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT  >  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: