Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,065,497 G→A Q232Q (CAG→CAA)  scpC → propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,065,4970GA96.3% 182.7 / NA 54Q232Q (CAG→CAA) scpCpropionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/33);  minor base C (2/0);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCTCGGCCGAT  >  NC_000913/3065253‑3065743
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
ggctgcGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGcc                                                                                                                                                                                                                                                   <  1:37145/250‑1 (MQ=255)
 gctgcGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCGGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCg                                                                                                                                                                                                                                                  <  1:275249/251‑1 (MQ=255)
         cGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTAGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTg                                                                                                                                                                                                                                          >  1:77653/1‑251 (MQ=255)
            agaaAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTc                                                                                                                                                                                                                                       <  1:347047/251‑1 (MQ=255)
                   gATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTc                                                                                                                                                                                                                                 <  1:182325/250‑1 (MQ=255)
                       aTTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGTCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCAGGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAATCCAACTTGCCCGACGCCGGTAATATGCTGGAGAAGCAAAATCCCAGGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTAtt                                                                                                                                                                                                                            <  1:259600/251‑1 (MQ=255)
                        ttGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGCGATTCCTGGCGCGCCACCGCGGCGCAATAGCGAGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTATCGCCGTCGTGGAAACCAACTTACCCGACCCCGGTAATATGCTGGATAAGCAAAATCCCATGGGCCAGCACAGTGCCGaaaaagtgtgaacgtttttatct                                                                                                                                                                                                                           >  2:97216/1‑229 (MQ=255)
                           aaCTCAATCACGATCGCGATCCGCGCGTGGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCTACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCATATTCCCATGTGCCAGCAAATGGCCGATAACGTGGTCACGTTCTTATTGCAg                                                                                                                                                                                                                        <  2:350047/251‑1 (MQ=255)
                           aaCTCAATAACTATCACGAACAACGCGTTGCAGAACTGGCGGATATTGTGATTCCTGCCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAACGATTGGAGCCGTCGTGGAAACCAAATTGCCCGCCGCCGGTACTATGCTGGATAAGCAAAATCCCATGTGACAACACATTGCCGTTAACCTGTTCACGTTCTTATTGc                                                                                                                                                                                                                          >  2:433735/1‑249 (MQ=255)
                              ccAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGaa                                                                                                                                                                                                                     <  2:324187/250‑1 (MQ=255)
                                        ccACGGTCCGCGCGTGGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCGTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAAGGGCGCATg                                                                                                                                                                                                           <  2:208356/250‑1 (MQ=255)
                                          aCGATCCGCGCGTTGCAGAACTGGCGGATATTGTGCTTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGCTCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGGAACCCACTTTCCCGCCGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAACTGCCCATAACGTGGTCCCGTTCTTATTGCAGGAAATGGCGCAgggg                                                                                                                                                                                                         >  2:384261/1‑251 (MQ=255)
                                           cGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTTCCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGc                                                                                                                                                                                                        <  1:32122/251‑1 (MQ=255)
                                              tCCGCGCGTTGCAGAACTGGCGGATATTGTGATTACTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCACTGGATCGCGTCGGAACCCGCTATGGGCAAAACGATCCGAAAAAGATTGTCGCCGTCGGGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGAAAAACAAAATCCCCTGTACCAGCAAATTGCCGCTAACGGGGTACCGTTCTTATTGCAGGAAATGGCGCCTGGGCGt                                                                                                                                                                                                      >  2:306901/1‑250 (MQ=255)
                                               ccGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTAt                                                                                                                                                                                                    <  1:350949/251‑1 (MQ=255)
                                                  cgcgTTGCAGAACGGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTGGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCTATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTcc                                                                                                                                                                                                 <  1:382975/251‑1 (MQ=255)
                                                        gCAGAACTGGCGGATATGGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTAGTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGa                                                                                                                                                                                           <  2:61787/251‑1 (MQ=255)
                                                         taaaaCTGGCGGATATTGTGATTACCGGCGCGCCACCGCGGCGCAATAGCGTGGCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCCAATCGATCCGCAAAAGCTTGTCGCCGACGGGGCAACCAACATGCCCGACGCCGGAAACATGCTGGATAAGCAAAATCCCATGTGCCATCAAATTGCCGATAAAGTGGTCACGCTCTCAATGCAGGAAATGGCGCATGGACGGATTccgcctgtc                                                                                                                                                                                          >  2:147812/4‑247 (MQ=255)
                                                            aaCTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAAttt                                                                                                                                                                                       <  1:407377/251‑1 (MQ=255)
                                                             aCTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGAAATTTc                                                                                                                                                                                      <  1:126580/251‑1 (MQ=255)
                                                                     atatTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTg                                                                                                                                                                              <  1:179682/251‑1 (MQ=255)
                                                                             aTTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTgg                                                                                                                                                                      >  1:74927/1‑251 (MQ=255)
                                                                                   ggCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTgg                                                                                                                                                                 <  1:56953/250‑1 (MQ=255)
                                                                                   ggCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGGTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAAGTTCTGCCGCTGCAAAGTGGCGTg                                                                                                                                                                  <  1:384261/249‑1 (MQ=255)
                                                                                     cgcgcCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCAGAAGCCGGTAATATGCTGGATACGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGAAATTTCTGCCGCTGCAAAGTGGCGTGGGc                                                                                                                                                               >  1:143589/1‑250 (MQ=255)
                                                                                      gcgcCACCGAGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCa                                                                                                                                                              >  2:430547/1‑250 (MQ=255)
                                                                                                ggcgCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTTCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGTGCAATAT‑Caataat                                                                                                                                                  <  1:74751/251‑1 (MQ=255)
                                                                                                             gtCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGgcgc                                                                                                                                     >  1:231688/1‑251 (MQ=255)
                                                                                                                        aTGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGa                                                                                                                           >  1:394097/1‑250 (MQ=255)
                                                                                                                           cAATGGAACGCGTCGGTACCCGCTATGTGCCAATCGCTCAGAACAAGCTTGACGCAGTCGTGGAAACAAACATGCCCGAAGCCGGTAATATGCTGGATAAGCAACATCCCATGTGCCAGCAACTTGCCGATAACATGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGCGGCGTGGGCAATATACAA‑AATGCGGTGATGGAGCGACTGGGGGAAAAc                                                                                                                       >  1:395406/1‑251 (MQ=255)
                                                                                                                                gATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGt                                                                                                                  <  1:243296/251‑1 (MQ=255)
                                                                                                                                        cGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATtcctc                                                                                                          >  1:299118/1‑251 (MQ=255)
                                                                                                                                         ggTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATtcctcc                                                                                                         <  1:212440/251‑1 (MQ=255)
                                                                                                                                             cccGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGttt                                                                                                     <  1:79305/251‑1 (MQ=255)
                                                                                                                                                gCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTatg                                                                                                  <  1:93327/251‑1 (MQ=255)
                                                                                                                                                          aaTCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCgg                                                                                        <  1:163384/251‑1 (MQ=255)
                                                                                                                                                              gATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACGTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATGGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGGGGGCAATAT‑CAATAATGCGGTAAGGGCTCGTGTTGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGt                                                                                    <  1:195603/251‑1 (MQ=255)
                                                                                                                                                               aTCCGAACAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGGAAAACCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTg                                                                                   >  1:173583/1‑251 (MQ=255)
                                                                                                                                                               aTCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTg                                                                                   <  1:257107/251‑1 (MQ=255)
                                                                                                                                                                       aaGATTGTCGCCGTCAGGGAAACCAACGTGCCCGACGCCGGTAGTATGCTGGATAAGCAAAATCCCATGTGCCAGGAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGAATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGa                                                                           <  2:4955/251‑1 (MQ=255)
                                                                                                                                                                                   gtcgtGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGTTGGTGCa                                                               <  1:171479/251‑1 (MQ=255)
                                                                                                                                                                                        ggAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTa                                                           <  1:113186/250‑1 (MQ=255)
                                                                                                                                                                                         gAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACt                                                         <  2:74927/251‑1 (MQ=255)
                                                                                                                                                                                          aaaCCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTg                                                        <  1:329772/251‑1 (MQ=255)
                                                                                                                                                                                                      ccGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCaa                                             >  1:136586/1‑250 (MQ=255)
                                                                                                                                                                                                        gACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAAt                                          <  1:435706/251‑1 (MQ=255)
                                                                                                                                                                                                                  atatGCTGGATAAGCAAAATCCCATGTGCCCACAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGAAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGGGCATTTATTGGAAACCGGCAAAATCAGCGGGGcc                                >  1:43661/1‑251 (MQ=255)
                                                                                                                                                                                                                  atatGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGCCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTTTGGGGGAAAACCCGGTACTTCCTCCGTTTATGATGGATTCGGAAGTGCTAAAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGcc                                >  1:286952/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGCATTTCTGCCGCTGCAAAGTGGCGTGGGCAACAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCCTTTACTGGAAACCGGCAAAATCAGCGGGGCCAgc                             >  2:80619/1‑251 (MQ=255)
                                                                                                                                                                                                                      gCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAAGTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAgc                             <  2:136586/250‑1 (MQ=255)
                                                                                                                                                                                                                        tggctAAGCAAAATCACATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CCATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTCTGATGTATTCGGCAGTGCTAAAGGCATCGGGTGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGcc                          >  2:449528/5‑251 (MQ=255)
                                                                                                                                                                                                                                      cccATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAAt            <  1:430547/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ccATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGAGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTTACAAtc           >  2:316918/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAAtct          >  1:355693/1‑251 (MQ=255)
                                                                                                                                                                                                                                               cagcaaATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCTCGGCCg    >  1:211783/1‑250 (MQ=255)
                                                                                                                                                                                                                                                agcaaATTGCCGATCACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCTCGGCCGAt  >  1:144008/1‑251 (MQ=255)
                                                                                                                                                                                                                                                agcaaATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCTCGGCCGAt  >  1:95418/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
TGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATAT‑CAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCTCGGCCGAT  >  NC_000913/3065253‑3065743

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: ‑