Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,065,449 G→A L216L (TTG→TTA 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,4490GA93.0% 158.9 / ‑5.3 58L216L (TTG→TTAscpCpropionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/27);  minor base C (2/2);  total (29/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.36e-01

CGGATGGTCGAGTCTGGTTAACCAGCGGGATCGGTAATGCGCCGACCTGGCTGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGG  >  NC_000913/3065202‑3065690
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   
cGGATGGTCGAGTCTGGTTAACCAGCGGGATCGGTAATGCGCCGACCTGGCTGCTGCGGTCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGAGCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTAccc                                                                                                                                                                                                                                                  <  1:38878/251‑1 (MQ=255)
cGGAGGGTCGAGTCTGGTTAACCAGCGGGATCGGTAATGCGCCGACCTGGCTGCTGCGGGCGAAGAAAGGGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTcccc                                                                                                                                                                                                                                                  <  2:341083/251‑1 (MQ=255)
     ggTCGAGTCGGTTTAACCAGCGGGATCGGTAATGCGCCGACCTGGCTGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATGGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGGCGCCGTCGTGGAAACCAACTTCCCCGACg                                                                                                                                                                                                                                              <  2:106169/250‑1 (MQ=255)
      gtctggTCTGGTTAACCAGTGGGATCGTGTATGCTGCGTCCGGGCTGCGGCGGGGTAAGAAAGTGATCAGTGAACGCAAGCACTAGCACGATCCGTGCGTGGCAGAACTGGCGGATAGTTTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAAGGGTTCGGGTCGGTACCCGCTATGTGCAGATCGATCCGAAAAAGATTGTCGCCGCGGTGGAAACCAACTTACCCGACGcc                                                                                                                                                                                                                                            <  2:319797/246‑1 (MQ=255)
        cGAGTCTGGTTAACCAGCGGGATCGGTAATGCGCCGACCTGGCTGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCgg                                                                                                                                                                                                                                          >  1:320684/1‑251 (MQ=255)
          aGTCTGGTTAACCAGCTGGATCGGTAATGCGCCGACCTGGCTGCTGCGGGCGAAGAAAGTGATCATAGAACTCAATCACTATCACGTTCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGTCGCCGGTa                                                                                                                                                                                                                                        <  2:242537/251‑1 (MQ=255)
              tGGTTAACCAGCGGGATCGGTAATGCGCCGACCTGGCTGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGCAAAAGATTGTCGCCGTCGTGGGAACCAACTTCCCCGACGCCGGtattat                                                                                                                                                                                                                                    >  1:47388/1‑251 (MQ=255)
                   aaCCAGCGGGATCGGTAATGCGCCGACCTGGCTGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTgg                                                                                                                                                                                                                               >  1:148263/1‑251 (MQ=255)
                        gCGGGATCGGTAATGCGCTGACCTGGCTGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAg                                                                                                                                                                                                                          <  1:365916/251‑1 (MQ=255)
                                    aaTGCGCCGACCTGACTGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATGGAGTGTCGCTCTTCCCTGCAATGGATCGCGGCGGTACCCGCACTGTGCAAATCGGTCCGGAAACAGATGGCGCCGTCGAGGAAACCACATTACCCGACGCCGGGCAAATGCTGGATacgaaacaagcca                                                                                                                                                                                                                >  2:313520/1‑238 (MQ=255)
                                       gcgcCGACCTGGCTGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCAtgtg                                                                                                                                                                                                            >  1:61035/1‑251 (MQ=255)
                                        cgcCGACCTGGCTGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGc                                                                                                                                                                                                           <  1:128438/251‑1 (MQ=255)
                                                     gctgcGGGCGAAGAAAGTGATCATTGAACTCAATCTCTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGa                                                                                                                                                                                            <  2:291574/251‑1 (MQ=255)
                                                                agaaAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCTCGCCCCCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAACTCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACACCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTc                                                                                                                                                                                  <  2:187485/251‑1 (MQ=255)
                                                                     gagATCATTGAACTCAATCACTATCTCGATCCGCGCGGTGCAGAACTGGCGGTTATTGTGATGCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGAGCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGAGAACGTGGTCACGtt                                                                                                                                                                            <  2:104026/249‑1 (MQ=255)
                                                                     gTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGtt                                                                                                                                                                            <  1:379275/251‑1 (MQ=255)
                                                                              gAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGc                                                                                                                                                                    <  1:212372/250‑1 (MQ=255)
                                                                                     aTCACTATCACGATCCGCGCGGAGCAGAACTGGCGGATATGGTTATGCCTGGCTCGCCACCTCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGCTCGCGTCGGTACCCGCAATGTGCAAATCGATCCGAAAAAGATTGTCGCCGGCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGAGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATg                                                                                                                                                            <  2:189258/251‑1 (MQ=255)
                                                                                      tCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGCTTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATgg                                                                                                                                                           >  1:400717/1‑251 (MQ=255)
                                                                                        aCTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGgcg                                                                                                                                                         <  1:288202/251‑1 (MQ=255)
                                                                                          tATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTAACCGCTATGGGCAAATCGATCCGCAAAAGATTGTAGCCGTCGCGGAAACAAAATTCCCCGACGACGGGAATATGCTGGTTAAGCAACATCCCATTTTCCAGCCGTTTGCCGCTAACGTGGCACCGTTTTTATTGCAGGAAATTgcga                                                                                                                                                        >  2:250001/1‑249 (MQ=255)
                                                                                          tATCACGATCAGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCCCCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGGTCGCGTCGGTACCCGCTATGGGCAACTCGATCCGAAAAAGGTTGTCGCCGGCGTGGAACACAACTTTCCCGACGCCGGAAATAAGACGGATAAGTAAAATACCCTGTGCCAGCAGGTTGCCGATAACGTGGTCACGTTCTTATTGCAGGGAATGgcgc                                                                                                                                                        >  2:515266/1‑250 (MQ=255)
                                                                                             cACGATCCGCGCGTGGCAGAACTGGCGGATATTGTGATTCCTGGAGCGCCACCGCGGCTCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTCCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATg                                                                                                                                                     <  2:148263/250‑1 (MQ=255)
                                                                                               cGATCCGCGCGTTGCAGAACTGGCGGATATTGTGTTTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCCATGGATCGCGTCGGTACCCCCTATGTGCAAATCGATCCAAAAAAGATTTGCGCCGGCGTGGAAACCAACTTACCCGACGCCGGAAAAATGATGGATAAACAAAATCACATGTGCCAACAGAATGCCGATAACGGAGTCACGTTCTTAtaacagggaattgccgcatggg                                                                                                                                                  >  2:328485/1‑230 (MQ=255)
                                                                                                   ccGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGACGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTAt                                                                                                                                              >  1:524277/1‑251 (MQ=255)
                                                                                                    cgcgcgTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGGACCCGCTATGGGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACAAAATTACCCGACGCCGGAAATATGCTGGCTAAGCAAAATCCCATTTGCCAACATATGTGCCGTAACGTGGGCACGTTATTATTGCTGGAAATGGCGCGTGGGcagatt                                                                                                                                             >  2:141872/1‑246 (MQ=255)
                                                                                                          ttGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCCATAGCGGGTCGATCTTCCATGCAATGGATCGAGTCGGTACCCGCCATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTTGAAACCAACTTACCCGACCCCGGTAATATGCTGGATAAGCAAAATCCCAAGTGCCAGCAGGTTGCCGATAACGTGGTCACGTTCTTATTGCCGGAATAGGCGCATGGGCGTATTccgccg                                                                                                                                       >  2:78314/1‑251 (MQ=255)
                                                                                                                    ggCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATTGCGCATGGGCGTATTCCGCCGGAATTTCt                                                                                                                               <  1:116345/249‑1 (MQ=255)
                                                                                                                    ggCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCt                                                                                                                               <  1:479027/249‑1 (MQ=255)
                                                                                                                           atTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGc                                                                                                                       >  1:24466/1‑250 (MQ=255)
                                                                                                                              gggagtCCTGGCGCGCCACCGCGGCCCAATAGAGTGTCGATCTTCCATGCAAGTGATCGCGTCGGTACCCGCTATGTGAAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGCGTATTCCGCCGGAATGTCTGCCGCTGCAAAGt                                                                                                                  <  2:203327/246‑1 (MQ=255)
                                                                                                                                  gtcgtggggcgcaccctcggcgCAATAGCGGGGCGAGCTTCCATGCAATGGATCGCGTCGGTACCCGAGATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATAGGCTGGATAAGCAAAATTCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCTCATGGGCGTATTCCGCCGGAATTTCTGCCGCAGCAAAGTggt                                                                                                               <  2:59007/234‑2 (MQ=255)
                                                                                                                                      tGGCGCGCCACCGAGGCGCAATAGCGGGTCGATCTTCCATGCAAGGGAACGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTATCGCCGTCGAGGAAACCAAATTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGCTAACGTGGTAACGTTCTTATTGCAGCCAAAGTCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTAGCGTgg                                                                                                           >  2:367855/1‑251 (MQ=255)
                                                                                                                                                    ggcgCAATAGCGTGTCGATCTTCCATGCACTGGATCGCGTCGGTACCCGCTATGTGCAAATCGCTCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCACATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGACGGAATTACTGCCGCTGCAAAGTGGCGTGGGCAATATCaataat                                                                                             >  2:256124/1‑251 (MQ=255)
                                                                                                                                                       gcAATCGCGTGTCGATCTTCCATGCAATGGATCTCGTCGGTACCCGCTATGTGCAAATCGATCCCAAAAAGAGTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATACCATTTGCCAGCATATTGCCGATAACGTGGTCATGTTCTTATTGCAGGCAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCATCTGCAAAGTGGCGTGTGCAATTTTAATAATGCg                                                                                          <  2:276047/251‑1 (MQ=255)
                                                                                                                                                              cGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTAACCGACGCCGGTAATATGCTGGATAAGCAAAATACCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATgg                                                                                   >  1:352645/1‑251 (MQ=255)
                                                                                                                                                                         tCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCt                                                                            <  1:526040/247‑1 (MQ=255)
                                                                                                                                                                          ccATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTAGTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTggggg                                                                       <  2:130717/251‑1 (MQ=255)
                                                                                                                                                                            aTGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGa                                                                      >  1:98540/1‑250 (MQ=255)
                                                                                                                                                                                 aTGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAAccc                                                                >  1:219065/1‑251 (MQ=255)
                                                                                                                                                                                    gATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCggg                                                             >  1:197061/1‑250 (MQ=255)
                                                                                                                                                                                       cgcgTCGGTACCCGCTATGTGCAACTCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAAt                                                          >  1:461894/1‑251 (MQ=255)
                                                                                                                                                                                          gTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCAGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGCATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATtcc                                                       >  1:25444/1‑251 (MQ=255)
                                                                                                                                                                                              gTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCg                                                   >  1:102785/1‑251 (MQ=255)
                                                                                                                                                                                                   cgcgatgggCAAATCGATACGAAAAAGATTGTCGACGGCGTGGAAACCAACTTACACGACGCCGGTAATATGCTGTATAAGCAAAATCCCAGGTGCCAGCAGAGTGCCGATAACGTGGTCACGTTCTTATTGCATGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCGGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGTTAATGGCGCGTCTGGGGGAAAACCCGGGAATTCCTCCGTTTat                                              <  2:25444/243‑1 (MQ=255)
                                                                                                                                                                                                       aTGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTatgatg                                          >  1:28884/1‑251 (MQ=255)
                                                                                                                                                                                                             aaaTCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCg                                    >  1:504514/1‑251 (MQ=255)
                                                                                                                                                                                                                   aTCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTg                              <  1:529227/251‑1 (MQ=255)
                                                                                                                                                                                                                     ccGAAAAAGATTGTCGCCGTCGTGGAATACAACTTACCCGACGCCGGGAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCt                            <  2:12482/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gCCGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTTTGATGTATTCGGAAGTGCTACAGGAATCggtggg             >  2:73987/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ccGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTgg                                                                                                           <  1:367805/156‑1 (MQ=255)
                                                                                                                                                                                                                                      ccGTCGTGGAAACCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGACGGAATTTCTGCCGCTGCAAAGTGGCGTgg                                                                                                           >  2:367805/1‑156 (MQ=255)
                                                                                                                                                                                                                                             ggAAACCAACTTACCCGACGCCTGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTa      <  2:504514/250‑1 (MQ=255)
                                                                                                                                                                                                                                               aaaCCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAg                                <  1:263949/222‑1 (MQ=255)
                                                                                                                                                                                                                                               aaaCCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAAACCGGTAATTCCTCCGTTTATGATGTATTCGGAAg                                >  2:263949/1‑222 (MQ=255)
                                                                                                                                                                                                                                               aaaCCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAAACCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTg   <  2:28884/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aaCCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCtt                                                                                                                                                                         <  1:227965/84‑1 (MQ=255)
                                                                                                                                                                                                                                                aaCCAACTTACCAGGCGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCtt                                                                                                                                                                         >  2:227965/1‑84 (MQ=255)
                                                                                                                                                                                                                                                 aCCAACTTACCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTgg  >  1:75309/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   
CGGATGGTCGAGTCTGGTTAACCAGCGGGATCGGTAATGCGCCGACCTGGCTGCTGCGGGCGAAGAAAGTGATCATTGAACTCAATCACTATCACGATCCGCGCGTTGCAGAACTGGCGGATATTGTGATTCCTGGCGCGCCACCGCGGCGCAATAGCGTGTCGATCTTCCATGCAATGGATCGCGTCGGTACCCGCTATGTGCAAATCGATCCGAAAAAGATTGTCGCCGTCGTGGAAACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGG  >  NC_000913/3065202‑3065690

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