Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,225,271 G→A V482I (GTA→ATA)  metH → cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,225,2710GA98.5% 237.7 / NA 67V482I (GTA→ATA) metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/37);  minor base C (1/0);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCT  >  NC_000913/4225036‑4225516
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
tGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACg                                                                                                                                                                                                                                        >  2:92181/1‑251 (MQ=255)
tGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACg                                                                                                                                                                                                                                        >  1:376156/1‑251 (MQ=255)
            cGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACGGGGACa                                                                                                                                                                                                                              <  1:132858/249‑1 (MQ=255)
            cGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACa                                                                                                                                                                                                                              <  2:389972/249‑1 (MQ=255)
            cGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAg                                                                                                                                                                                                                             <  1:164521/250‑1 (MQ=255)
                  cGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACTGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATa                                                                                                                                                                                                                      <  1:247898/251‑1 (MQ=255)
                     aGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTc                                                                                                                                                                                                                   <  1:201111/251‑1 (MQ=255)
                     aGCGGCGATGGAGCGTTTTCTCAAACTGATTGGCGGTAAACTGGATAACGCACGAGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCAGTGAAAAAGGTCTGAAGTGTATCCCGGGCAAAGGCATTGTTAACTCTCTCTCGAGGAAAGACGGCGTCGATTCCTTTATCCATCACGCGAAATTGGTGCGTAGCTACGGTGCGGCAGTGGTGATACTGGCCCTTGACGAAACGGGAGAGGACGATACTc                                                                                                                                                                                                                   >  2:20181/1‑251 (MQ=255)
                      gcggcgATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCAGTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTcg                                                                                                                                                                                                                  <  2:533665/251‑1 (MQ=255)
                       cggcgATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTcgc                                                                                                                                                                                                                 >  1:364352/1‑251 (MQ=255)
                        ggcgATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGACCAGGGACAGGCCGATACTcgcg                                                                                                                                                                                                                <  2:352614/251‑1 (MQ=255)
                               gTGCGTTTTCTCATGTTGATTGCCGGTGAACCGGATATCGCGCGCGTGCCGATTATGATCGACTCCTAAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATTACTCTAAATGGTGTCGTCGCTACTGTGCGGCAGTGGTTACAATGGCCTTTGACGATCAGGGACAGGCTGATACTCTCGCCCGGaa                                                                                                                                                                                                         <  2:266437/251‑1 (MQ=255)
                                 gCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGaaaa                                                                                                                                                                                                       <  2:201882/251‑1 (MQ=255)
                                          cAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGGCGCTACGGTGCGGCAGTGGTGATAATGGCCTGTGACGAACAGGGACAGGCCGATACGCTCGCACGGAAAATCGAGAtt                                                                                                                                                                                               <  2:76272/250‑1 (MQ=255)
                                          cAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGGATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCTTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAtt                                                                                                                                                                                               <  1:282771/250‑1 (MQ=255)
                                            aTCTGATTACCGGTGAACCGAATATCGCCCGAGTGCCGATTATGATCGACTCATCAAAATGGGACGTCATTGAACCCGGTCTGACGTGTATCCAGGGCAACGGCATTGTTACCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGCTGATAATGGCCCTTGACGAACAGGGACAGCCCGATACTCGCGCACGGAACATCGAGATTTGc                                                                                                                                                                                            >  2:201855/1‑251 (MQ=255)
                                              cTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAttt                                                                                                                                                                                              >  1:166458/1‑247 (MQ=255)
                                              cTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAttt                                                                                                                                                                                              <  2:166458/247‑1 (MQ=255)
                                                 aTTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCg                                                                                                                                                                                       <  1:315429/251‑1 (MQ=255)
                                                  ttGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCgg                                                                                                                                                                                      >  1:106064/1‑251 (MQ=255)
                                                    gCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGc                                                                                                                                                                                    <  2:97108/251‑1 (MQ=255)
                                                           aaCCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCTGATACTCGCGCACTGAAAATCGATATTTGCCGTCGGGCGTACaaa                                                                                                                                                                             <  1:406181/251‑1 (MQ=255)
                                                              cGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATc                                                                                                                                                                          >  2:411696/1‑251 (MQ=255)
                                                                     cgctcgcGAGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTTGTGATAATGGCCTTTGACGAACAGGGACAGGCCGAGACTCGCGCACGGAAAATCGAGATGTGCCGGCGGGCGTACAAAATCCTCAcc                                                                                                                                                                    <  2:230873/250‑1 (MQ=255)
                                                                           cgTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAgg                                                                                                                                                             >  1:496845/1‑251 (MQ=255)
                                                                            gTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGt                                                                                                                                                            >  2:105493/1‑251 (MQ=255)
                                                                              gCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTg                                                                                                                                                          <  2:350853/251‑1 (MQ=255)
                                                                                cGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGc                                                                                                                                                        >  1:37508/1‑251 (MQ=255)
                                                                                cGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGc                                                                                                                                                        >  1:429245/1‑251 (MQ=255)
                                                                                    tATGATCGACTCCTCAACATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTcc                                                                                                                                                    >  1:348764/1‑251 (MQ=255)
                                                                                     aTGATCGACTCCTCAAAATGGGACTTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTccc                                                                                                                                                   <  1:144336/251‑1 (MQ=255)
                                                                                          cGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGcc                                                                                                                                                >  2:211705/1‑249 (MQ=255)
                                                                                               cctcAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCaga                                                                                                                                             >  1:104270/1‑247 (MQ=255)
                                                                                                    aaaTGGGACGGCATTGAAAAAGTTCTGAAGTGTAACAAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGCCGATGCCTTTAGCCATCACGCGAAATTGTTGCGTCGCTACGGGGCGGAAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatcat                                                                                                                                    <  2:13359/251‑1 (MQ=255)
                                                                                                           aCGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGAc                                                                                                                             <  1:92181/251‑1 (MQ=255)
                                                                                                            cGTCATTGCAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGACATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCCAGGAAAATCGAGAttt                                                                                                                                                                                              >  2:462902/1‑185 (MQ=255)
                                                                                                            cGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAttt                                                                                                                                                                                              <  1:462902/185‑1 (MQ=255)
                                                                                                                   gAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACAt                                                                                                                     >  1:478106/1‑251 (MQ=255)
                                                                                                                   gAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACAt                                                                                                                     >  1:415644/1‑251 (MQ=255)
                                                                                                                            cTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGt                                                                                                            >  2:200173/1‑251 (MQ=255)
                                                                                                                             tGAAGTGCATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCgggc                                                                                                           >  2:494891/1‑249 (MQ=255)
                                                                                                                                aGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGGGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCa                                                                                                        <  1:484958/251‑1 (MQ=255)
                                                                                                                                    tATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCtt                                                                                                                                                                                                                                             >  1:85373/1‑114 (MQ=255)
                                                                                                                                    tATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCtt                                                                                                                                                                                                                                             <  2:85373/114‑1 (MQ=255)
                                                                                                                                       ccAGGGCCAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGc                                                                                                  >  2:155535/1‑250 (MQ=255)
                                                                                                                                       ccAGGGCACAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGc                                                                                                  >  2:483348/1‑250 (MQ=255)
                                                                                                                                        cAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCAt                                                                                                <  2:15606/251‑1 (MQ=255)
                                                                                                                                          gggCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTg                                                                                              >  2:446704/1‑251 (MQ=255)
                                                                                                                                              aaaGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGa                                                                                             >  2:136191/1‑248 (MQ=255)
                                                                                                                                               aaGGCATTGTTATCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGTCTTTGACGAACATGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAagag                                                                                         <  2:335619/251‑1 (MQ=255)
                                                                                                                                               aaGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAagag                                                                                         <  2:377290/251‑1 (MQ=255)
                                                                                                                                                                aTCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTAcgcg                                                                          <  1:174185/249‑1 (MQ=255)
                                                                                                                                                                           agagGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCCTTGAAGAGCACAACAACTACGCGCAGGACTTTATCg                                                             >  2:504500/1‑251 (MQ=255)
                                                                                                                                                                              gggCGTGGATGCATTTATCCTTCACGCGACATTGTTGCGTCGCTACGGTGCGGCAGTGGTGCTAATGGCCTTTGCCGACAAGGGACAGGCCGATACTCGCACACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGCAGCGGTTGGCTTCCCGCCAGAAGATCTCATCTTCGCCCCAAACATCTTCGCGGTCGCACCGGGCATTGTAGAGACCAACAACTACACGAACGACTTTATCGgcg                                                          >  2:259585/1‑251 (MQ=255)
                                                                                                                                                                              gggCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCGCCGAAGAGGTTGGCTTCCCGCCAGAAGATAGCATCTTCGACCCAAACAGCGTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACGTTATCGgcg                                                          <  2:281931/251‑1 (MQ=255)
                                                                                                                                                                              gggCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGgcg                                                          >  2:433716/1‑251 (MQ=255)
                                                                                                                                                                                         ccTTTATCCATCACGCGACATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGCAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACCACAACTACGCGCAGGCCTTTATCGGCGCGTGTGAAGa                                                >  2:336125/1‑250 (MQ=255)
                                                                                                                                                                                            ttATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTGCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCTGCGCGTGTGAAGACa                                              <  1:416581/249‑1 (MQ=255)
                                                                                                                                                                                              aTCCATCACGCGAAATTGTTGCGTCGCGACGGTGCGGCAGGGGTGATAATGGCCTTTGACGAACAGGGAAAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCGTCGACCCAAACATCTTCGCGGTCGCAACTGGCAGTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCaa                                          <  2:37508/251‑1 (MQ=255)
                                                                                                                                                                                                ccATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCTCCGAAGAGGTTGGCTTCCCGCCAGATGATATCATCTTCGTCCCAAACATCTGGGCGGTCGCAACTGGCATTGATGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAAc                                        <  1:529458/251‑1 (MQ=255)
                                                                                                                                                                                                ccATCACGCGAAATTGTAGCGTCGCTTCGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAAc                                        <  2:147155/251‑1 (MQ=255)
                                                                                                                                                                                                    cacgaataaTTGTGGCGTCTCTAAGGTGCTGCAGTGGTGATAATTGCCTTTCAGGAACAGGTACAGGCAGATACTCGCGCACGGAAAATAGAGATTTGCCGTCGGGAGTACAAAATCCTCACCGAAGAGGTTGGCTTCCTGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACGACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGa                                    <  2:301886/244‑1 (MQ=255)
                                                                                                                                                                                                    cACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAAcgcg                                     >  1:416491/1‑250 (MQ=255)
                                                                                                                                                                                                          aaaTTGTTGCGTCGCTACGGTGCGGCAGTGGTGAGAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCTCAACAACTACGCGCATGACTGTATCGGCGCTTGTGAAGACGTCAAACGCGAACTGcc                              <  2:524477/251‑1 (MQ=255)
                                                                                                                                                                                                                 tgCGTCGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCAcgc                        <  1:411696/250‑1 (MQ=255)
                                                                                                                                                                                                                      cGCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGa                   <  2:283829/250‑1 (MQ=255)
                                                                                                                                                                                                                       gCTACGGTGCGGCAGTGGTGATAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGAtt                 <  1:105493/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tggtgaTAATGGCCTTTGACGAACCGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCGTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACAGCCGCACGCGCTGATTTCCGgcggcgacact  >  2:289697/1‑246 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
TGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTCAAAATGGGACGTCATTGAAAAAGGTCTGAAGTGTATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACATCAAACGCGAACTGCCGCACGCGCTGATTTCCGGCGGCGTATCT  >  NC_000913/4225036‑4225516

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: