Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,085,928 G→A intergenic (‑106/‑88) fdoG ← / → fdhD formate dehydrogenase O subunit alpha/sulfurtransferase for molybdenum cofactor sulfuration

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,085,9280GA97.6% 144.9 / NA 42intergenic (‑106/‑88)fdoG/fdhDformate dehydrogenase O subunit alpha/sulfurtransferase for molybdenum cofactor sulfuration
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/25);  minor base C (1/0);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCGCTGGTC  >  NC_000913/4085684‑4086153
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 
gggTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCACAATTgg                                                                                                                                                                                                                             >  1:378327/1‑251 (MQ=255)
 ggTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGc                                                                                                                                                                                                                            <  2:330004/251‑1 (MQ=255)
        cgGGTGCGTAGCAGTTTATACTTCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCGGTGCCGCTGCCGTGGGGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATGCGCTTGCTGCGTGTGCAGGTAGTGGTAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAAttt                                                                                                                                                                                                                     <  1:259741/251‑1 (MQ=255)
         gggTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTc                                                                                                                                                                                                                    <  2:316884/251‑1 (MQ=255)
             gcgcAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCtc                                                                                                                                                                                                                <  2:307813/251‑1 (MQ=255)
              cgcAGCAGTTTATACTGCCGGGTCGCCGCGAGTCCTACGCGGGGTGCAAACCCCAGTGCCGCTGCCGTGGGGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGATATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCAGCAAAGGCCCAAAATTGGCTCAATTTCCCtcc                                                                                                                                                                                                               <  2:62806/251‑1 (MQ=255)
                       ttATACTGCCGGGTTTCCGCGAGTTCTACGCTGGGTGCAAAACCCAGTGCGGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGa                                                                                                                                                                                                       <  2:188359/250‑1 (MQ=255)
                             tGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAgg                                                                                                                                                                                                    >  2:336060/1‑247 (MQ=255)
                             tGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAgg                                                                                                                                                                                                    <  1:336060/247‑1 (MQ=255)
                              gCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGtt                                                                                                                                                                                                 <  1:82493/249‑1 (MQ=255)
                                                cTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATg                                                                                                                                                                             >  2:227469/1‑251 (MQ=255)
                                                                         cGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGtt                                                                                                                                                    <  1:326809/251‑1 (MQ=255)
                                                                         cGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGtt                                                                                                                                                    <  1:41352/251‑1 (MQ=255)
                                                                         cGCTGCCGTGGTGCCCGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTTCATGGATTGCTCCTTGGTTCGACATTGTCACGTTCCATTTTCATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTTCATCCAAGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGtt                                                                                                                                                    <  1:79391/251‑1 (MQ=255)
                                                                               cGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAgaga                                                                                                                                              <  1:223400/251‑1 (MQ=255)
                                                                                  ggtgCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGaaa                                                                                                                                            <  2:320472/250‑1 (MQ=255)
                                                                                    ggcctgccATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTTTTCTTATCTGGTTAAGAGAAAGTg                                                                                                                                         <  2:371344/250‑1 (MQ=255)
                                                                                       ctgccATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGaa                                                                                                                                       >  1:323222/1‑250 (MQ=255)
                                                                                         gccATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGaaa                                                                                                                                      <  2:147826/249‑1 (MQ=255)
                                                                                         gccATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAAGACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGaaa                                                                                                                                      <  1:225904/249‑1 (MQ=255)
                                                                                                cGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGCGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGCAAAAAccac                                                                                                                               >  1:256751/1‑249 (MQ=255)
                                                                                                               tAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAAt                                                                                                                  <  1:342077/247‑1 (MQ=255)
                                                                                                               tAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATc                                                                                                                 <  2:323222/248‑1 (MQ=255)
                                                                                                                aaaGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGaa                                                                                                              >  1:384710/1‑250 (MQ=255)
                                                                                                                    aaCTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGaaa                                                                                                             <  1:392037/247‑1 (MQ=255)
                                                                                                                    aaCTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGaaa                                                                                                             >  2:392037/1‑247 (MQ=255)
                                                                                                                        gCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGa                                                                                                       <  1:239611/249‑1 (MQ=255)
                                                                                                                            tCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAat                                                                                                  >  1:338845/1‑250 (MQ=255)
                                                                                                                                         aTGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCg                                                                                    <  2:338845/251‑1 (MQ=255)
                                                                                                                                                      tGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGtt                                                                        >  2:136719/1‑250 (MQ=255)
                                                                                                                                                      tGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGt                                                                         <  2:378327/249‑1 (MQ=255)
                                                                                                                                                           cGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATggcg                                                                  <  1:5813/251‑1 (MQ=255)
                                                                                                                                                                 tGTCACGTCCCCTTTACATTCGATTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGCAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGCTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTCTATGGTTAAGAGAAAGGGAAAAAACTACAGCGCACAGAAATCGAAAATGTGACAACTATCAAAGGTGTTCGTCAACTCGAGTTATGGCGGCGtgg                                                            >  2:78551/1‑250 (MQ=255)
                                                                                                                                                                  gTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGtgat                                                           >  2:176085/1‑251 (MQ=255)
                                                                                                                                                                            aTTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTacaac                                                 <  1:176085/251‑1 (MQ=255)
                                                                                                                                                                                  aTTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGCGATGATTTACAACACCCAc                                           >  1:78838/1‑251 (MQ=255)
                                                                                                                                                                                    tCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAAAACCCCCgg                                         >  2:208390/1‑251 (MQ=255)
                                                                                                                                                                                          gctgcGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGtg                                                             <  2:316371/225‑1 (MQ=255)
                                                                                                                                                                                          gctgcGCGAGCAGGGAGTGGGAGTTATTTTTCTTTGAGGAAGGGGCCGCAAAGGACCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGCAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGtg                                                             >  1:316371/1‑225 (MQ=255)
                                                                                                                                                                                                tgtgCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTc                             >  1:153087/1‑251 (MQ=255)
                                                                                                                                                                                                                    tttttCTTTGCGGAAGGGCCCGCAAAGGGCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGGGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTgc         >  1:115298/1‑251 (MQ=255)
                                                                                                                                                                                                                           ttGCGGAAGGGGCCGCAAAGGTCCAAAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGAGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACAACCACGGGTAGATGAGGGCGCGGAAGAAGTTCCCGCTGCTCTGGTc  >  2:92707/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 
GGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACAGAATACCATAATGTTGGTGTGTGTGTTCTTATCTGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCGCTGGTC  >  NC_000913/4085684‑4086153

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