Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,064,898 G→C 12.8% S33S (TCG→TCC scpC → propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,064,8980GC12.8% 96.6 / 2.5 55S33S (TCG→TCCscpCpropionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (28/16);  new base C (3/5);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

TATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGGCGCGGTGAGTTTCGTTGACCTGCATTTGAG  >  NC_000913/3064677‑3065122
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 
tATTCAGGGGATGCGCCGCGCGGCGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGGTCGCCGaa                                                                                                                                                                                                                              >  2:2826/1‑227 (MQ=255)
  gcaaGGAGAGGCGCCGCGCGGTGTAAGACAGCAAAAATTACGAGGAAGGGATGAACGTTTCCCTCGAAAAACGAAAACCTAATTGCGTTGGTCATTAATCCCGGCGAACGAAGGAGTAAAAATGAAAACTCAGTGAAAAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGc                                                                                                                                                                                                                           <  2:173537/225‑1 (MQ=255)
          aTGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCCAAAGCCTTAccc                                                                                                                                                                                                                    >  2:611983/1‑227 (MQ=255)
          aTGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGAAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGAGGCATTTAGCGGCTTTACCCCGGCGGGGTCGCCGAAACCCCTAccc                                                                                                                                                                                                                    >  2:418051/1‑227 (MQ=255)
           tGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGCAATTATCCACCATAACGACATGGTGGCATTTAGCGGCTTTCCCCCGGCGGGTTACCCGAAACCCCTACCCAc                                                                                                                                                                                                                  >  1:552774/1‑228 (MQ=255)
           tGCGCCGCGCGGTGGATGACAGCGAAGCTTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGACATTAATCCCTGCGAACGAAGGAGTCAAAATGGTAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCACAACGACATGGTGGCATTTAGCGGCTTTACCCCGTCAGGTTCTACCAAACCCCTCCCCa                                                                                                                                                                                                                   >  2:645252/1‑227 (MQ=255)
                                 cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGACCGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGGTGGCCGAAAGACCCACCCCACCCCATGTCCCCCAGAGCCAAc                                                                                                                                                                                            >  1:482675/1‑228 (MQ=255)
                                 cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGACCGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGCAAGACATTATCCAGCAAAATGCCTGGTGGCAATTTACCGGGCTTACCCCGCCGGGGTCTCCCAAAACCCTAACCCCcccgctagtcccccggccaaca                                                                                                                                                                                            >  1:21270/1‑210 (MQ=255)
                                 cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGCGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGAATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCAACCGCAATTAACCGCAGGGCTAAc                                                                                                                                                                                            >  2:438734/1‑228 (MQ=255)
                                 cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGGTCGCCGCAAGCCCTACCCACCGCGATTGCCCGCAGGGCTaag                                                                                                                                                                                            >  2:350706/1‑227 (MQ=255)
                                 cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCACGGCGGGTTCTCCGAAAACCCTACCCACCGCGATTGCCCGCAAGGCTAAc                                                                                                                                                                                            >  2:336144/1‑228 (MQ=255)
                                 cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGAATTTAGCGGCTTTACCCCGGCGGGTTTGCCGAAAGCCCCACCCCCCGCGACTGCCCGCAGAGCTaa                                                                                                                                                                                             >  2:406500/1‑227 (MQ=255)
                                 cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCAAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTAACCACCCACATTGCCCGccagactaac                                                                                                                                                                                            >  1:181601/1‑219 (MQ=255)
                                 cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATACAGCATAACGACATGGGGGCATTTAGCGGCTTTACCCCAGCGGGTTCCCCGAAAGCCCTCCCCACCGCGATTGCCCGccacgcaaac                                                                                                                                                                                            >  2:321690/1‑219 (MQ=255)
                                   aaGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAAccg                                                                                                                                                                                          >  2:226181/1‑226 (MQ=255)
                                   aaGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTAGCCGCAAGCCCTACCCACCGCGGTTGCCCGCAGAGCTCACg                                                                                                                                                                                           >  2:574500/1‑227 (MQ=255)
                                                aGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGCATGCCCGCCGAGCTAACGAACAGCATGAGGc                                                                                                                                                                              >  2:455797/1‑227 (MQ=255)
                                                 ggatgtaaccttttCATCAAAAAACGAAACCATATTTCGTTGGGTCAGTACACCCTGAAAACAAGGAAGTAAAAATGGAACCACAGGGAAAAAGGAGGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCa                                                                                                                                                                            <  1:604777/217‑1 (MQ=255)
                                                                                ccatagtgttgttgcatatAACCCCCGCGACAGAAGGAAAAAAAAGAGAACACCAGTGAAAAGGGTTGGCCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACgg                                                                                                                                              <  1:611983/210‑1 (MQ=255)
                                                                                 tAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACggg                                                                                                                                             >  2:571549/1‑227 (MQ=255)
                                                                                         ttGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGAAGGGTGCggaaa                                                                                                                                     >  2:410479/1‑223 (MQ=255)
                                                                                          tGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGTCGGGTTCGCCGAAAGCCCTACCCACCGCGATTTCACGCAGGGCTCACGCACAGCATGAGGCCAAAAAGCCGTATCAAATTGCCCTTCTGACGGGTGAGTCAATc                                                                                                                                   >  2:169643/1‑228 (MQ=255)
                                                                                          tGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGCCGGGGGCGTCAATc                                                                                                                                   >  2:371059/1‑228 (MQ=255)
                                                                                                  tAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCCATTGCCCGCAGGGCTACAGAAAAGCATGAGGCCCAAAAACCGGATCAAAATCGCCATCTGACggggggggcaaaaagccccg                                                                                                                            >  1:497935/1‑211 (MQ=255)
                                                                                                       ccTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTTCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTTACGGGTGCGTCAATCAGCCCCGCCGCt                                                                                                                       >  2:499532/1‑228 (MQ=255)
                                                                                                          gCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGGATCAAATTCGCCTTCTGACGGGTGCGCCACTCAGCTCCGCCCCTg                                                                                                                      >  2:327335/1‑226 (MQ=255)
                                                                                                                aaagaaaaaaaaaGGGAAACCCGGGGGACAGGGAGTGCCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTa                                                                                                              <  2:482675/220‑1 (MQ=255)
                                                                                                                                            aagaggtacccgaaaagaaaGCGGCAGAAATTATCCAGCCGATAGACATGGGGGCATTTGGGCGCTTTGCCCGGGCGGGTGCGCCGAAATCCATGCCCACCGGGGTTGCCCGAAGAGCTGACGAACAGCATGAGGCAAAAAAGCCGTTTAAAATTCGCCTTCTGACGGGTGCGGCAATCAGCGCCGACGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGc                                                                                  <  1:438734/210‑1 (MQ=255)
                                                                                                                                               ggtgccccccctaTGAAGCGGCAGAAATTATCCAGCATCACGACCGGGTGTCATTTAGCGTCGTTTCCCGGGCGGGTTACCCGAAAGCCCTCCCCGCCCGTATTCCCCGCAGAGCTAACGAACAGCAGGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTAACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTg                                                                               <  1:2826/216‑1 (MQ=255)
                                                                                                                                                   cccccccAATGAGGCGGCAAAAATTATCCAGCATACCGACATGGTGTCATTTAGCGGCTTTACCCCGGGGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGcc                                                                           <  1:571549/223‑1 (MQ=255)
                                                                                                                                                   cacccgaaaTGAAGGGGAAAAAATTTCCCATCATAGCGACATGGGGGCATTTGGCGGCCTGACCCCGGCGGTTTACCAGAAAGCCCCACCCCCCGCAATTGCCCCCCGAGTTAACGAACAGCATGGGGCCAAAACGCCGTACCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGcc                                                                           <  1:227137/221‑1 (MQ=255)
                                                                                                                                                    cccctcAATGAGGGGGGAGAAATGTCCCCGCATAAAGACAGGGTGCCATTTAGCGGCTTTCCCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGGGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGcc                                                                           <  1:321690/222‑1 (MQ=255)
                                                                                                                                                    cccGCCAATGAACGGGCAGAAATTCTCCGGCAGAACCACATGGTGGCCTTTGGCGGTTTTACCCCGGGCGGTTCGCCGAAAGCCCCACCCACCGCGATTGCCCGCAGCGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGcc                                                                           <  1:350706/226‑1 (MQ=255)
                                                                                                                                                    cccGCCAAAGACGGGGCAAAAATGACCCAGCAGAACGAAAGGGGGGAATGTAGCGGCTTTACCCAGGCGGGTGCGCAAAAAGCCCTCCCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAGGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCTCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGcc                                                                           <  1:523668/226‑1 (MQ=255)
                                                                                                                                                           aTGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGGTCGCCGAAAGCCCTACCCCCCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCCAACAGCCGTATCCAATTCGCCTTCTCACGGGAGCGGCAATCAGCGCCGCCGCTGCCGCTGGTCTTTCTATCCCCGCGGCTGTTTCCCGGCGTGCGCCCAATCaaa                                                                  >  1:279588/1‑228 (MQ=255)
                                                                                                                                                             taaggtggaaaaaTTAACCATGATAAGGGTATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCCCCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATACCAAAc                                                                 <  1:226181/217‑1 (MQ=255)
                                                                                                                                                                   gCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGAGGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTcc                                                          >  2:441637/1‑228 (MQ=255)
                                                                                                                                                                        aaTTATCCCGCCTAACGTATTGGTGGCTTTTAGAGCTTTTCCCCCGGAGGGTTCCCCGAAACCCCTACCCACCGCGATTGCCCGCAGATCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTCCGGttt                                                     <  1:264494/228‑1 (MQ=255)
                                                                                                                                                                               caagaaaaCGAGAGGGTGGCTTTAAGGGGCTTTCCCACGGGGGATGTCCCAAAAGCCCTCCCCACGGGGATTGCCCGCAGGGCTAACGAACAGCATGAGGCCAAAAAGCCTTATCAAATGCGCCTTCTGACGGGTGCGTCAATAAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTCCGGTTTACGTaa                                               <  1:657681/221‑1 (MQ=255)
                                                                                                                                                                                 gCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGAGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCCTATCAAACATCGTCCGGTTTACGTaaaaa                                            >  2:490868/1‑228 (MQ=255)
                                                                                                                                                                                  cATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTTCCCGCCGAGCTAACGCACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCCGCGCCCCCGCTGACGATGTACTTTCTGACCCCGATGCTGTTTCCTGGCGTGCGCCATACCCAACATCGCCCGGTTTCCCTaaaaaa                                           >  1:539676/1‑227 (MQ=255)
                                                                                                                                                                                  cATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGGGCGTCAATCAGCGCCGCCGCTGACCATGGACCTTCTTACGCCGGAGCCGTTTCCCGGCGGGCGCCAATTTCAAAAAACTCCGGGTTAAGaaaaaaa                                           >  1:304522/1‑227 (MQ=255)
                                                                                                                                                                                  cATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCCCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGAt                                                                                               >  2:181090/1‑177 (MQ=255)
                                                                                                                                                                                   aaaaCGACATGGTGGCTTTTAGCGGCTTTCCCCGGGCGGGTTCGCCGAAAGCCCTACCCACCGGGATTGCCCGAAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTCCGGTTTACGTaaaa                                             <  2:497935/223‑1 (MQ=255)
                                                                                                                                                                                   aTAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCCCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATg                                                                                              <  1:181090/177‑1 (MQ=255)
                                                                                                                                                                                            tggtggCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCCAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGGGCGCCATATCCAACATCGTACGGGTTACGGAAAAAGAAAAATCAgg                                 >  1:175667/1‑228 (MQ=255)
                                                                                                                                                                                            tggtggCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGc                                                                                                                                                        >  2:334713/1‑110 (MQ=255)
                                                                                                                                                                                            tggtggCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTCCGGTTTACGTAAAAAAATCAATCAgg                                 >  1:390265/1‑228 (MQ=255)
                                                                                                                                                                                             ggtggCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGcc                                                                                                                                                       <  1:334713/110‑1 (MQ=255)
                                                                                                                                                                                              gtggCATTTAGGGGCTTTCCCCCGGGGGGTTAGCCGAAAGCCCTACCCACCGCGATTGCCCGCGGGGCTAACGAACAGCAGGAGGCAAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGgc                               <  1:495463/228‑1 (MQ=255)
                                                                                                                                                                                              gtggCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAggg                                <  2:390265/227‑1 (MQ=255)
                                                                                                                                                                                                    ctttGGCGGCTTTCCCCCGGCGGGTTCGCAGAAAGCCCTCACCACCCGGATTGCCCGAAGGGCTAACGAACAGCAGGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGGCGCGGt                          <  2:246507/227‑1 (MQ=255)
                                                                                                                                                                                                          cGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGGGCGCCATCTCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGGCGCGGGGAGttt                    >  2:174983/1‑228 (MQ=255)
                                                                                                                                                                                                                 cccccGGGGGGTTCGCCGAAAGCCCCACCCCCCGCGTTTGCCCGCAGGGCTAACGAACAGCATGAGGCCAAAACGCCCTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGGCGCGGTGAGTTTCGTTGAc             <  1:174983/227‑1 (MQ=255)
                                                                                                                                                                                                                            ttCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCCTGAGGCCCAAACGCCGTATCCAATTCGCCTTCTGACGGGTGCGGCAATCAGCGCCGCCGCTTACGATGTACTTTCTGTCGCCGATTCCGTTTCCTGTGGGGCGCCATATCAAACATCGTCCGGCTTTAGTAACAAGCGCAATCCGGGAGGGGCGAGATTCGTGGACCCCCATTTGAg  >  1:102691/1‑228 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 
TATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTTGGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCCATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGGCGCGGTGAGTTTCGTTGACCTGCATTTGAG  >  NC_000913/3064677‑3065122

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