Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,487 T→G A101A (GCA→GCC)  ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,4870TG98.2% 163.9 / ‑3.8 56A101A (GCA→GCC) ytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (29/25);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

GCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGT‑TCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCC  >  NC_000913/4624243‑4624718
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
gcGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa                                                                                                                                                                                                                                    <  1:139829/251‑1 (MQ=255)
 cGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa                                                                                                                                                                                                                                    >  1:42438/1‑250 (MQ=255)
 cGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa                                                                                                                                                                                                                                    >  1:348039/1‑250 (MQ=255)
       gtgcatcactacggcGATTGCCATTGAAAGCAGCATCATCATGCGTAAAATGTTAGTGGTGTTATCCGCCTGTTGTGCTACTGTGGCGGTCGTATTGGTGTCGAGCGGCAGGCGCAGATAGCCGTGCGGTCCTTTTTTCCCGGCAATTGTTTCGACAATCTGCTGGTTAAAATAGCCGCCCGTTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCAGCTCGACGCTTTCGCCAGAACGGGCGATAAGATCGc                                                                                                                                                                                                                             <  2:342863/237‑1 (MQ=255)
           aGCATTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCTTAAAATGTTAGTGGTGTTATCCCCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCTAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCTCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTg                                                                                                                                                                                                                         <  1:272932/251‑1 (MQ=255)
               ctacgccggttgctcgTGAAAGCAGCAACAGCAGGCGTAAAATGTTAGTGGGGGTAGCATCCTGTTGGGCTTCGTTGGCTGGCGTATGGGTGTTGCGCGTCAGGCGCGGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAACCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGt‑t                                                                                                                                                                                                                      <  2:375496/234‑1 (MQ=255)
                   gccGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCAt                                                                                                                                                                                                                >  2:567033/2‑251 (MQ=255)
                                            aTCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCt                                                                                                                                                                                       <  1:390533/251‑1 (MQ=255)
                                                          gTTAGTGGAGTTATCCACCTGTTGGGCTTCGGTGGCGCGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGCGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTCAAATAGACGCCCGCTTTTTTACAGTCGCGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCCTACACGCCCGCGTCGTGGATACGGCTTTCCTCCggttact                                                                                                                                                                         >  2:340132/1‑245 (MQ=255)
                                                             aGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCAAGGATACGGCTTTCATCCGTTAACTGAt                                                                                                                                                                      >  2:512957/1‑251 (MQ=255)
                                                                  tgtTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGCCGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATGCGGCTTTCATCCGTTAACTGCTCGAgg                                                                                                                                                                 >  2:463936/1‑251 (MQ=255)
                                                                     tATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGGGGATACGGCTTTCATCCGTTacttgatgtaggatc                                                                                                                                                              >  2:328791/1‑237 (MQ=255)
                                                                     tATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATc                                                                                                                                                              >  1:98986/1‑251 (MQ=255)
                                                                       tCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGc                                                                                                                                                            >  1:31257/1‑251 (MQ=255)
                                                                              gTTGGTCTTCGGGGGCGAGCGTATGGGTGTCGGGGGTCAGGCGCAGATAGCCGAGGGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGAGCGCCTGaa                                                                                                                                                       <  2:42438/249‑1 (MQ=255)
                                                                                 gggcttgggggtgtgacGTAGGGGGGTCGAGGGTCAGGCCCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCGTCTGGTAAAAATAGCCGCCCGCTTTTTTCCCGGCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGGCGAGGAGACGGGTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGt                                                                                                                                                  <  2:98428/235‑1 (MQ=255)
                                                                                 gggCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGt                                                                                                                                                  >  1:98428/1‑251 (MQ=255)
                                                                                   gCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGttt                                                                                                                                                <  1:512957/251‑1 (MQ=255)
                                                                                      tCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGttttt                                                                                                                                              <  1:201613/250‑1 (MQ=255)
                                                                                       cGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTc                                                                                                                                             <  1:218311/250‑1 (MQ=255)
                                                                                               aGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCAt                                                                                                                                                                                                                <  1:242688/175‑1 (MQ=255)
                                                                                               aGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCAt                                                                                                                                                                                                                >  2:242688/1‑175 (MQ=255)
                                                                                                      gggTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGT‑TCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGtc                                                                                                                             >  1:226254/1‑250 (MQ=255)
                                                                                                       ggTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTACGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTAc                                                                                                                            <  1:328791/251‑1 (MQ=255)
                                                                                                         tgtCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGc                                                                                                                          >  1:148078/1‑251 (MQ=255)
                                                                                                                 gtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCgg                                                                                                                  >  1:58007/1‑251 (MQ=255)
                                                                                                                    aGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGc                                                                                                               >  1:1217/1‑251 (MQ=255)
                                                                                                                                    actatcagcgttcagtatcAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGaa                                                                                                                                                       <  1:487136/177‑1 (MQ=255)
                                                                                                                                    actatcagcgttcagtatcAATTGGTTAGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGaa                                                                                                                                                       >  2:487136/19‑195 (MQ=255)
                                                                                                                                      ggTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGCGGCTCGCCTGACTGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGGGCAAGGTTCAGCGtttcctg                                                                                             >  2:329231/1‑245 (MQ=255)
                                                                                                                                            tttttCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTGCCGGCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCAGCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTGCCTGAcgcg                                                                                        <  2:31257/250‑1 (MQ=255)
                                                                                                                                                 cccggCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAggg                                                                                  >  1:384792/5‑251 (MQ=255)
                                                                                                                                                 cccggCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAgg                                                                                   <  1:567033/246‑1 (MQ=255)
                                                                                                                                                  ccggCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCCgggg                                                                                 >  2:308691/4‑250 (MQ=255)
                                                                                                                                                    ggCAATTGTGTCGACATTCTGCTGGTTAAAATAGCCGCCAGGTTTTTTACCGCGGAGCGCCAGCCGGTCGCGCACGTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑GCGTCATACACGCCCGCGTCGAGGATACGGCTTTCAGCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGc                                                                               <  2:98986/250‑1 (MQ=255)
                                                                                                                                                       aaTTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCggg                                                                            >  1:425255/1‑251 (MQ=255)
                                                                                                                                                       aaTTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCggg                                                                            <  1:330109/251‑1 (MQ=255)
                                                                                                                                                               cGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACCTTCATGGTTACCTGACGCGCCCGGGTGCGGGCCAGttct                                                                    >  2:419664/1‑251 (MQ=255)
                                                                                                                                                                               aaaaTAGCCGCCCGCTTTTTAACCGCGGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTAc                                                    <  2:150537/251‑1 (MQ=255)
                                                                                                                                                                                     cccgcccgcTTTTTTACCGGCGAGCGCCAGCCGGTCGCGCACTGCGACGTCTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCAGACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTc                                              <  2:410702/250‑1 (MQ=255)
                                                                                                                                                                                         cccgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGCTACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTCCGTTGGCGCtg                                          >  2:664578/4‑251 (MQ=255)
                                                                                                                                                                                         cccgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGAGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCCGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGATGTCGCtg                                          >  2:621571/4‑251 (MQ=255)
                                                                                                                                                                                           cgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCTCTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGGCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGGTACCTGACGCGCCAGGGTGCGGGTCAGTTCTTCCAGCTGTGGATTACGTTGTCGCtggt                                        <  2:348039/250‑1 (MQ=255)
                                                                                                                                                                                            gtttgtttACCGGCGAGTGCCCGCTGATCGTGCACTTCGACGCTTTCGCCAGAACCGGCGATAAGATCGCCT‑GTGTCGTCAGACCCGCCCGCGTCGAGGATACGGCTTTCATCCTGTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGTCGAGTCGGTACTCATCAGCGGTGCAACGTGCATCTTTACCTGTCGCGCCAGGGTTCGTGCCAGTTCTTCCAGCATTGGATTActttgtctctggga                                       <  2:114783/246‑14 (MQ=255)
                                                                                                                                                                                              ttttttCCCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCATAGCGGGCGTTTAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGTCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGTtggtgg                                      <  2:432368/250‑1 (MQ=255)
                                                                                                                                                                                                ttttACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGtt                                    <  1:463936/250‑1 (MQ=255)
                                                                                                                                                                                                         cGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGGGGTTTTGAATaatc                          >  2:305276/1‑249 (MQ=255)
                                                                                                                                                                                                             cgcCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATg                      >  1:522389/1‑251 (MQ=255)
                                                                                                                                                                                                                cAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGCATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCCACGTTCAGCGTTACCTCACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACg                   >  2:225354/1‑251 (MQ=255)
                                                                                                                                                                                                                   ccGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTc                >  1:491351/1‑251 (MQ=255)
                                                                                                                                                                                                                   ccGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCAACAGCGGAGCAACGTTCAGCGTTACCTGTCGCGCCAGGGAGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTTACTAAACCCTGACGCTc                >  1:454807/1‑251 (MQ=255)
                                                                                                                                                                                                                      gTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCt             <  2:491351/251‑1 (MQ=255)
                                                                                                                                                                                                                               ttCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGGTTACGTTGTCGCTGGTGGTTTTGACTCAACCATGACGCTCCCTGCATCAgcg    >  2:317253/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAgc     >  1:194562/1‑247 (MQ=255)
                                                                                                                                                                                                                                  gACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGcc  >  1:62972/1‑250 (MQ=255)
                                                                                                                                                                                                                                  gACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGT‑TCGTCATACACGCCCGCGTCGAGGATACGGCTGTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGGACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAgc     <  2:194562/247‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
GCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGT‑TCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCC  >  NC_000913/4624243‑4624718

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