Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,785,913 C→T T218T (ACC→ACT gpmM → 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,785,9130CT98.2% 192.5 / NA 56T218T (ACC→ACTgpmM2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/29);  minor base G (0/1);  total (26/30)
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

ACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGC  >  NC_000913/3785667‑3786155
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
aCTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCAGTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATActgccg                                                                                                                                                                                                                                                <  1:359389/251‑1 (MQ=255)
   ggCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGCTACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCAGTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTg                                                                                                                                                                                                                                             <  1:205320/251‑1 (MQ=255)
   ggCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGtt                                                                                                                                                                                                                                              <  1:139444/250‑1 (MQ=255)
   ggCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTAGACGGACGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAACGGCCGCGTAGCATCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGAAGGGTGAGTTCCAGGCCGCTACTGCCGTTg                                                                                                                                                                                                                                             >  1:253967/1‑251 (MQ=255)
          gAACGCGGCGCAGAAAAAATCTACCTGCACGCATCCCTTGACGGTCGCGACAACCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACTCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGTGTGAGTTCCAGGCCGATACTGCCGTTGCTGGttt                                                                                                                                                                                                                                      <  1:107959/251‑1 (MQ=255)
           aaCGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTg                                                                                                                                                                                                                                     >  1:67226/1‑251 (MQ=255)
                     aGAAAAAAGCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTGCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCtt                                                                                                                                                                                                                           <  2:278536/251‑1 (MQ=255)
                      gAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGACCAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACAGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCtt                                                                                                                                                                                                                           >  2:66010/1‑250 (MQ=255)
                       aaaaaaTCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTAGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTAt                                                                                                                                                                                                                         >  2:95167/1‑251 (MQ=255)
                         aaaaTCTACCTGCACGCATTCATTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGAGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGAAGGCAGCTTATGc                                                                                                                                                                                                                       >  2:383157/1‑251 (MQ=255)
                                aCCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTcgct                                                                                                                                                                                                                  <  2:275578/249‑2 (MQ=255)
                                     cACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGaaa                                                                                                                                                                                                            >  2:130203/1‑250 (MQ=255)
                                         cATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCGCGCGACGaaaacgac                                                                                                                                                                                                       <  1:122332/251‑5 (MQ=255)
                                         cATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGACCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGaaaacgac                                                                                                                                                                                                       <  1:321877/251‑5 (MQ=255)
                                          aTTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGACATTTGCTGCGCAGGGAAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGaaaacgacg                                                                                                                                                                                                      >  1:297514/1‑246 (MQ=255)
                                                     gTCGCGACACCCCTCCGCGCTGTGCTGAATCCTCGCTGAAAAAATTTGAAGATAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTTAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTTAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGaat                                                                                                                                                                                           <  2:208729/251‑2 (MQ=255)
                                                           acacCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGAcc                                                                                                                                                                                     <  2:122191/251‑1 (MQ=255)
                                                               cccgccgCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTa                                                                                                                                                                                 <  1:414288/250‑1 (MQ=255)
                                                                    cgcgCAGTGCTAAATACTCGCTGAAACACTTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGACGCGTAGCGTACATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAACACGACGGGTTCGTGAAAGCGACCGTTAttcgg                                                                                                                                                                            >  2:382749/1‑247 (MQ=255)
                                                                          gTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGaa                                                                                                                                                                      <  1:392443/251‑1 (MQ=255)
                                                                                      cGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGcc                                                                                                                                                              >  1:310455/1‑247 (MQ=255)
                                                                                                 tttgaagagAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATgg                                                                                                                                                >  1:116837/4‑250 (MQ=255)
                                                                                                   tgaagagAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGaa                                                                                                                                              <  2:124029/249‑1 (MQ=255)
                                                                                                    gaagagAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGTCGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGa                                                                                                                                               >  1:434434/1‑248 (MQ=255)
                                                                                                     aagagAAATGTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGAc                                                                                                                                           <  2:123910/251‑1 (MQ=255)
                                                                                                        agaaaTTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGc                                                                                                                                        <  2:349473/249‑1 (MQ=255)
                                                                                                           aaTTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGAt                                                                                                                                     >  1:7983/1‑251 (MQ=255)
                                                                                                              ttGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGc                                                                                                                                   <  2:118697/250‑1 (MQ=255)
                                                                                                                gctgcgcAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGcgct                                                                                                                                >  1:86738/4‑248 (MQ=255)
                                                                                                                     gcagggtaaaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGAttt                                                                                                                           <  1:431139/244‑1 (MQ=255)
                                                                                                                        gggtaaaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTc                                                                                                                         >  1:353062/5‑250 (MQ=255)
                                                                                                                             aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAAc                                                                                                                   >  1:104220/1‑251 (MQ=255)
                                                                                                                                gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTc                                                                                                                <  2:104220/251‑1 (MQ=255)
                                                                                                                                gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTc                                                                                                                >  1:301032/1‑251 (MQ=255)
                                                                                                                                 ccGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTcc                                                                                                               >  2:10273/1‑251 (MQ=255)
                                                                                                                                       tgcttCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCt                                                                                                          <  1:182789/246‑1 (MQ=255)
                                                                                                                                           tCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGaaaacgacgagttc                                                                                                                                                                                                 <  2:301454/159‑11 (MQ=255)
                                                                                                                                           tCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGaaaacgacgagttc                                                                                                                                                                                                 >  1:301454/1‑149 (MQ=255)
                                                                                                                                             catcatTCGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCAGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGCCGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACcgcg                                                                                                   >  2:373996/1‑251 (MQ=255)
                                                                                                                                                catTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCGGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGGCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACcgcgcg                                                                                                 <  2:373766/250‑1 (MQ=255)
                                                                                                                                                           actacGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACt                                                                                     <  1:268832/251‑1 (MQ=255)
                                                                                                                                                            ctacGCGATGGACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTc                                                                                    >  2:335382/1‑251 (MQ=255)
                                                                                                                                                                  gATGGACCGCGATAACCGTTGGGATCGCTTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGcat                                                                              <  2:336798/251‑3 (MQ=255)
                                                                                                                                                                     ggACCGCGATAACCGTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATTGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCg                                                                           <  1:58041/251‑1 (MQ=255)
                                                                                                                                                                                   gTTGGGATCGCGTTGAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTAGTCGTGCAGAAGGTCAGCCAGATGCGGCAAGGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAGCGCAGATTTc                                                             <  1:69855/251‑1 (MQ=255)
                                                                                                                                                                                           cgcgTTGAAAAAGCCTATGCCCTGCTGACGCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCt                                                      <  1:375823/250‑1 (MQ=255)
                                                                                                                                                                                                 gAAAAAGCCTATGACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCTTcgcgcg                                               <  1:59556/251‑1 (MQ=255)
                                                                                                                                                                                                    aaaGCATATGACCCGCTGACCCTGGCGCAGGGTGTGTTCCAGGCCGATACTGCCATTCCTAGCTTACAGGCAGCTTATGCCCACGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCTTCCCGCGTa                                             >  2:205816/1‑250 (MQ=255)
                                                                                                                                                                                                             gACCTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCTTCGCGCGTAAGAAAGTGGt                                   >  2:367932/1‑251 (MQ=255)
                                                                                                                                                                                                               ccTGCTGACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCTTCGCGCGTAAGAAAGTGGTTa                                 <  1:402336/251‑1 (MQ=255)
                                                                                                                                                                                                                  gctgACTCTGGCGCAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCTTCGCGCGTAAGAAAGTGGTTAACg                              >  1:56504/1‑251 (MQ=255)
                                                                                                                                                                                                                              gcAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCTTCGCGCGTAAGAAAGTGGtt                                  >  1:44168/1‑235 (MQ=255)
                                                                                                                                                                                                                              gcAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACGCGTGCATTCGTAAACGCAGATTTCGACGGCTGCGCGCGTAAGAAAGTGGtt                                  <  2:44168/235‑1 (MQ=255)
                                                                                                                                                                                                                               cAGGGTGAGTTCCAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGAt                 >  2:240519/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggtgAGTTCCAGGCCGATACGGCCGTTGCTGGTGTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCt              <  2:67226/248‑1 (MQ=255)
                                                                                                                                                                                                                                          ccAGGCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAAt      <  1:130203/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gCCGATACTGCCGTTGCTGGTTTGCAGGCAGCTTATGCTCGCGACGAAAACGACGAGTTCGTGAAAGCGACCGTTATTCGTGCAGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCGCTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGCGAAATCACTCGTGCATTCGTAAACGCAGATTTCGACGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATAcgc  >  1:345565/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
ACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGC  >  NC_000913/3785667‑3786155

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