Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,319,404 A→C 22.9% V68G (GTG→GGG)  phnJ ← carbon‑phosphorus lyase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,319,4040AC22.9% 116.4 / 2.1 57V68G (GTG→GGG) phnJcarbon‑phosphorus lyase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (23/17);  new base C (2/14);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 2.71e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCTCAGATTAGCCATGGTTCTGCTCCTGTTGCAGACGTTTGAGT  >  NC_000913/4319190‑4319635
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       
cGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCa                                                                                                                                                                                                                  >  1:592313/1‑238 (MQ=255)
            cTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATATCCGTCGCATCGTCCGTACGTTCCGGGGTGTTTACCCCGGTCACGCGCTTGAAGAACTTGCCAATCGACACGGCGTTGGTGTGCTCGTCCGCGCCCTGGTCAAACCCCTTCACCCCTCCGCTTTCGCCAATCaacccgccgg                                                                                                                                                                                                     >  2:740126/1‑230 (MQ=255)
                      cgcAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTGGCCAATTACGCTGGCGGTAAGCCGAAt                                                                                                                                                                                           >  1:304574/1‑239 (MQ=255)
                         aGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGGTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCCCCTTCAGCACTGCGCTTTCGCCAATCACGCTGGGGGTAAGCGGAAGgcc                                                                                                                                                                                        >  2:443833/1‑239 (MQ=255)
                                    gATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTccc                                                                                                                                                                              <  1:658361/238‑1 (MQ=255)
                                        ggCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCCTCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATACCCCTGACGGTAAGCTGAAGGCCGCCGGTTCCCCacc                                                                                                                                                                          >  1:526470/1‑236 (MQ=255)
                                                               aTCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATgcc                                                                                                                                                                                        <  1:238536/201‑1 (MQ=255)
                                                               aTCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATgcc                                                                                                                                                                                        >  2:238536/1‑201 (MQ=255)
                                                                         ggcggcgGTTGGGGGTTGCGGTGCCGGGTCTGGTTACCCGTCGCATCGTCCGTACCTTCCGTGGTGTTTTCCCCGGTCACGCGCTTGAAGAGGTTGCGAATCGCCACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCCCCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACgg                                                                                                                                        <  2:592313/238‑1 (MQ=255)
                                                                         cgcggcgTTTGGGGGTTCGGGGTGCGGGTCGGGATACCCTCCGCGTCGGCCGGACGTTCCGGGTTGTTTACCCCGGTCACGCGCTGGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCGCCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACgg                                                                                                                                        <  2:931957/238‑1 (MQ=255)
                                                                          gcgggctttCGGGGATGCGGGGACGCGGTTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACgg                                                                                                                                        <  1:555596/232‑1 (MQ=255)
                                                                            ggcgTTTCGGGGATGCGGTGACGGGTCTGGATAACCGGCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACgg                                                                                                                                        <  1:416053/236‑1 (MQ=255)
                                                                               gTTTGGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTTTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCt                                                                                                                                   <  1:740126/238‑1 (MQ=255)
                                                                                tttCGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCCCCTGa                                                                                                                                 >  1:204369/1‑239 (MQ=255)
                                                                                        aTGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCCTCGGCATCTCGCGCCCGCCAAACGGCACCCGTTAACCggg                                                                                                                         >  1:806478/1‑239 (MQ=255)
                                                                                        aTGCGGTGACGCGTCTGGATAACCGTCACATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCCCCCAAAGGCACCTGGTAACCcgg                                                                                                                         >  1:798184/1‑239 (MQ=255)
                                                                                           gCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAg                                                                                                                                                                           >  1:260742/1‑187 (MQ=255)
                                                                                           gCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAg                                                                                                                                                                           <  2:260742/187‑1 (MQ=255)
                                                                                            cGGTGACGCGTCTGGATAACCGGCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGcgcccg                                                                                                                                                <  1:603163/213‑1 (MQ=255)
                                                                                                 gcgctttcgtaTAACCGTGCCTTTGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCa                                                                                                                 <  2:204369/229‑1 (MQ=255)
                                                                                                  cgcgTCTGGATTACCTTTCGCTCGTTCGTTCGTTCCGGTGGGTTTGCCCGGGCCCCGGCCTTGAAGAAGTTGCAGATCGACCCGGCGTTGGTGTGTTCATCCGGGCCCTGGTCAATCCCCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCAc                                                                                                                <  2:284053/239‑1 (MQ=255)
                                                                                                   gcgtttgcgtAACCCTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCCCCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCAc                                                                                                                <  1:992438/229‑1 (MQ=255)
                                                                                                   gcgttgggATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCCCGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCAc                                                                                                                <  1:64689/232‑1 (MQ=255)
                                                                                                   gcgTTTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCAcc                                                                                                               <  1:671144/239‑1 (MQ=255)
                                                                                                               ccGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCCTCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCCCCGCTTTTAAGCTg                                                                                                   >  1:335030/1‑239 (MQ=255)
                                                                                                               ccGTCGCATCGTCCGTACGTTCCGTGGGGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGGACGGCATCGCCATCCCGCCCCCCCCAAACGGCACCTGATAACCGGGGGACCCCCCCCCTTTTAAGagg                                                                                                   >  2:923705/1‑237 (MQ=255)
                                                                                                                        tCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCCCCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCgcggcg                                                                                           <  1:507759/239‑1 (MQ=255)
                                                                                                                          gTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACAGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGAGGCGGCGCCGGa                                                                                         >  1:549519/1‑239 (MQ=255)
                                                                                                                            ccGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGAt                                                                                        <  1:800050/238‑1 (MQ=255)
                                                                                                                             ctttctttCCGGGGTGTTTGCCCCGGTCCCGCCCTTGAAGAATTTGCGAATCGACCCGGCGTTGGGGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCCCGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCa                                                                                      <  2:761870/233‑1 (MQ=255)
                                                                                                                                      cGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTCCGGCATCGGCATCTCGCGCCCGCCAAACGGCCCCTGTTAACCCGGGATCGCCCCCCCTTTTAAGATGGCGCGGCGGGTCATGCGTTTgg                                                                             >  1:453897/1‑239 (MQ=255)
                                                                                                                                                cccccGGTCACGCTGGTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCCCCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCtcgtc                                                                   <  1:850494/238‑1 (MQ=255)
                                                                                                                                                     ggTCACGCGCTTGAAGAAGTTGCGAATCGACCCGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCCCCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGt                                                              <  1:929645/239‑1 (MQ=255)
                                                                                                                                                     ggTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGCTCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGGGGt                                                              >  1:617709/1‑239 (MQ=255)
                                                                                                                                                     ggTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCCCGTCGAGGt                                                              >  1:945204/1‑239 (MQ=255)
                                                                                                                                                     ggTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGCCATCTCGGCCCCGCCAAACGGCACCTGATAACCCGGGATGACACCCGCTTTTAAGATTGCGCGGGGGATCACGCGTTTTGGCTGCCCGCCggggg                                                              >  2:638513/1‑238 (MQ=255)
                                                                                                                                                      gTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAg                                                                                                      >  1:888057/1‑198 (MQ=255)
                                                                                                                                                                   agaagTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAg                                                <  1:843481/239‑1 (MQ=255)
                                                                                                                                                                         tcccaatcctcacGGGCTTGGGTGTTTCCTCCGCGCCCTGGGCAATCACCTTCGCCACGGCGCTTTCGCCAATCCCCCTGGCGGTAAGCTGAATGCCCCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCt                                           <  2:798184/228‑1 (MQ=255)
                                                                                                                                                                         tGCGCATCTACCCGGGGTTGGTGGGTCCACCCGCGCCCTGGTCAACCCCCGTCAGCCCGTCGCTTTCGCCACTCCCGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGGCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCt                                           <  2:860273/238‑1 (MQ=255)
                                                                                                                                                                         tGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATgcc                                                                                                                                                                                        <  1:849022/97‑1 (MQ=255)
                                                                                                                                                                         tGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATgcc                                                                                                                                                                                        >  2:849022/1‑97 (MQ=255)
                                                                                                                                                                         tGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGt                                                              >  1:709270/1‑219 (MQ=255)
                                                                                                                                                                         tGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGt                                                              <  2:709270/219‑1 (MQ=255)
                                                                                                                                                                         tGCGAACCGGCACGGCGTTGGTGGTGTCAACCGCGCCATGGTTAATCCCCTTCAGCCCGTGGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCt                                           <  1:980606/238‑1 (MQ=255)
                                                                                                                                                                              ccgcgcccGGGGTTGGTTGTGTCCTCCGCGCCCTGGTCAAGCCCCTTCAGCACGTCGCTTTCGCCATTCACGCTGGCGGGAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGAATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCAGCTCAGAtt                                     <  2:453897/232‑1 (MQ=255)
                                                                                                                                                                                  acacGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGGTCATGCGATTGGTCTGCTCGTCGAGGTAGGCCAAGGTGTAGACGCCCAGCTTAGcc                                 >  1:240677/1‑239 (MQ=255)
                                                                                                                                                                                   cacGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCCTCGGCCTCTCGCGCCCGCCAAACGGCACCTGCTAACCGGGGGTCGCCACCGCCTTTAAGATGGCGCGCGGGACATTGCGTTTGGTCTGCTCGTCGCGGGGGGGAAAGTTGGTGCCCCCCCAATTAGCCa                                >  2:528727/1‑239 (MQ=255)
                                                                                                                                                                                      ggCGTTGGTGGTGTCATCCGAGCCCTGGTTAATCCCCTTCAGCACGTCGCTTTAGCCAATCACGCGGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCTCAGATTAGCCATgg                             <  1:107131/239‑1 (MQ=255)
                                                                                                                                                                                                gtgtCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCCTTCGTTTGGTCTGGCCGTCGGGGTAGGCAAAGTTGTAGCCGCCCAGATTAGCCCTGGTTCTGGTCCt                   >  1:706552/1‑239 (MQ=255)
                                                                                                                                                                                                  gtCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGGCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCCCAGATTAGCCATGGTTCtttccccgt                 >  1:375272/1‑231 (MQ=255)
                                                                                                                                                                                                      tCCGCCCCCTGGTCAATCCCCTTCACCCCTTCGCTTTCGCCAATCCCGGTGGCGGTAAGCTTAATGCCGCCTGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCTCAGATTAGCCATGGTTCTGCTCCTGTTGCa             <  1:509183/239‑1 (MQ=255)
                                                                                                                                                                                                      gccccgcCCTGGTCAAGCCCCTTCAGCACGGCCCTGTCGCCAATCGCGCTGGGGGTAAGCTGAATGCCGCCGGTACCCTAGCCGGCCGGCATCGGCATCGCGCGCCCGCCAAACGGCCCCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCTCAGATTAGCCATGGTTCTGCTCCTGTTGCa             <  2:806478/235‑1 (MQ=255)
                                                                                                                                                                                                        cgcgCCCTGGTCAATCCCCTTCAGCCCGGCGCTTTCGCCAATCCCGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCTCAGATTAGCCATGGTTCTGCTCCTGTTGCAGa           <  1:528727/239‑1 (MQ=255)
                                                                                                                                                                                                          cgcCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCCCCTTTAAGATGGCGCGGCGGATCATGCGTTTGGGCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCCCCCGGTTTAGCCATGGGTCTGCTCCTGTTGCAGACg         >  1:269087/1‑239 (MQ=255)
                                                                                                                                                                                                                ggTCAATCACCTTCAGCACGTCGCTTTCGCCAATCCCGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCTCAGATTAGCCATGGTTCTGCTCCTGTTGCAGACGTTTGAg   <  1:970007/239‑1 (MQ=255)
                                                                                                                                                                                                                  tCATTCCCCTTCAGCCCGTCGCTTTCGCCAATCGCGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCTCAGATTAGCCATGGTTCTGCTCCTGTTGCAGACGTTTGAGt  <  1:638513/238‑1 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       
CGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTCGCATCGTCCGTACGTTCCGTGGTGTTTACCCCGGTCACGCGCTTGAAGAAGTTGCGAATCGACACGGCGTTGGTGGTGTCATCCGCGCCCTGGTCAATCACCTTCAGCACGTCGCTTTCGCCAATCACGCTGGCGGTAAGCTGAATGCCGCCGGTTCCCCAGCCGTACGGCATCGGCATCTCGCGCCCGCCAAACGGCACCTGATAACCGGGGATCGCCACCGCTTTTAAGATGGCGCGGCGGATCATGCGTTTGGTCTGCTCGTCGAGGTAGGCAAAGTTGTAGCCGCTCAGATTAGCCATGGTTCTGCTCCTGTTGCAGACGTTTGAGT  >  NC_000913/4319190‑4319635

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