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,7720CA96.2% 154.4 / NA 53G155G (GGC→GGAppnNnucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (28/23);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATCATGCCAGA  >  NC_000913/2926530‑2927017
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
gAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTc                                                                                                                                                                                                                                               >  1:181146/1‑251 (MQ=255)
      aTTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATATCGACAACTCCGTCCCCATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCACATATGGTGGTCTGCTGGGGCGGTCACTCAATTa                                                                                                                                                                                                                                          >  2:199846/1‑250 (MQ=255)
          atctttcggcagagGCTTTTGTCGAAGGGCGAATTATTCGCGCATTGCAGGCCAACTTGTTCGCGGTCCCGCGTGACATTCTCTCCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATATCGACAACTCAGTCCACATCTCTAACCTGGTCTTTTCCATCTTGCGTAAAGCTCGCGCGCTGCATGTGGGTGAGGCGCAAACTAAGGTGGTCTGCTGGGGAGGTCACTCacttcacgc                                                                                                                                                                                                                                      >  2:558943/11‑242 (MQ=255)
            cgtccGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTCCGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGaaaa                                                                                                                                                                                                                                   <  2:34423/248‑1 (MQ=255)
              tccGGAAGAGGCTTTTGTCGATGGGCGAATGATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCAGCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAAc                                                                                                                                                                                                                                  <  2:184040/249‑1 (MQ=255)
                     gagGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTAttt                                                                                                                                                                                                                          >  1:307369/1‑251 (MQ=255)
                             tGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCACATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGccc                                                                                                                                                                                                                  >  2:35222/1‑251 (MQ=255)
                             tGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTGCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGcc                                                                                                                                                                                                                   <  1:37691/250‑1 (MQ=255)
                                cGATGGGCGAATTATTCGCGGTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtc                                                                                                                                                                                                               <  2:407066/251‑1 (MQ=255)
                                 gATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcg                                                                                                                                                                                                              >  2:565138/1‑251 (MQ=255)
                                    gggCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcgcgt                                                                                                                                                                                                           >  2:542811/1‑251 (MQ=255)
                                        gAATTATGCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTACGTCCACATCACTAACCTTGCCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGACGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAACACGAGTATTTGTATGCCcgtcccttcggg                                                                                                                                                                                                       >  2:571680/1‑243 (MQ=255)
                                          attattCGCGCTTTGCAGGCCAACTTGTTAGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGaa                                                                                                                                                                                                      >  1:473670/1‑250 (MQ=255)
                                            tattCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGaacc                                                                                                                                                                                                    >  2:223453/1‑248 (MQ=255)
                                            tattCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAAcc                                                                                                                                                                                                   >  1:560465/1‑251 (MQ=255)
                                             attCTCGCTTTGCAGGCCAACCTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCa                                                                                                                                                                                                  <  2:307369/251‑1 (MQ=255)
                                                   gttttgatggccaatgtGTTCGCGGGCCTGGGTGACATTCTCTTCGTTGACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCCGAAGCACGAAAGCTCTGGCCACATCACTAACCTGTTCTGTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCGCTCAATTAACGAAAGCGATGATGTGTATGCCCGTCGTGTCGGAAGCCAGCTggg                                                                                                                                                                                            <  2:371713/235‑1 (MQ=255)
                                                   gcTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTggg                                                                                                                                                                                            >  1:284979/1‑251 (MQ=255)
                                                         cAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAACGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCg                                                                                                                                                                                      >  2:483458/1‑251 (MQ=255)
                                                                      ccgcgGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGGTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTTTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGATTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAatat                                                                                                                                                                          <  2:247891/249‑1 (MQ=255)
                                                                             ccTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCCCCg                                                                                                                                                                  >  2:114044/1‑251 (MQ=255)
                                                                             ccTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCg                                                                                                                                                                  <  1:483458/251‑1 (MQ=255)
                                                                                              tCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTTAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGc                                                                                                                                                  <  2:287753/250‑1 (MQ=255)
                                                                                                     cGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGa                                                                                                                                            >  2:279834/1‑249 (MQ=255)
                                                                                                          aaaTCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCg                                                                                                                                     >  1:419072/1‑251 (MQ=255)
                                                                                                               cATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGa                                                                                                                                 >  1:55782/1‑250 (MQ=255)
                                                                                                                aTAACACCGTTCGTTTTCCCAACCTGAATCTCGCCAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGCTGGAAGCGCCGCTGaa                                                                                                                                >  2:238263/1‑250 (MQ=255)
                                                                                                                      ccgttctggttcccaaccgttaTCTCGGCAATTCCGTCCACAGCATTAACCTGGTCTTTTCCAGCTTGCGGAACGCTCGCGCGCTGCATGTGGGTGAATCGCTAAATATGGTGGTCTGCTTGGGAGGTCACTCAATTAACGAAAACGAGGAGTTGTATGCCCGTCGTGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGTAAGCGCCGAGGAAAGGtgct                                                                                                                         <  2:395910/230‑1 (MQ=255)
                                                                                                                                      ccTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACgcgc                                                                                                         >  1:412580/1‑251 (MQ=255)
                                                                                                                                               cGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCagcagc                                                                                                   <  1:565138/248‑1 (MQ=255)
                                                                                                                                                   aaCTCCGTCCCCATCACTAACCTGGCCCTGTCCAGCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGAGCCAAATATGGTGGTCTGCTGGGGAGGTCACGCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCGTGGCCTGCGTGAGCTGAATATCTGCACCTGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCTGACACGCGCAGCAGAGTTACaa                                                                                            <  2:459316/251‑1 (MQ=255)
                                                                                                                                                        cGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCCATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGGGCGATGGAAGCGCCGATGAAAGGGGCTGCGGTCGGGCACGCGCAGCAGCGTTACAAAGAc                                                                                        >  2:279463/1‑250 (MQ=255)
                                                                                                                                                           ccACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGt                                                                                     <  1:122049/250‑1 (MQ=255)
                                                                                                                                                            cacaTCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTAGGACACGCGCAGCAGCGTTACAAAGACAGTCg                                                                                   >  2:170758/1‑251 (MQ=255)
                                                                                                                                                             acaTCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGt                                                                                  >  1:151669/1‑251 (MQ=255)
                                                                                                                                                              caTCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTGGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCGGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGtt                                                                                 <  2:55782/251‑1 (MQ=255)
                                                                                                                                                              caTCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGCATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGtt                                                                                 >  2:528088/1‑251 (MQ=255)
                                                                                                                                                                          ggTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATg                                                                      <  1:13202/250‑1 (MQ=255)
                                                                                                                                                                             cTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATTTGGGTGACTCGCCAAATATGGGGGTCTGCTGGGGAGGTCACTCAATTAACGAAACCGAGTATTTGGCTGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCCGCACCGGCTGTGGTCCGGGAGCGATGGGAGCGCCGCTGACAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACCACGACAGGCGTTTTCTTGGtctgta                                                                    >  2:105636/1‑244 (MQ=255)
                                                                                                                                                                                 tCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGc                                                               <  1:35222/250‑1 (MQ=255)
                                                                                                                                                                                  ccATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCATCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGcc                                                              <  1:209319/250‑1 (MQ=255)
                                                                                                                                                                                        tGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCtgtg                                                                                                                                                            <  1:310281/150‑1 (MQ=255)
                                                                                                                                                                                        tGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCtgtg                                                                                                                                                            >  2:310281/1‑150 (MQ=255)
                                                                                                                                                                                           gTAACGCTCGCGCGTTGCATGTGGGTGAAGCGACAAATAGGTGGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATAGCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATc                                                    <  2:419072/251‑1 (MQ=255)
                                                                                                                                                                                             aaCGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTc                                                           <  1:126408/242‑1 (MQ=255)
                                                                                                                                                                                             aaCGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAAACAGTCGTTTTATTGGTATGACAGAGCCGTcagatcggg                                                   >  2:126408/1‑242 (MQ=255)
                                                                                                                                                                                                    gcgcgcTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGaa                                           <  1:376045/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTc                         >  1:549037/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ggtggtCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACa                                                                                             >  1:4933/1‑170 (MQ=255)
                                                                                                                                                                                                                                   ggtggtCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACa                                                                                             <  2:4933/170‑1 (MQ=255)
                                                                                                                                                                                                                                      ggtCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGAtcatca         <  1:114044/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tgGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGCATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCACCGAATTGATCATCATGcccg   >  1:98091/1‑248 (MQ=255)
                                                                                                                                                                                                                                               ggggAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGGGCTGCGGGCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCCTCATGCCAgt  >  2:155166/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
GAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATCATGCCAGA  >  NC_000913/2926530‑2927017

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: