Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 147,213 2 bp→TT coding (246‑247/903 nt) rpnC → recombination‑promoting nuclease RpnC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913147,2130CT97.3% 108.2 / NA 38H82H (CAC→CATrpnCrecombination‑promoting nuclease RpnC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (13/23);  minor base G (0/1);  total (13/25)
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.94e-01
*NC_000913147,2140CT97.3% 112.2 / NA 38Q83* (CAA→TAA) rpnCrecombination‑promoting nuclease RpnC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (13/23);  minor base G (0/1);  total (13/25)
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.94e-01

GCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCAT  >  NC_000913/146974‑147446
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                        
gCACCTAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCACGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGcc                                                                                                                                                                                                                                <  2:425893/251‑1 (MQ=255)
          ccACCCCGCAAGACGCGGTATTTAAACAATGTTTAATGCATTCGGAGGCGGCTCTCCACTTTCTGAAGATACGTTTGACAGTTGTATTACGCGAACTTTGTGACCTCAACACGCTTTATTTAGAGTCGGGTAGTGTCATTGAACAGAGCCTGAAAGTACACATCTCCGATGTGCGCTATTCCGTTCAATTGCTGGGCATTCCCAGGTATCAGCATGGGGTGATTTAACATGAAAGCCATCCGGATAAg                                                                                                                                                                                                                         <  2:177211/248‑1 (MQ=255)
            acacCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGCACATTAAAGCAAGCCGGATACGaaaa                                                                                                                                                                                                                     >  2:453094/1‑250 (MQ=255)
                    tGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGCACTTTGTGACCTCAACACCCTTCATTTAGAGTCGGGGAGATTCATTGAAGAGAGCCTGAAAAGACCCAGCCCGGACGTGCTCTATTACGTGCCAATGCAGGGCCATACCGGTTATCTGCATGTTGTGATTGAACATTAAAGCTAGCCGGATAAGAAAATGGCCTTTc                                                                                                                                                                                                            >  1:305975/1‑251 (MQ=255)
                    tGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGGGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTCTCTGCATGTTGTGATTGAACATTAaagcaagc                                                                                                                                                                                                                                 >  2:361473/1‑230 (MQ=255)
                    tGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAaagcaagc                                                                                                                                                                                                                                 <  1:361473/230‑1 (MQ=255)
                     gACGCGGTATTTAAACAATTTTGAATGCATGCGGAGACGGCTCGCTACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACGCGCTTCATTTAGAGGCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCg                                                                                                                                                                                                           <  2:303697/251‑1 (MQ=255)
                      aCGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGc                                                                                                                                                                                                          <  1:126229/251‑1 (MQ=255)
                        gcgGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCat                                                                                                                                                                                                        <  1:28434/251‑1 (MQ=255)
                                tAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGtt                                                                                                                                                                                                >  1:327745/1‑251 (MQ=255)
                                 aaaCAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGtt                                                                                                                                                                                                <  2:162587/250‑1 (MQ=255)
                                           tAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAgc                                                                                                                                                                                      <  1:108106/250‑1 (MQ=255)
                                            aaTGCAGGCGGATACGGCTCGCATCTTTTTGGAGATATCTTCGCCAGTGGAAGTGCGCGATCTTTGTGACCTCAACACGCTTATATTAGAGTCGGGGAGTTTTATTGAAGATAGCCTGAAAGGACACAGCGCGGACGTGTTCTTTTCCGTGCAAATGCAGGGCAATCCCGGTTAACTGCATGTTGTGATTGGACATTAAAGCAAGCCTGACAAGAAAATGGCCGGTCGCATGATGTTTTCTTCTATAgcc                                                                                                                                                                                     <  2:309535/250‑1 (MQ=255)
                                             atgcatgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGCAGAGAGCATGAAAGGACACAGCAAGGAAGTGCTCTATTCCGTGCAAATGCAGGGCAAACCCGGCTATCTGCATGTTGTGATTGAACATTAAAGCAAACCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAgccgc                                                                                                                                                                                   >  1:432404/1‑251 (MQ=255)
                                               gcatgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCa                                                                                                                                                                                 <  2:456303/251‑1 (MQ=255)
                                                 atgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGACAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGGGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATg                                                                                                                                                                               >  1:280638/1‑251 (MQ=255)
                                                           ggCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATc                                                                                                                                                                     <  1:57542/251‑1 (MQ=255)
                                                               cgcgACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGAAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTg                                                                                                                                                                   <  1:306936/249‑1 (MQ=255)
                                                                                       ccAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGGTATCTGCATGTTGTGGTTGAACATTAAAGCAAGACGGATAAGAAAATGGCCTTTCGCATGATGAGTTCTTCTATAGCCGCCATGCAACGGCATCTGGAGGATGACCACGATAAGCTGCCGCt                                                                                                                                         >  1:157198/1‑251 (MQ=255)
                                                                                          gTGGAATTACGCGACCTTTGTGACCTCAACACGCTTCAGTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCTTGCAAATGCAGGTCAATCCCGGTTATCTGCATGTTGTGAGTGAACAGTAAATCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGAACCGGCATCTGGAGGCTGACCACGATAAGCTGCCTCtggt                                                                                                                                      <  2:355655/251‑1 (MQ=255)
                                                                                          gTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCTCtggt                                                                                                                                      <  1:310677/251‑1 (MQ=255)
                                                                                                                aCCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAg                                                                                                                <  1:397901/251‑1 (MQ=255)
                                                                                                                            ttCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCg                                                                                                                                                                          >  1:2389/1‑181 (MQ=255)
                                                                                                                              cATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACAcc                                                                                                  >  1:414207/1‑251 (MQ=255)
                                                                                                                                         cGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTa                                                                                        <  1:187287/251‑1 (MQ=255)
                                                                                                                                         cGGGGAGTGTCTTTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTa                                                                                        <  2:335739/251‑1 (MQ=255)
                                                                                                                                                           agagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGat                                                                      >  1:216627/1‑251 (MQ=255)
                                                                                                                                                                 tGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAACTGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGttt                                                                >  2:30051/1‑251 (MQ=255)
                                                                                                                                                                    aaGGACACAGCACGCACCTGCTCTATTCCGCGCCCCTGCTGGGAAATCCAGCTTACCTGCAAGTTGTGATTGAACATTACAGCAACCCGGATAACAAAATGCCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGCGGTGCCGATACTGTTTTATCAGGCCGAGGCCACACCTTCTCCGCTATCAATGTGCTGGTTTGATATGTTTTa                                                              >  2:55804/1‑250 (MQ=255)
                                                                                                                                                                               ccGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTTCTGGTTTGATATGTTTTACTCGCCGGAGCt                                                  <  1:402366/250‑1 (MQ=255)
                                                                                                                                                                                              tCCGTGCAAATGCAGGGCAATCCCGGTTAGCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCGACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTGTACTCGCCGGAGCTGGCGCGACGCGTCt                                    <  2:39142/250‑1 (MQ=255)
                                                                                                                                                                                               ccGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCta                                   <  1:89909/250‑1 (MQ=255)
                                                                                                                                                                                               ccGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCta                                   <  1:304074/250‑1 (MQ=255)
                                                                                                                                                                                                   gCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCGGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAAc                               <  1:453094/250‑1 (MQ=255)
                                                                                                                                                                                                                   cccGGTTATCTGCAGGTTGTGAGTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTGCGCATGATGCGTTAGTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGTTGTTTCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCGATGTGCTGGTTTGATATGTTTTACTCGCCGGGGCTGGCGGGATGGGTCTATAACCGTCCGTTCCCGCtggt              <  2:42535/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCtggtg             <  1:101305/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGAGGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCa      <  2:74712/248‑1 (MQ=255)
                                                                                                                                                                                                                               cATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCAt  >  1:75442/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                        
GCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCAT  >  NC_000913/146974‑147446

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