Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,750,073 G→T 5.4% intergenic (+13/+42) 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,0730GT 5.4% 256.9 / 4.3 74intergenic (+13/+42)yfjD/grpEUPF0053 family inner membrane protein/heat shock protein
Reads supporting (aligned to +/- strand):  ref base G (36/34);  new base T (2/2);  total (38/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 = 8.30e-01

CAAAACGACGGTTCGGTGATTATCGATGGCACCGCCAACGTGCGGGAAATCAACAAAGCCTTTAACTGGCATCTACCGGAAGATGATGCCCGCACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGC  >  NC_000913/2749803‑2750344
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                               
cAAAACGACGGTTCGGTGATTATCGATGGCACCGCCAACGTGCGGGAAATCAACAAAGCCTTTAACTGGCATCTACCGGAAGATGATGCCCGCACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAGGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGNNNGTTGCGGAGAAAAGAAAACGGCCCGGCAt                                                                                                                                                                                                                                                                             >  3:157268/1‑275 (MQ=255)
 aaaaCGACGGTTCGGTGATTATCGATGGCACCGCCAACGTGCGGGAAATCAACAAAGCCTTTAACTGGCATCTACCGGAAGATGATGCCCGCACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGc                                                                                                                                                                                                                                                                               >  2:480444/1‑272 (MQ=255)
     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                                                                                                                                                                                                                                                                  <  8:185290/141‑1 (MQ=255)
                                                                                                                                                 cccGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCgg                                                                                                                                                                                                                                                                  >  7:185290/1‑141 (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                                                                                                                                                                                                                               <  6:312379/71‑1 (MQ=255)
                                                                                                                                                                                                                                                          cGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCttt                                                                                                                                                                                                                               >  5:312379/1‑71 (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)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                               
CAAAACGACGGTTCGGTGATTATCGATGGCACCGCCAACGTGCGGGAAATCAACAAAGCCTTTAACTGGCATCTACCGGAAGATGATGCCCGCACGGTTAATGGCGTCATTCTTGAGGCACTGGAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCATTGCGATGGCCTGATGCACATTCGGGTCCAGTGGGACGTTAGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGC  >  NC_000913/2749803‑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: