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→CAAscpCpropionyl‑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

TGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCTCGGCCGAT  >  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)
                                                                                                ggcgCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTTCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGTGCAATATCaataat                                                                                                                                                  <  1:74751/251‑1 (MQ=255)
                                                                                                             gtCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGgcgc                                                                                                                                     >  1:231688/1‑251 (MQ=255)
                                                                                                                        aTGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGa                                                                                                                           >  1:394097/1‑250 (MQ=255)
                                                                                                                           cAATGGAACGCGTCGGTACCCGCTATGTGCCAATCGCTCAGAACAAGCTTGACGCAGTCGTGGAAACAAACATGCCCGAAGCCGGTAATATGCTGGATAAGCAACATCCCATGTGCCAGCAACTTGCCGATAACATGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGCGGCGTGGGCAATATACAAAATGCGGTGATGGAGCGACTGGGGGAAAAc                                                                                                                       >  1:395406/1‑251 (MQ=255)
                                                                                                                                gATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGt                                                                                                                  <  1:243296/251‑1 (MQ=255)
                                                                                                                                        cGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATtcctc                                                                                                          >  1:299118/1‑251 (MQ=255)
                                                                                                                                         ggTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATtcctcc                                                                                                         <  1:212440/251‑1 (MQ=255)
                                                                                                                                             cccGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGttt                                                                                                     <  1:79305/251‑1 (MQ=255)
                                                                                                                                                gCTATGTGCAAATTGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTatg                                                                                                  <  1:93327/251‑1 (MQ=255)
                                                                                                                                                          aaTCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCgg                                                                                        <  1:163384/251‑1 (MQ=255)
                                                                                                                                                              gATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACGTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATGGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGGGGGCAATATCAATAATGCGGTAAGGGCTCGTGTTGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGt                                                                                    <  1:195603/251‑1 (MQ=255)
                                                                                                                                                               aTCCGAACAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGGAAAACCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTg                                                                                   >  1:173583/1‑251 (MQ=255)
                                                                                                                                                               aTCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTg                                                                                   <  1:257107/251‑1 (MQ=255)
                                                                                                                                                                       aaGATTGTCGCCGTCAGGGAAACCAACGTGCCCGACGCCGGTAGTATGCTGGATAAGCAAAATCCCATGTGCCAGGAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGAATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGa                                                                           <  2:4955/251‑1 (MQ=255)
                                                                                                                                                                                   gtcgtGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGTTGGTGCa                                                               <  1:171479/251‑1 (MQ=255)
                                                                                                                                                                                        ggAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTa                                                           <  1:113186/250‑1 (MQ=255)
                                                                                                                                                                                         gAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACt                                                         <  2:74927/251‑1 (MQ=255)
                                                                                                                                                                                          aaaCCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTg                                                        <  1:329772/251‑1 (MQ=255)
                                                                                                                                                                                                      ccGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCaa                                             >  1:136586/1‑250 (MQ=255)
                                                                                                                                                                                                        gACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAAt                                          <  1:435706/251‑1 (MQ=255)
                                                                                                                                                                                                                  atatGCTGGATAAGCAAAATCCCATGTGCCCACAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGAAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGGGCATTTATTGGAAACCGGCAAAATCAGCGGGGcc                                >  1:43661/1‑251 (MQ=255)
                                                                                                                                                                                                                  atatGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGCCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTTTGGGGGAAAACCCGGTACTTCCTCCGTTTATGATGGATTCGGAAGTGCTAAAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGcc                                >  1:286952/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGCATTTCTGCCGCTGCAAAGTGGCGTGGGCAACATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCCTTTACTGGAAACCGGCAAAATCAGCGGGGCCAgc                             >  2:80619/1‑251 (MQ=255)
                                                                                                                                                                                                                      gCTGGATAAGCAAAATCCCATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAAGTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAgc                             <  2:136586/250‑1 (MQ=255)
                                                                                                                                                                                                                        tggctAAGCAAAATCACATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCCATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTCTGATGTATTCGGCAGTGCTAAAGGCATCGGGTGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGcc                          >  2:449528/5‑251 (MQ=255)
                                                                                                                                                                                                                                      cccATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAAt            <  1:430547/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ccATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGAGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTTACAAtc           >  2:316918/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cATGTGCCAGCAAATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAAtct          >  1:355693/1‑251 (MQ=255)
                                                                                                                                                                                                                                               cagcaaATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCTCGGCCg    >  1:211783/1‑250 (MQ=255)
                                                                                                                                                                                                                                                agcaaATTGCCGATCACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCTCGGCCGAt  >  1:144008/1‑251 (MQ=255)
                                                                                                                                                                                                                                                agcaaATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCTCGGCCGAt  >  1:95418/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
TGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCTCGGCCGAT  >  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: