Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,325,837 C→T intergenic (‑96/+562) yjdN ← / ← yjdM PF06983 family protein YjdN/zinc ribbon domain‑containing protein YjdM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,325,8370CT89.8% 119.0 / NA 51intergenic (‑96/+562)yjdN/yjdMPF06983 family protein YjdN/zinc ribbon domain‑containing protein YjdM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/16);  minor base G (3/2);  total (33/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAAC  >  NC_000913/4325601‑4326068
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                       
gTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGCTGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACAGTTg                                                                                                                                                                                                                           >  2:137558/1‑251 (MQ=255)
     cGGCGCTGGCCTTCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTTCGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGATGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGGGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCgg                                                                                                                                                                                                                       <  2:113560/250‑1 (MQ=255)
        cgcTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGGCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATgc                                                                                                                                                                                                                   <  1:55365/251‑1 (MQ=255)
         gcTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATCATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGCCTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGGAGTCGGCACTGTTGCCGGATgcg                                                                                                                                                                                                                  >  2:286919/1‑251 (MQ=255)
             tCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGt                                                                                                                                                                                                              >  1:57431/1‑251 (MQ=255)
                 gcgcTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGCTAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTaagc                                                                                                                                                                                                          >  2:104283/1‑249 (MQ=255)
                           tGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                                                                 <  2:76122/250‑1 (MQ=255)
                                     ccGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACgg                                                                                                                                                                                        >  1:302313/1‑249 (MQ=255)
                                       gAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGGAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggg                                                                                                                                                                                    >  1:272496/1‑251 (MQ=255)
                                       gAAGTTGATTTTATAGAGCAGTTCCGCGCACCACGTACGTTGATAATAGGCAATCGCGACGGAACCGTTACCGGCAAAAGAGAGGTAGGGCCTTAACGGCATGCTGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGGGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAAAGCCTTATCCGGCCTACGGAggg                                                                                                                                                                                    >  2:361771/1‑251 (MQ=255)
                                              aTTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                                                                                             <  1:76394/251‑1 (MQ=255)
                                                 ttATAGAGCAGGTCCGCGCCCAACGTACGTGGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCTAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGTGTGCGGGAAt                                                                                                                                                                           <  1:265760/250‑1 (MQ=255)
                                                        gCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGATTTGTAGGc                                                                                                                                                                   >  2:288171/1‑251 (MQ=255)
                                                           gatCCGCGCCCACCGTACGTTGATCATAGGCAATCGCGTCGGAACAGTTACCGGCACAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGAAGTTCGAACTTTATACAGATGAGAGGCGGCCTGATAAGACGCGAGAGCGTCGCATCAGGGAGTCGGCACTGTTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTTGGCCTg                                                                                                                                                                >  2:86780/3‑251 (MQ=37)
                                                             cccGCGGCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGGGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                                                                                                              <  2:124827/250‑1 (MQ=255)
                                                                    ccAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                                                                                                       >  1:132705/1‑251 (MQ=255)
                                                                         gtacgtTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCa                                                                                                                                                    >  1:1897/1‑249 (MQ=255)
                                                                            cgtTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGAAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCAATGTGGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCgg                                                                                                                                               >  1:103200/1‑250 (MQ=255)
                                                                                gATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGc                                                                                                                                            <  1:140853/250‑1 (MQ=255)
                                                                                gATAATAGGCAATCGAGGCGGCCCAGTTACCGGCAACAGAGAGGTAGGGACTTAACGGCATGAAGGTGCCCTCAGTTAAGAGAAGACAGTTCAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTAGCATCCGGTAGACGGCACTGTTGCCGGATGAGGCGGAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAATTGTAGGCCTGATCAGAcccgacaacgtaga                                                                                                                                            >  2:335168/1‑237 (MQ=21)
                                                                                  taataGGCAATCGCGTCGGAACAGTTACAGGCAAAAGAGAGGCAGGGACTTAACGGCATGATGGTGACATCAGTTAAGAGAAGCCAGTTAAGTGTAGGTCGAAATTTATACAGATGAGAGGCGGACTGATAAGACGCGCGAGCGTCGCATCCGGAAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGATGTGTCGGCCTGATAAGACCCACAAGCACCGCATc                                                                                                                                         >  2:77417/1‑251 (MQ=25)
                                                                                        gCAATCGCATCGAAACAGTTACAGGCCAAAGAGAGGTAGGGACTTAACGGCATGATGGGGACCTCCGTTAAGAGCAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGCGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGAGGGGGCGGGGATTTGTAGGCCTGATAAAACGCGCAAGCGTCGCATCCGGc                                                                                                                                     >  2:184206/1‑249 (MQ=255)
                                                                                         cAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGt                                                                                                                                  <  1:350403/251‑1 (MQ=255)
                                                                                          aaTCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTc                                                                                                                                 <  1:184206/251‑1 (MQ=255)
                                                                                               cgTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAc                                                                                                                            >  1:236842/1‑251 (MQ=255)
                                                                                                           tACCGGCAAAAGAGAGGTAGGGACTTAACGTCATGATTGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTTCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCTGTCGGCACCATTGCCGGAt                                                                                                                 <  1:382563/250‑1 (MQ=255)
                                                                                                                   aaaGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTa                                                                                                         >  1:162802/1‑250 (MQ=255)
                                                                                                                    aaGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAAACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGCATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCCAGCAGTCGGAACCATTGCCGGATgcggcggaa                                                                                                        >  2:350140/1‑247 (MQ=25)
                                                                                                                     agagagGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAAc                                                                                                      <  1:303463/251‑1 (MQ=255)
                                                                                                                                  cTTAACGGCATGATGGAGACCTCAGTTAAGAGAAGCCAGATCAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTCGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGAGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGGAAACGCCTTACCCGGcc                                                                                         >  2:138487/1‑251 (MQ=18)
                                                                                                                                     aaCGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTGCGAAATTTATACAGATGAGAGGCGGCCTGAGAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                      <  2:57431/251‑1 (MQ=255)
                                                                                                                                      aCGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTAAACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGGAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                                                                     <  2:309228/251‑1 (MQ=255)
                                                                                                                                            tgatgGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGt                                                                               >  1:284573/1‑251 (MQ=255)
                                                                                                                                                             aaGAGAAGCCAGCTAAGTGTAGTTCGCAATTTATACAGATGAGAGGCGGCCAGATAAGACGCCCGAGCGTCGCATCAGATAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTCGGCCTGCTAAGAAGCGCAAGCGTCGCACCAGGCAGTCGGCACCATTGCCGGCTGCGGCGCAAAAGCCTAATCCGGCCTACGGAGGGGGCGGGCATTTGTAGGcc                                                              >  2:47552/1‑251 (MQ=17)
                                                                                                                                                              agagAAGACAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGCTGCGGCGTAAACGCCTTATCCGGCATACGGAGGGTGCGGGAATTTGTAGGCCt                                                             >  1:345466/1‑251 (MQ=255)
                                                                                                                                                                       aGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGCTACGACGCGCGAGCGTCGCATCAGGAAGTCGGCACTGTTGCCGGCTGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCCTCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGAGGGCGCGGGCATTTGTAGGCTTGATTAGAcg                                                    >  2:50858/1‑251 (MQ=17)
                                                                                                                                                                       aGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcg                                                    >  1:42837/1‑251 (MQ=255)
                                                                                                                                                                            aGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCCCCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGCTAAGACGCGCa                                                >  2:190882/1‑250 (MQ=255)
                                                                                                                                                                                 aGTTCTAAATTGTTACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCTGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGTCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAATCGTCGCATCAGGCAGTCGGCACCATTGCCGTATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                          <  1:238027/251‑1 (MQ=14)
                                                                                                                                                                                         aTTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg                                  >  1:64887/1‑251 (MQ=255)
                                                                                                                                                                                               aCAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGAGCGGGAATTTGTAGGCCTGATAAGAAGCGCAAGCGTCGCATCAGGCAGTCg                            >  1:23847/1‑251 (MQ=255)
                                                                                                                                                                                                    tGAGAGGAGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGTTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAcc                       >  2:265533/1‑251 (MQ=255)
                                                                                                                                                                                                        agGCGGACTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCTACGGAGGGTGCGGGAAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTg                   <  1:104283/251‑1 (MQ=255)
                                                                                                                                                                                                           cggcCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCg                >  1:58565/1‑251 (MQ=255)
                                                                                                                                                                                                             gcCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCCTCCGGCCGTAGGCACCGTTGCCgg               >  2:280900/1‑250 (MQ=255)
                                                                                                                                                                                                                 gATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcg          <  1:190882/251‑1 (MQ=255)
                                                                                                                                                                                                                  aTAAGCCGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGACTTATCCGGCCTCCGGAGGGTGCGGGACTTTGTAGACCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATACGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGACAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgggg         >  2:197092/1‑248 (MQ=14)
                                                                                                                                                                                                                    aaGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTCAACGCCTTATCCGGCATAAGGAGGGTGCGGGCATTTGTAGGCATGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGCATTTGTAGGACTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGCTGCCGGATgcggcg       >  1:253720/1‑251 (MQ=14)
                                                                                                                                                                                                                    aaGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACCCGCACGCGTCGCATCAGGCAGGCGGCACCGTTGCCGGTTgcggcg       >  2:110458/1‑251 (MQ=255)
                                                                                                                                                                                                                    aaGACGCGCGAGCGTCGCATCAGGGAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggc        <  1:335168/250‑1 (MQ=255)
                                                                                                                                                                                                                         gcgcgAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAAc  <  1:288171/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                       
GTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAAC  >  NC_000913/4325601‑4326068

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