Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,783,749 A→G 7.1% H240R (CAC→CGC)  fadK → short chain acyl‑CoA synthetase, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,783,7490AG 7.1% 158.5 / 2.7 56H240R (CAC→CGC) fadKshort chain acyl‑CoA synthetase, anaerobic
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (2/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

TAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTTTGCGGCGGAACCACAATCCCCAAAAAAGTGGCGCGTGAATGCCAGCAGCGCGGCATTAAATTATTAAGTGTTTATGGTTCCACAGAAAG  >  NC_000913/1783449‑1784046
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         
tAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCa                                                                                                                                                                                                                                                                                                           <  2:168779/301‑1 (MQ=255)
   tCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGc                                                                                                                                                                                                                                                                                                        >  2:215902/1‑301 (MQ=255)
    ccTGCCGCTGCAAAATCAGCTTCCACAGCTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCa                                                                                                                                                                                                                                                                                                       <  1:215659/301‑1 (MQ=255)
       gtcGCTGCAAAATCAGCTTCCCCAACTACTGCAAATTGGCGGCGTGGACAAACTGGCTCCCGCCACCTCGTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACg                                                                                                                                                                                                                                                                                                    <  1:522377/299‑1 (MQ=255)
             gCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCttt                                                                                                                                                                                                                                                                                              >  2:410836/1‑301 (MQ=255)
             gCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGAAGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCttt                                                                                                                                                                                                                                                                                              >  2:419339/1‑301 (MQ=255)
                 aaTCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGc                                                                                                                                                                                                                                                                                          <  2:184281/301‑1 (MQ=255)
                  aTCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGTTGCTGTTTACCTCCGGACCCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGTTGTCTTTATGCTTCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCa                                                                                                                                                                                                                                                                                         >  2:17690/1‑301 (MQ=255)
                      gCTTCCACAACTACAACAAATGTTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGAAAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGc                                                                                                                                                                                                                                                                                     <  2:508176/301‑1 (MQ=255)
                           cacaACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAAc                                                                                                                                                                                                                                                                                >  2:453722/1‑301 (MQ=255)
                                  ccccttcacgacgcttttccgatctcTGTCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCGCGCAACGGGCTTTCTGCATGGCGTAACGGCACCa                                                                                                                                                                                                                                                                         <  1:341624/276‑1 (MQ=255)
                                   caacaaATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCaa                                                                                                                                                                                                                                                                        >  2:203910/1‑300 (MQ=255)
                                    aacaaATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCAtt                                                                                                                                                                                                                                                                       >  2:527831/1‑301 (MQ=255)
                                            gTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTg                                                                                                                                                                                                                                                               >  2:70416/1‑301 (MQ=255)
                                             ccGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTgg                                                                                                                                                                                                                                                              <  1:453158/300‑1 (MQ=255)
                                                gCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGCATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCGCGCAACGGGCTTTCTGCATGGCGTAACGCCAACATTCTTAATTGccgc                                                                                                                                                                                                                                                           >  1:25975/1‑301 (MQ=255)
                                                    ggACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGACCCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGCCTGAATCTGACCTGGCAGGATGCCTTTATGCTGCCTGCCCCACTTGGTCACGCAACGGGCTTTCTGCATGCCGTAACGGAACCATTCTTAATTGccgcccgc                                                                                                                                                                                                                                                       >  1:50019/1‑301 (MQ=255)
                                                           cTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCGCGCAACGGGCTTTCTGCATGGCGTAACGGCAACaagatcggaagagcacacgtctgaaa                                                                                                                                                                                                                                                >  2:342087/1‑276 (MQ=255)
                                                            tGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGtgtt                                                                                                                                                                                                                                               <  2:564050/301‑1 (MQ=255)
                                                               cTCCTGTCACCTTTTTCTTCTCATTAAGGCAGCTTTCCGCCGACAATACCTCACTGACCACGGCGATAACGACCAACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTgttgtt                                                                                                                                                                                                                                            <  1:615313/301‑1 (MQ=255)
                                                                      cACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATAtt                                                                                                                                                                                                                                     <  1:94379/301‑1 (MQ=255)
                                                                       aCCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATAttt                                                                                                                                                                                                                                    >  2:60473/1‑301 (MQ=255)
                                                                              cccTCTCATTAAGTCGGATTATCGCCGACAATACCTCACTGACCACGGCGATACCGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGGCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTcc                                                                                                                                                                                                                             <  1:410273/301‑1 (MQ=255)
                                                                                                     gCCGACAATGCCTCGCTGACCACGGCGATAACTACCCACGGCGATTAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTg                                                                                                                                                                                                      <  1:60405/301‑1 (MQ=255)
                                                                                                                cTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGt                                                                                                                                                                                           >  2:232179/1‑301 (MQ=255)
                                                                                                                  cacgggtcACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTg                                                                                                                                                                                         <  1:70326/294‑1 (MQ=255)
                                                                                                                                  aaCGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCCCGgc                                                                                                                                                                         >  1:191272/1‑301 (MQ=255)
                                                                                                                                         cACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTATCGCCAGTGAGCGGGCTTATTGCGAGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGCCGCTCGCAGCGCGTTGTTAGACATTTTCAATCCAGATCCGTTTTTCCCGCAGCTTGAGCAACAACGCTGCACCCGTAAGcacagaccaaacc                                                                                                                                                                  >  1:450941/1‑289 (MQ=255)
                                                                                                                                            ggtgATGAATTGGCTGCGGTGTTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCGccgagcatacaggtgcaacgctgc                                                                                                                                                               <  1:231936/298‑23 (MQ=255)
                                                                                                                                                  gAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCt                                                                                                                                                         <  2:390566/301‑1 (MQ=255)
                                                                                                                                                      tAGTTGCTGTGCTGTTTACCTCCTGAACCGAGGGTCTGTCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGa                                                                                                                                                     <  2:108162/301‑1 (MQ=255)
                                                                                                                                                                      tACCTTCGGAACCGAGGGTCTGCCAAGGGGCGTGATGCTAACGCATAATAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGTGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTa                                                                                                                                     <  1:573121/301‑1 (MQ=255)
                                                                                                                                                                      tACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGAAACGCCGTTTGTCTATGATCTTTTGAATGTACTa                                                                                                                                     >  2:405640/1‑301 (MQ=255)
                                                                                                                                                                        ccTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTaga                                                                                                                                   >  2:547531/1‑301 (MQ=255)
                                                                                                                                                                        ccTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCGATGATCTTTTGAATGTACTaga                                                                                                                                   >  2:450323/1‑301 (MQ=255)
                                                                                                                                                                         ttCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTagag                                                                                                                                  <  2:191515/300‑1 (MQ=255)
                                                                                                                                                                             ggAACCGGGGGTCTGCCAAAGGGCGTGATGCTAACGTATAACTATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAaac                                                                                                                              <  1:203667/301‑1 (MQ=255)
                                                                                                                                                                                   gAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCGCGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCcg                                                                                                                        <  2:26026/301‑1 (MQ=255)
                                                                                                                                                                                      ggTCTGCCAAATGTCGTGATGCTAACGCATAACAATTTTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCgg                                                                                                                     <  2:451505/301‑1 (MQ=255)
                                                                                                                                                                                             cAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTc                                                                                                              >  2:121866/1‑301 (MQ=255)
                                                                                                                                                                                                gctcttccgatcTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAgc                                                                                                           <  2:493982/290‑1 (MQ=255)
                                                                                                                                                                                                   gggtggtgcttatgctgtaaaaTGTTCTCGCTAGGTGGCGGGTTTATGGCGCGTGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCGGCGCAACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCt                                                                                                        <  1:418776/281‑1 (MQ=255)
                                                                                                                                                                                                    ggTGATGCTTACGCATAACAATTTTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTGATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTg                                                                                                       <  1:546858/300‑1 (MQ=255)
                                                                                                                                                                                                          gCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGtttt                                                                                                 >  1:253728/1‑300 (MQ=255)
                                                                                                                                                                                                             aaCGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGCTCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCaacagatcggacacgc                                                                                              >  1:493336/1‑286 (MQ=255)
                                                                                                                                                                                                                   tAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCAGACCTTTCAGCGCTGCGTTTCTTTCTTTgt                                                                                        >  1:581279/1‑300 (MQ=255)
                                                                                                                                                                                                                     aCAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTTTgctg                                                                                      >  2:487216/1‑299 (MQ=255)
                                                                                                                                                                                                                          tttttggacggtgggggggCGTATTGGGCGCGAGTGAATCTGGCCTGGCAGGATGTCTTTATGATGCCGGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAAGTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTTTGCGGCGGaa                                                                                 <  1:17639/285‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ggtggcgtgATGTCTTTTTGGTGCCTGCGCGTCGTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTGAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCGGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTGTGCGGCGGAACCACAATCCCCAAAAAAGTGGCGCGTGAATGCCAGCAGCGCgg                                      <  1:494246/293‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTTTGCGGCGGAACCACAATCCCCAACAAAGTGGCGCGTGAATGCCAGCAGCGCGGCa                                    >  2:557677/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                         gCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGGGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACATTTCAGCGCTGCGTTTCTTTCTTTGCGGCGGAACCACAATCCCCAAAAAAGTGNCGCGTTAATGACAGCAACACGGCAtt                                  >  2:304223/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                             taTGTCTTTATGATGCTTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTTTGCGGCGGAACCACAATCCCCAAAAAAGTGGCGCGTGAATGCCAGCAGCGCGGCATTAAat                              <  2:50087/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tgtctgttggatgCTTGTGCCTCTTTGTCACGCAACGGTCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTTTGCGGCGGAACCACAATCCCCAAAAAAGTGGCGCGTGAATGCCAGCAGCGCGGCATTAAatta                            <  2:12793/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tGTCTTTATGATGCCTGCGCCACTTGTTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTTTGCGGCGGAACCACAATCCCCAAAAAAGTGGCGCGTGAATGccgcgctgctggcattcacgcg                            <  2:54292/301‑21 (MQ=255)
                                                                                                                                                                                                                                                                               tGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTTTGCGGAGGAACCACAATCCCCAAAAAAGTGGCGCGTGAATGCCAGCACCCCGGCATTAAattt                            >  1:12793/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                           gCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTTTGCGGCGGAACCACAATCCCCAAAAAAGTGGCGCGTGAATGCCAGCAGCGCGGCATTAAATTATTAAGTGTTTAt                >  2:44828/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                         gTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGCGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCTCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGACACAACCCGCGGACATTTCAGCGCTGCGTTTCTTTCTTTGCGGCCGACCCACAATCCCCAAAAACGTGGCGCGAGCATGCCAGCAGCGCGGCATTACATTAACAAGTGTTTATCGATCCACAaaaaa  >  1:365810/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         
TAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCACCTCTTCCCTCTCATTAAGTCAGATTATCGCCGACAATACCTCACTGACCACGGCGATAACGACCCACGGCGATGAATTAGCTGCGGTGCTGTTTACCTCCGGAACCGAGGGTCTGCCAAAGGGCGTGATGCTAACGCATAACAATATTCTCGCCAGTGAGCGGGCTTATTGCGCGCGACTGAATCTGACCTGGCAGGATGTCTTTATGATGCCTGCGCCACTTGGTCACGCAACGGGCTTTCTGCATGGCGTAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCTTTTGAATGTACTAGAGAAACAACCCGCGGACCTTTCAGCGCTGCGTTTCTTTCTTTGCGGCGGAACCACAATCCCCAAAAAAGTGGCGCGTGAATGCCAGCAGCGCGGCATTAAATTATTAAGTGTTTATGGTTCCACAGAAAG  >  NC_000913/1783449‑1784046

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: