Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,926,772 C→A G155G (GGC→GGA ppnN → nucleotide 5'‑monophosphate nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,7720CA98.4% 185.2 / NA 61G155G (GGC→GGAppnNnucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (32/28);  minor base G (1/0);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATCATGCC  >  NC_000913/2926529‑2927014
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
ggAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCAc                                                                                                                                                                                                                                              <  2:127628/250‑1 (MQ=255)
    cTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGGCATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCAGTCAACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAAt                                                                                                                                                                                                                                         >  1:424852/1‑251 (MQ=255)
           aTCCCCCGGAAGATGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACAGCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGaa                                                                                                                                                                                                                                  <  2:361171/251‑1 (MQ=255)
               cccGGAAGAGGCTTTTGTCGATGGGTGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCCAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACg                                                                                                                                                                                                                              <  1:897/251‑1 (MQ=255)
                      gagGCTTTTGTCGCTGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTATTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTAttt                                                                                                                                                                                                                       >  1:435814/1‑251 (MQ=255)
                       agGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTg                                                                                                                                                                                                                      >  1:240285/1‑251 (MQ=255)
                                   aTGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAACCGAGTATTTGTATGCCcgtcgc                                                                                                                                                                                                          >  2:225601/1‑251 (MQ=255)
                                    tGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcgcg                                                                                                                                                                                                         >  2:75680/1‑251 (MQ=255)
                                                         gCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTg                                                                                                                                                                                     <  1:406170/250‑1 (MQ=255)
                                                               cAACTTGTTCGCGGTCCTGCGTGACAGTCTCTTCGTTGACGGGCAAATCCATAACTCCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTATCCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGATGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGc                                                                                                                                                                              <  1:383158/251‑1 (MQ=255)
                                                                          cgGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGGCCACATCACTAACCTGGTCTTTTCCATCGTGCGTAACGCTCGCGCGCTGCAGGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGATGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGc                                                                                                                                                                   <  2:391793/251‑1 (MQ=255)
                                                                                  cGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCtg                                                                                                                                                           >  1:143549/1‑251 (MQ=255)
                                                                                   gTGACATTCTGTTCGTTTACAGGCCAATCCATAACACCGTTCGTTTTCCCCACCTGCCTCTCGACAACTCCGTCCACCTCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTATCGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCAGGGCCTGCGTGAGATGACTATCTGCACCGGCtgg                                                                                                                                                          >  2:437001/1‑250 (MQ=255)
                                                                                    tGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCtgtg                                                                                                                                                         >  1:169606/1‑251 (MQ=255)
                                                                                       cATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTc                                                                                                                                                      <  2:143549/251‑1 (MQ=255)
                                                                                           ctctTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCggg                                                                                                                                                  >  2:261778/1‑251 (MQ=255)
                                                                                                cGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCg                                                                                                                                              <  1:258577/250‑1 (MQ=255)
                                                                                                cGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGACACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGa                                                                                                                                             <  1:252169/251‑1 (MQ=255)
                                                                                                 gTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGAt                                                                                                                                            <  2:65841/251‑1 (MQ=255)
                                                                                                  tttACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGTCAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATg                                                                                                                                           >  2:121153/1‑251 (MQ=255)
                                                                                                     aCGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTACCCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCCAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGa                                                                                                                                         >  2:73444/1‑250 (MQ=255)
                                                                                                      cGGGCAAATCCATAACCCCGTTCGTGTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCAGCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACGCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAg                                                                                                                                       <  1:327735/251‑1 (MQ=255)
                                                                                                      agggCAAATCCATAACACAGTTCGTTTTCCCCACCTGAATCTCGACACCTCAGTCCACATCACTAACCTGGTATTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGGGGTCAGTCAATTAACGAAAACGAGAATTTGAATGCCAGACGCGTAGGAAACGAGCTGGGCCTGCGTGTGCTGGATATCTGCACCGGCTGGGGTCCGGCCGCGATGGacg                                                                                                                                       >  2:169998/2‑249 (MQ=255)
                                                                                                            aaTCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGc                                                                                                                                                                   <  1:318008/217‑1 (MQ=255)
                                                                                                            aaTCCATAACACCGTTCGTTTTCCCAACCTGAATCTAGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGGGAGCTGAATATCTGc                                                                                                                                                                   >  2:318008/1‑217 (MQ=255)
                                                                                                                     cacCGTTCGTTTTCCCAACCTGAATCACGACAACTACGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCACATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAACCCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGCCTGTGGTCCTGGAGCGATGGCAGCGCCGATGAACGGtg                                                                                                                        >  2:139094/1‑251 (MQ=255)
                                                                                                                       ccGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGtgct                                                                                                                      <  1:75680/251‑1 (MQ=255)
                                                                                                                         gttcgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGtgctgc                                                                                                                    >  2:120961/1‑251 (MQ=255)
                                                                                                                         gttcgttTTCCAAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGATCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGtgctgc                                                                                                                    >  1:110002/1‑251 (MQ=255)
                                                                                                                            cgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAACCCAGCTGGGCCTGCGTGAGCt                                                                                                                                                                             <  2:347427/191‑1 (MQ=255)
                                                                                                                            cgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCt                                                                                                                                                                             >  1:347427/1‑191 (MQ=255)
                                                                                                                                 tCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTATCCATCTTGCGTAACGCTCGCGCGCTGAATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAACGGTGCTGCGGTCGGac                                                                                                            >  2:78656/1‑251 (MQ=255)
                                                                                                                                  cccAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGac                                                                                                            <  1:261733/250‑1 (MQ=255)
                                                                                                                                      aCCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGaa                                                                                                                                        >  1:457760/1‑218 (MQ=255)
                                                                                                                                      aCCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGaa                                                                                                                                        <  2:457760/218‑1 (MQ=255)
                                                                                                                                        cTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACgcgcc                                                                                                     >  2:86662/1‑250 (MQ=255)
                                                                                                                                            aTCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcg                                                                                                                                                                                                           >  1:281220/1‑145 (MQ=255)
                                                                                                                                            aTCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcg                                                                                                                                                                                                           <  2:281220/145‑1 (MQ=255)
                                                                                                                                            aTCTCGACAACTACGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAACCGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCgcagc                                                                                                   >  1:391714/1‑249 (MQ=255)
                                                                                                                                             cctcGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTGGCGTAACGCTCGCCCGCCGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACTCGCagcagc                                                                                                <  2:82710/250‑1 (MQ=255)
                                                                                                                                                            ccACAGCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCTCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCTTTACAAAGACAGt                                                                                  <  2:398581/250‑1 (MQ=255)
                                                                                                                                                                    ccgaccttgttgtttacATCGTGCGTAACGCTTGCACGCTGCATGTGAGTGAAGCGCCAAATATGGTGGTCTTCTGGGGAGGTCACGCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCGCCTGCTGTGGTCCGGGAGCGATGGAAGCTCCGATGAAAGGTGTTGCGGTCGGACACGCACTGCAGCGTTACAAAGACATTCGTTTTAtt                                                                         <  2:296198/235‑1 (MQ=255)
                                                                                                                                                                     tAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATTTGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCTGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTg                                                                        >  2:386340/1‑251 (MQ=255)
                                                                                                                                                                      aaCCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTgg                                                                       >  2:257837/1‑251 (MQ=255)
                                                                                                                                                                       aCCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGt                                                                      >  1:149719/1‑251 (MQ=255)
                                                                                                                                                                                  tCCATCTTGAGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCACTTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACACAGACAGTCGTTTTATTGGTATGACAGAGcc                                                           >  2:419936/1‑251 (MQ=255)
                                                                                                                                                                                       cTTGCGTAACGCTCGCGCGCTGCATGTGGGTGCAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAACCGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGACTATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGAGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTCTTGGCATGACAGAGCCGTCg                                                       >  2:412213/1‑250 (MQ=255)
                                                                                                                                                                                            gTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATc                                                 <  1:181140/251‑1 (MQ=255)
                                                                                                                                                                                              aaCGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGATCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATcg                                                <  2:110002/250‑1 (MQ=255)
                                                                                                                                                                                                    cgcgcgCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGa                                         <  1:225601/251‑1 (MQ=255)
                                                                                                                                                                                                        cgcTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCg                                     >  1:347106/1‑251 (MQ=255)
                                                                                                                                                                                                                  gggTGAAGCGCCAAATATGGAGGTCTGCTGGGGAGGTCACTCAATTAACGCAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGAGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGc                           >  2:44349/1‑251 (MQ=255)
                                                                                                                                                                                                                   ggTGAAGCGCCAAATATGGTGGTCTGCTGGGTAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGGCGTGTTATTGGTATTGCAGAGCCGTCTATTATCGCCGCTGAACCTCCTAACCCGCt                          <  2:391714/251‑1 (MQ=255)
                                                                                                                                                                                                                          cgcCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAAc                   >  2:185590/1‑251 (MQ=255)
                                                                                                                                                                                                                           gcCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGTTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACg                  <  1:437707/251‑1 (MQ=255)
                                                                                                                                                                                                                           gcCAAATATGGTGGTCTGCTGGGGAGGTCACGCAATTAACGAAAACGAGTAGTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCGGAATATCTGCCCCGGCTGTGGTCCGGTAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACg                  <  1:386340/251‑1 (MQ=255)
                                                                                                                                                                                                                            ccAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGa                 <  1:84594/251‑1 (MQ=255)
                                                                                                                                                                                                                            ccAAATATGGTGGGCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGAGCAGAAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGACGACGATTATCGCCGCTGACCCGGCTAACCCGCTGGTCAACg                  >  2:43996/1‑250 (MQ=255)
                                                                                                                                                                                                                                  atGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGat           >  1:91514/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gtCTGCTGGGTAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAGCCAGCTGGGCCTGCGTGAGTTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCATCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGAGcatcat     <  1:266729/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tgGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATCATGcc  >  1:378978/1‑248 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
GGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATCATGCC  >  NC_000913/2926529‑2927014

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