Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 421,256 A→G T91A (ACC→GCC)  proY → putative transporter ProY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913421,2560AG97.1% 236.1 / NA 68T91A (ACC→GCC) proYputative transporter ProY
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (38/28);  minor base T (1/1);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCAT  >  NC_000913/421011‑421494
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
gtgggtaaaGTCCCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACAGTGCCGGc                                                                                                                                                                                                                                            <  1:448733/244‑1 (MQ=255)
gggggtaaagcccccGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGc                                                                                                                                                                                                                                            <  1:215397/239‑1 (MQ=255)
gTGGGCTAAGTCCCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGc                                                                                                                                                                                                                                            <  1:25669/251‑1 (MQ=255)
gTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGc                                                                                                                                                                                                                                            >  1:133481/1‑251 (MQ=255)
gTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGc                                                                                                                                                                                                                                            >  1:542451/1‑251 (MQ=255)
       aaGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCt                                                                                                                                                                                                                                     >  1:234794/1‑251 (MQ=255)
          cctccGCCACATACGCTTTATGGCACTGGGTTCAGCAATGGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACt                                                                                                                                                                                                                                  <  2:601855/248‑1 (MQ=255)
           aCCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACTGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACtt                                                                                                                                                                                                                                 >  1:555853/1‑250 (MQ=255)
            cccGCCACATACGCGTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGc                                                                                                                                                                                                                                <  2:234794/251‑1 (MQ=255)
             ccGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACt                                                                                                                                                                                                                                  >  1:601855/1‑248 (MQ=255)
                    aTACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGACAGACGAAATCAAAATGGCCGGTCAGAGCGTGGTGTTGGCCTATATTCTCGGAGGTCTCGCGGCGTATATCATTATGCGTGCGCTGGCGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCTGGTTAACTTGCCGGCTGGACCTCCTGCTTTGTAAt                                                                                                                                                                                                                       >  2:529339/1‑251 (MQ=255)
                               ggCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCGGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTc                                                                                                                                                                                                             <  1:96491/250‑1 (MQ=255)
                               ggCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCg                                                                                                                                                                                                            >  1:20022/1‑251 (MQ=255)
                               ggCACTGGGTTCAGCAATTGGCACAGGGCTGTTTTAAGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATCTTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGGCCTACTGCTTTGAAATCCTTATTGTCg                                                                                                                                                                                                            >  2:60996/1‑251 (MQ=255)
                                       gTTCAGCAATTGGCAACGGGCTGTTTTACGGTTCGGAAGACGCCATCAAAATGGCCGGTCTGAGCGTGTTGATGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGACGCCAGCTCTTTCTCGCGTTATGATCAGGAAAACCTCGGCCCGATGGCAGGTTACATTGACGGCTGGACCTACGGCTTTGAAATCCTTATTGTCGCCATCGcc                                                                                                                                                                                                    >  2:47714/1‑251 (MQ=255)
                                          ccgCAATTGTCCCCGGGCTGTTTTCCGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGTGTGTTTTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCTCTGGAAAGCCTCGGCCCGCGGGCAGGTGACATTTCCGCCGGGACCTACTGCTTTGAAATCCTTTTTGTCACTAGCGCCGAt                                                                                                                                                                                                 <  2:350693/249‑1 (MQ=255)
                                            gCAATTGGCACAGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGAtgt                                                                                                                                                                                               >  1:464430/1‑251 (MQ=255)
                                                    cACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGACGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGAAAAGATCGGTCCATAACCGGGCCGCCAGCTTTTTCTCGCGTTATGCCCAGGAAAACTTCGGCCCGCGGCCAGGTAACTTTGCCGGCGGGACCTCTGCCTTGGAAATCCTTTTTTTCGCCATCGCCGAGTTGACCGCtt                                                                                                                                                                                       >  1:363134/1‑251 (MQ=255)
                                                                   cGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGCAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATgg                                                                                                                                                                        >  2:170331/1‑251 (MQ=255)
                                                                       tCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGAGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGgtgt                                                                                                                                                                    >  1:545091/1‑251 (MQ=255)
                                                                        cGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATAACATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGATTTGAAATCCTTTTTGTCGCCATAGCCGTTGTGACCGCTTTTGGTTTCTATTTggggttt                                                                                                                                                                   >  2:365832/1‑247 (MQ=255)
                                                                         ggCAGACGACATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTCCTGCTTTGACATCCTTATTGTCGCCATCGCCGCTGTGACAGCTTTTGGTATCTCTATGGgtgttt                                                                                                                                                                  >  2:508654/1‑249 (MQ=255)
                                                                          gCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGGCATACTGCTTTGAAATCCTTATTGTCGCCATCGCCGGTGTGACCGCTTTTGGTATCTATATGGGTGTCTg                                                                                                                                                                 >  2:602894/1‑251 (MQ=255)
                                                                              aCGCCATCAAAATGGCCTGTCCGCGCGTCTTGTTGGCCTATATTATCGGAGATATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTTCCGGCTGGACCTACTGCTTTGAACTCCTTATTGTCGCCATCGCCGATGTGACAGCTTTTGGTATCTATATGGGTGTCTGGTTc                                                                                                                                                             >  2:56788/1‑251 (MQ=255)
                                                                              aCGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGTACCTACTGCTTTGAAATCCTTATTGTCTCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTc                                                                                                                                                             <  1:365446/251‑1 (MQ=255)
                                                                                       aaaTGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTg                                                                                                                                                    <  1:145112/251‑1 (MQ=255)
                                                                                        aaTGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGTGTGTCTTGTTCCCGACGGTGc                                                                                                                                                   <  1:481517/251‑1 (MQ=255)
                                                                                               ggTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTg                                                                                                                                            <  1:444998/251‑1 (MQ=255)
                                                                                                    gAGCGTGTTGTTGGCATATATTATCGGTGGTATCGCGGCGTCTCTATTTATGCGTGCGCTGGGGGCAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGCCCTCATGCTTTGAAATCCTTATTGTCGCCATCGCCGCTGTGCCCGATTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGACTCACTGGAttt                                                                                                                                       >  2:597771/1‑251 (MQ=255)
                                                                                                            tgttgGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACt                                                                                                                                <  1:137247/250‑1 (MQ=255)
                                                                                                                     atatTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCgtggtg                                                                                                                      >  1:179453/1‑251 (MQ=255)
                                                                                                                          atCGGTGGTTTCGCGGCTTATATCATTAGGCTTGCGCTGGGGGAAATGGCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCATGAAATCCGCGGCCCTCTGGCAGGTTACATTGCCGTCTGGGCCGACTGCTTTGAAATCCTGGTTGTCGCCATCGCCGGTGTTACCGCTTTTGGGATCTATATGGTTGACTGGGTCCCCACTGTGCCGCTCTGGATTTGGGTACTGAGCGTGGTGCTGa                                                                                                                  <  1:290343/250‑1 (MQ=255)
                                                                                                                              gtggtATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCAtt                                                                                                             <  2:332348/251‑1 (MQ=255)
                                                                                                                              gtggtATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGGTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGTTTTGGGGACTGAGCGTGGTGCTGATCAtt                                                                                                             >  1:495435/1‑251 (MQ=255)
                                                                                                                                ggtATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTg                                                                                                           >  2:139270/1‑251 (MQ=255)
                                                                                                                                       cggcgTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTg                                                                                                           >  1:231706/1‑244 (MQ=255)
                                                                                                                                       cggcgTATATCATTATGCGTGCGCGGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTg                                                                                                           <  2:231706/244‑1 (MQ=255)
                                                                                                                                             atatCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGAGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCtg                                                                                              >  1:248188/1‑251 (MQ=255)
                                                                                                                                                    tATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCCCTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCg                                                                                       >  2:325879/1‑251 (MQ=255)
                                                                                                                                                            cgcTGGGGGAAATGTCGGTACATATCCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAGCCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTa                                                                               <  2:84558/251‑1 (MQ=255)
                                                                                                                                                                gggggAAATGTCGGTACATAAACCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCg                                                                           >  2:365081/1‑251 (MQ=255)
                                                                                                                                                                  gggAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTACCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGt                                                                         >  1:476334/1‑251 (MQ=255)
                                                                                                                                                                   ggAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTAt                                                                              >  2:200573/1‑245 (MQ=255)
                                                                                                                                                                   ggAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTAt                                                                              <  1:200573/245‑1 (MQ=255)
                                                                                                                                                                   ggAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCGGATGAGCGTGAAGGTATTCGGTg                                                                        <  2:39798/251‑1 (MQ=255)
                                                                                                                                                                     aaaTGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTgg                                                                                                                                     >  2:421351/1‑188 (MQ=255)
                                                                                                                                                                     aaaTGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTgg                                                                                                                                     <  1:421351/188‑1 (MQ=255)
                                                                                                                                                                      aaTGCCGGTACATACCCCGGCCGACATCTCTTTCTCGCGTTATGCCCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGTCTGGACCTACTGCTTTGAAATCCTTATTGTCTCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGCTCCCGACGGTGCCGCACTGGATTTGGGTACGGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGTTGAGc                                                                     <  2:41528/251‑1 (MQ=255)
                                                                                                                                                                       aTGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCt                                                                    <  1:523052/251‑1 (MQ=255)
                                                                                                                                                                                 aTAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTgg                                                          >  1:499574/1‑251 (MQ=255)
                                                                                                                                                                                   aaCCCGGCCGCCCGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCCGGACCTACTGCTTTGAAATCCTTCTTGACGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGTGCGTGAAGGTATTCGGTGAGCTGGACTTCTGGtt                                                        >  1:387858/1‑251 (MQ=255)
                                                                                                                                                                                    aCCCGGCCGCCCGCTCTTTATCGCGGTATGCGAAGGAAAACCTCGGCACGCCGGCAGGTTACATTGCCGGCTGGAGCTCCTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGCCGGGGCGGCCCTGGATTTGGGTACTGAGGGTGGTGGTGATCATTTGGGCCGTAAACCTGATGAGAGTGAAGGGATTCGGTGAGCTGGAATTCTGGTtc                                                       >  2:391925/1‑251 (MQ=255)
                                                                                                                                                                                        ggCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGt                                                   <  1:30011/251‑1 (MQ=255)
                                                                                                                                                                                         gccgccAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGtt                                                  <  1:508654/251‑1 (MQ=255)
                                                                                                                                                                                          ccgccAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTTCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGTTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTc                                                 <  1:396887/251‑1 (MQ=255)
                                                                                                                                                                                                 ctctTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTaa                                            >  1:391875/1‑249 (MQ=255)
                                                                                                                                                                                                   ctTTCTCGCGTTATGCGCAGGACAACCTCGGCCCGCTGGCAGGTTCCATTGCCGGCTGGACCTACTGCTTTGAACTCCTTATTGTCGCCATCGCCGATGTGACAGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAAACTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTtct                                                      >  2:244474/1‑237 (MQ=255)
                                                                                                                                                                                                   ctTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGTTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTtct                                                      <  1:244474/237‑1 (MQ=255)
                                                                                                                                                                                                          gcgTTATGCGCAGGCAAACCTCGGCCCGCTGGCAGGTTAATGTGCCGGCTGGACCTGCTGCTTTGAACTCGTGATTGTCGCCATAGCCGATGAGACCGCTTTTGGTCTATATGTGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGCGCGTGGTGCTGATAATTTGCGCCGTAAACCTGATGAGCGAGAAGGTATTCGGTGATGTGGATGTCTGGTTCTCGTTGTTTACAGTCGccaccg                                 >  2:375228/1‑250 (MQ=255)
                                                                                                                                                                                                          gcgTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACTGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGTTGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTAAAGTCGccacc                                  >  2:413699/1‑250 (MQ=255)
                                                                                                                                                                                                                   gcAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGACGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTAAAGTCGCCACcatcatcatc                         >  2:429874/1‑250 (MQ=255)
                                                                                                                                                                                                                   gcAGGAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTAAAGTCGCCACcatcatcatc                         >  2:140753/1‑250 (MQ=255)
                                                                                                                                                                                                                      ggAAAACCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCCCTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTAAAGTCGCCACCATCATCATCATg                      >  2:65456/1‑250 (MQ=255)
                                                                                                                                                                                                                         aaaCCTCGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATGTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCGCGTTCTTTAAAGTCGCCACCATCATCATCATGATTg                  <  1:353186/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGGCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGcc             <  1:597771/250‑1 (MQ=255)
                                                                                                                                                                                                                                ggCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGt                                                                                <  2:170046/182‑1 (MQ=255)
                                                                                                                                                                                                                                ggCCCGCTGGCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGt                                                                                >  1:170046/1‑182 (MQ=255)
                                                                                                                                                                                                                                         gCAGGTTACATTGCCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAAGTCTGGTTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGGTTCGGcat  <  1:424581/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
GTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCAT  >  NC_000913/421011‑421494

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