Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,750,078 T→G 5.6% intergenic (+18/+37) yfjD → / ← grpE UPF0053 family inner membrane protein/heat shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,750,0780TG 5.6% 242.0 / 4.8 72intergenic (+18/+37)yfjD/grpEUPF0053 family inner membrane protein/heat shock protein
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base G (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01

CGACGGTTCGGTGATTATCGATGGCACCGCCAACGTGCGGGAAATCAACAAAGCCTTTAACTGGCATCTACCGGAAGATGATGCCCGCACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGC  >  NC_000913/2749808‑2750344
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                          
cgacgGTTCGTGGATTATCGATGGACCCGCCAACGTGCGGGAAATCAACAAAGCCTTTAACTGGCATCTACCGGAAGATGATGCCCGCACGGTTAATGGCGTCATGCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTc                                                                                                                                                                                                                                                                           <  4:480443/272‑1 (MQ=255)
    ggttcgggtttgttCGATGGCACCGCCAACTTTCGGGAAACCAACAAAGCCTTTAACTGGCATCTACTGGAAGTTGATGCCCGACCTGTTAACGGCGTCATTCTTGAGGCACTGGAAGCGTTCCCTGTCGCAGGAACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGTAAAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCa                                                                                                                                                                                                                                                                    <  4:528056/262‑1 (MQ=255)
    ggTTCGGTGATTATCGATGGCACCGCCAACGAGCGGGAAATCAACAAAGCCTTTAACGGGCATCTACCGGAAGATGATGCCCGAACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCa                                                                                                                                                                                                                                                                    <  2:249637/275‑1 (MQ=255)
                       gCACCGCCAACGTGCGGGAAATCAACAAAGCCTTTAACTGGCATCTACCGGAAGATGATGCCCGCACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGa                                                                                                                                                                                                                                                 <  2:248932/275‑1 (MQ=255)
                                                caaAGCCTTTAACTGGCATCTACCGGAAGATGATGCCCGCACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCNCGCCGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGGNGAGTAACGAAAACGGCCCGGCATCCGCAGCCCGGGCCGGGAAGTANTNAGAAAGCAGAAATAAACCttttagttta                                                                                                                                                                                                                        >  3:42407/1‑269 (MQ=255)
                                                               gCATCTACCGGAAGATGATGCCCGCACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCATGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTAAAGTAACCAtt                                                                                                                                                                                                         >  1:87323/1‑274 (MQ=255)
                                                                         gAAGATGATGCCCGCACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTTGCTACCGTAACCATCGCCGCACGaa                                                                                                                                                                                               >  2:476194/1‑275 (MQ=255)
                                                                                       cACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCAtt                                                                                                                                                                                  >  4:333571/1‑274 (MQ=255)
                                                                                                    gTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCa                                                                                                                                                                                                           >  4:341085/1‑236 (MQ=255)
                                                                                                     tCATTCTTAAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCtt                                                                                                                                                                   <  2:119834/275‑1 (MQ=255)
                                                                                                      cATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTc                                                                                                                                                                  >  2:462587/1‑275 (MQ=255)
                                                                                                      cATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTCGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTc                                                                                                                                                                  >  2:495596/1‑275 (MQ=255)
                                                                                                        ttCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATcgc                                                                                                                                                                                                       <  2:341086/236‑1 (MQ=255)
                                                                                                            tGAGGAACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTGCGACCATTCAGCGTATAACCCTTCTGCATa                                                                                                                                                            <  2:531794/275‑1 (MQ=255)
                                                                                                            ggAGGCACTTGAAGAGTGCCCTGTCGCAGGCACCGGAGTGCGTATTGGTGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTGTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATa                                                                                                                                                            <  4:124950/274‑1 (MQ=255)
                                                                                                             gAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATa                                                                                                                                                            >  4:64336/1‑274 (MQ=255)
                                                                                                                  aCTGGAAGAGATCCCTGACGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGCCGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGcc                                                                                                                                                       <  1:150800/274‑1 (MQ=255)
                                                                                                                   cTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGTGCCGTGAATTATTACGAAAGCAGACATTAAGCTTTTGCTTTCGCTACAGTAACCATCTCCGCACGAAGCCTACAACCTTTCAGCGACTAAACCTTCTGCATacatccca                                                                                                                                                     >  4:529491/1‑268 (MQ=255)
                                                                                                                   cTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCt                                                                                                                                                                                                                           >  4:199488/1‑205 (MQ=255)
                                                                                                                   cTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCa                                                                                                                                                     >  2:415024/1‑275 (MQ=255)
                                                                                                                       aaGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCg                                                                                                                                                                                                                       <  2:199489/205‑1 (MQ=255)
                                                                                                                              accggTCCCAGGGACCGGGGGGCGTATTGGCGAGTACGATATCGATATTCTCGACGTAAAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTAcc                                                                                                                                           <  4:415023/270‑1 (MQ=255)
                                                                                                                                 tGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCCGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTgga                                                                                                                                       >  2:432579/1‑274 (MQ=255)
                                                                                                                                     gCAGGCACCCGCGTGCGTATTGGCGAGTACAATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAGCCGCTGCGCGAGAGTGTGGCGGAGTAACGAAATCGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTc                                                                                                                                                                                 <  2:63505/229‑1 (MQ=255)
                                                                                                                                            cccGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCgg                                                                                                                                                                                                                                                                  >  7:185290/1‑141 (MQ=255)
                                                                                                                                            cccGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCgg                                                                                                                                                                                                                                                                  <  8:185290/141‑1 (MQ=255)
                                                                                                                                             ccGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAAc                                                                                                                                                                       >  4:73965/1‑231 (MQ=255)
                                                                                                                                               atgtgcgtATTGGGGATAACAAGTACGATATTCTCGACGGAAAGGACAATATGATTAAGAAGGTAAAAGTTGTTCCTGGGAACCCGCTGCTGGGGGGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTGTTGCTTTCGCTACAGTCACCATCGCCGCACGAATCGTACGACCATTCAGCGTATCACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGa                                                                                                                         <  4:267111/273‑1 (MQ=255)
                                                                                                                                                 gtgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCtt                                                                                                                                                                   <  2:73966/231‑1 (MQ=255)
                                                                                                                                                   ggtgaTTGGGGAGTAGGATATGGATATTGTCGACGTAGAGGAAAATATGATTAAGAAGGTAAAAGTTTTGCCTGTGAAACCGCTGCGAGAGAGTGTGGCGGAGTAACGAAAACGGCCGGGCATTAGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCGTTCGCTACAGTAACCATCGGCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTc                                                                                                                      <  2:356821/270‑1 (MQ=255)
                                                                                                                                                            cGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCg                                                                                                            <  2:333572/275‑1 (MQ=255)
                                                                                                                                                            cGAGTACGATATCGATATTCTCGACGTAAAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCg                                                                                                            <  1:42407/275‑1 (MQ=255)
                                                                                                                                                                taggaaatcgttatTCTCGACGTACCGGACAAGATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGAATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCGTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATg                                                                                                         <  4:476193/268‑1 (MQ=255)
                                                                                                                                                                tACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCCCAACGTCATCAGATTCCACCATTTCAATgg                                                                                                        >  4:226627/1‑275 (MQ=255)
                                                                                                                                                                                tcGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCgg                                                                                        <  2:226628/275‑1 (MQ=255)
                                                                                                                                                                                  gACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGt                                                                                      >  4:100225/1‑275 (MQ=255)
                                                                                                                                                                                           gACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCg                                                                                                                                                                                                                                                                          <  8:36757/86‑1 (MQ=255)
                                                                                                                                                                                           gACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCg                                                                                                                                                                                                                                                                          >  7:36757/1‑86 (MQ=255)
                                                                                                                                                                                            aCAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGc                                                                                                                                                                                                                                                                >  8:618113/1‑95 (MQ=255)
                                                                                                                                                                                            aCAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGc                                                                                                                                                                                                                                                                <  7:618113/95‑1 (MQ=255)
                                                                                                                                                                                            aCAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAatt                                                                                                                                                                                                                                                        >  6:113801/1‑103 (MQ=255)
                                                                                                                                                                                            aCAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAatt                                                                                                                                                                                                                                                        <  5:113801/103‑1 (MQ=255)
                                                                                                                                                                                              aaTATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTcc                                                                                                                     >  4:25808/1‑232 (MQ=255)
                                                                                                                                                                                                tatGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAATGTGGCGGAGTAACGAAAACGGCGCGGCATTCGCATGCAGGGCCGTGAATTATTACGAATGCAGAAATTAAGCTTTTGCTTTCTCTACAGTAACCCTCGCCGCACGACTCGTACGACCATTCAGCGTATCACCCTTCTGCATAATGCCCAGTCCGTTACCTGGCGCCACGTCATTCGATTCCACCATTTCGATGGCCTGATGCACCTTCGGGTCCAGTgggcatt                                                                          >  4:96339/1‑269 (MQ=255)
                                                                                                                                                                                                  tGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTccacca                                                                                                                 <  2:25809/232‑1 (MQ=255)
                                                                                                                                                                                                    aTTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAg                                                                                                                                                                                                                                                                   <  8:570647/84‑1 (MQ=255)
                                                                                                                                                                                                    aTTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAg                                                                                                                                                                                                                                                                   >  7:570647/1‑84 (MQ=255)
                                                                                                                                                                                                     ttAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGaa                                                                                                                                                                                                                                                <  8:884713/102‑1 (MQ=255)
                                                                                                                                                                                                     ttAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGaa                                                                                                                                                                                                                                                >  7:884713/1‑102 (MQ=255)
                                                                                                                                                                                                     ttAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCNTCGCCGCACGCATCGTACGACCATTCAGCGTATCACCCTTCTGCTTCATCCCCAGTCCTTTCCCTGGCGCAACCNCATCAGATTCCCCCCTTCCGATGGCCTGTTNNACATTCTCTTCCAGCGTGACGTTAttt                                                                     >  3:62165/1‑273 (MQ=255)
                                                                                                                                                                                                     ttAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTTAGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCAACTTCGGGTCCAGGGGGACGTTAGTTTc                                                                   >  4:348891/1‑275 (MQ=255)
                                                                                                                                                                                                     ttAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGttt                                                                    >  1:101727/1‑274 (MQ=255)
                                                                                                                                                                                                        aGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGTTTCgg                                                                 >  2:5974/1‑274 (MQ=255)
                                                                                                                                                                                                          cAGGTAAAAGGTTTTCCTGTGAAACCGCTGCGCGAGAGGGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGTTTCGGCg                                                               <  2:348892/274‑1 (MQ=255)
                                                                                                                                                                                                             gTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACg                                                                                                                                 >  1:26479/1‑205 (MQ=255)
                                                                                                                                                                                                                 aaGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGtcat                                                                                                                             <  3:26479/205‑1 (MQ=255)
                                                                                                                                                                                                                     ttttCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGccgg                                                                                                                                                                                                                                                                  >  8:1063084/1‑65 (MQ=37)
                                                                                                                                                                                                                     ttttCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCgg                                                                                                                                                                                                                                                                  <  7:1063084/68‑1 (MQ=37)
                                                                                                                                                                                                                       ttCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCATGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGCTTCCACCATTGCGATGGCCTGCTGCACATTCGGGTCCAGTGGGACGTTAGTTTCGGCGATCACTTTCACGc                                                  >  2:107788/1‑274 (MQ=255)
                                                                                                                                                                                                                                  ccgtgtgcggggagTGTGGAGTAGAGACANCAACGGCCCGTCATTCGCATGCAGGNCCGTGAATTATTACGAAAGCAGAAATTAATCTTTTGCTTTCGCTACAGTAACCATCGCCTCACGAATCGTACGACCATTCNGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGTTTCGGCGATCACTTCAACGCCAAACTTACGca                                      <  3:101727/264‑1 (MQ=255)
                                                                                                                                                                                                                                  acctggggggggagTTTGGCGGAGGAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAAGTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGTTTCGGCGATCACTTCAACGCCAAACTTACGca                                      <  4:44338/264‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCa                                                                                   >  4:385879/1‑227 (MQ=255)
                                                                                                                                                                                                                                         gcgAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTgg                                                                               <  2:385880/227‑1 (MQ=255)
                                                                                                                                                                                                                                                 gtgGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTa                                                                                                                                                                                                                                     >  7:601740/1‑69 (MQ=255)
                                                                                                                                                                                                                                                 gtgGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTa                                                                                                                                                                                                                                     <  8:601740/69‑1 (MQ=255)
                                                                                                                                                                                                                                                  tgGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGCTCCAGTGGGACGTTAGTTGCGGCGATCACTTCCACNCCAAANNCACGCACAACAACCAACAACg                       >  3:14690/1‑274 (MQ=255)
                                                                                                                                                                                                                                                     cGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCttt                                                                                                                                                                                                                               >  5:312379/1‑71 (MQ=255)
                                                                                                                                                                                                                                                     cGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCttt                                                                                                                                                                                                                               <  6:312379/71‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGt             <  2:100226/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cccGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGc  >  2:538954/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                        cGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTc                                                                                                                                                                  >  6:419036/1‑113 (MQ=255)
                                                                                                                                                                                                                                                                        cGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTc                                                                                                                                                                  <  5:419036/113‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                          
CGACGGTTCGGTGATTATCGATGGCACCGCCAACGTGCGGGAAATCAACAAAGCCTTTAACTGGCATCTACCGGAAGATGATGCCCGCACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGC  >  NC_000913/2749808‑2750344

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