Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,412 C→T I30I (ATC→ATT deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,4120CT98.4% 196.2 / NA 64I30I (ATC→ATTdeoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (48/14);  minor base G (0/1);  total (48/16)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCT  >  NC_000913/4617181‑4617651
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
aaCCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGAAATTGCTCTGTGTCATCAGGcc                                                                                                                                                                                                                                >  2:545218/1‑250 (MQ=255)
  ccAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCaa                                                                                                                                                                                                                              >  1:76264/1‑250 (MQ=255)
   cAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGCGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCaaa                                                                                                                                                                                                                            >  2:283307/1‑250 (MQ=255)
         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:316608/1‑251 (MQ=255)
                                                                                                                               aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                                                >  1:120103/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)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
AACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCT  >  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: