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,4870TG94.1% 140.6 / ‑0.7 51A101A (GCA→GCCytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (22/26);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

CAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGA  >  NC_000913/4624241‑4624733
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
cagcGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCg                                                                                                                                                                                                                                                     >  2:140155/1‑250 (MQ=255)
   cGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATaa                                                                                                                                                                                                                                                 <  2:391370/251‑1 (MQ=255)
            cAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCtt                                                                                                                                                                                                                                        >  1:237649/1‑251 (MQ=255)
                   acgccGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGGCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTc                                                                                                                                                                                                                                 >  2:439504/4‑251 (MQ=255)
                                  gaaagcagcAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTAGGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTGTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTc                                                                                                                                                                                                                  <  1:47798/248‑1 (MQ=255)
                                   aaggcagcaACATCAGGCGAAAAATGTTAGTGGTGTTATCCACCTGTTGGACTTCGGAGGCGAGCGTATTGGAGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGATTTTCCCGGCAATTAGTTCGACAATCTGCTGGTTAAAATATCCACCCGCGTTTTTACCGTCGATCGCCACCCGGTCGCGCACTTCGCGACTTTCGCCCGAGCGGGCGATAAGATCACCTTGtacgaaataccagcccgactcg                                                                                                                                                                                                                 >  2:72363/4‑230 (MQ=255)
                                         gcaACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGCGGAt                                                                                                                                                                                                            >  1:358684/1‑250 (MQ=255)
                                             cATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTAGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGc                                                                                                                                                                                                       >  1:35358/1‑251 (MQ=255)
                                                cAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCttt                                                                                                                                                                                                    <  1:391751/251‑1 (MQ=255)
                                                     gTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATcc                                                                                                                                                                                               >  1:341866/1‑251 (MQ=255)
                                                            gTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGTCAGAACGGGCGATAAGATCGCTTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAAc                                                                                                                                                                                         <  1:406633/250‑1 (MQ=255)
                                                                ggggtgTTATCCGCCTGTTGGGCTTCGGTGGCGGGCGTTTGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGTTCCGTTTTCCCCGGCAATTGGTTCGACAACCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATc                                                                                                                                                                                    <  2:261567/249‑1 (MQ=255)
                                                                       tATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGAt                                                                                                                                                                              <  1:480212/250‑1 (MQ=255)
                                                                       tATCCACCGGTTGGGCTCCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGAt                                                                                                                                                                              <  2:165367/250‑1 (MQ=255)
                                                                                gTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAAt                                                                                                                                                                     <  1:99725/249‑1 (MQ=255)
                                                                                    ggCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGtt                                                                                                                                                                >  1:47189/1‑251 (MQ=255)
                                                                                      cTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGtttt                                                                                                                                                              >  1:481573/1‑251 (MQ=255)
                                                                                            tGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATcc                                                                                                                                                        >  2:129227/1‑251 (MQ=255)
                                                                                               cGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGCTCCGGCTTTCATCCGTTCACTGATCGAGGATCGGCTGGATGCgatgttaatcccgc                                                                                                                                                     >  2:365576/1‑238 (MQ=255)
                                                                                               cGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGcc                                                                                                                                                                                                                                          >  2:310494/1‑166 (MQ=255)
                                                                                               cGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGcc                                                                                                                                                                                                                                          <  1:310494/166‑1 (MQ=255)
                                                                                                    gTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTc                                                                                                                                                <  1:241790/251‑1 (MQ=255)
                                                                                                     tATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCCGCAGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGCTCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTCATCAAGGATCGCCGGCATGCGTTTTTCATCCggcggagcg                                                                                                                                               >  2:68209/1‑244 (MQ=255)
                                                                                                               gagcgtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCTCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCa                                                                                                                                     <  1:120471/247‑1 (MQ=255)
                                                                                                                 gcgtCAGGCGCAGATGGCCGAGCGGTCCTGTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTGAAAATAGCCGCCCTCTTTTTTACCGGCGAGCGCCAGCCGGTCGCGCTCTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGc                                                                                                                                   <  2:381792/249‑1 (MQ=255)
                                                                                                                     cAGGCGCGGATGGCCCAGGGGTCCGGTTTTCCCGCCAAGTGGGCCGACAATCTGCTGGTTAAAATAGCCGCCCGTTTTTTTACCGTGGAGCGCCAGCCGGTCGCGCACTTCGACGCGTTCGCCAGAACGGGCTATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTGTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTg                                                                                                                               <  2:287073/251‑1 (MQ=255)
                                                                                                                       ggCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCg                                                                                                                                                                                                                                   <  2:362891/149‑1 (MQ=255)
                                                                                                                       ggCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCg                                                                                                                                                                                                                                   >  1:362891/1‑149 (MQ=255)
                                                                                                                       agctcAGATAGCCGAGCGGGCCGTTTTTCCCGGAAATTGGTTCGACAATCGGCTGGTTAAAATAGCCGCCCGCTTTTGTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGGCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCa                                                                                                                             <  2:99426/247‑1 (MQ=255)
                                                                                                                           cAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCa                                                                                                                             >  1:99426/1‑247 (MQ=255)
                                                                                                                                aGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGAGCAACGTGCAGc                                                                                                                    >  2:108130/1‑251 (MQ=255)
                                                                                                                                       cGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGACCGCTTAAGCGGTGCAACGTTCCGCGTTACCt                                                                                                             >  2:418253/1‑251 (MQ=255)
                                                                                                                                             gTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGAcgcg                                                                                                       <  1:304369/251‑1 (MQ=255)
                                                                                                                                                                     aaaTCTGTTGGGTAAAATAGCCGCCCGCTTTTTTACCGTCGGGCGCCAGCCGGTCGCGCACTTTGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTcc                                                                                <  2:132392/250‑1 (MQ=255)
                                                                                                                                                                           gctgGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAATGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAg                                                                              <  2:47189/247‑1 (MQ=255)
                                                                                                                                                                           gctgGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAg                                                                              <  1:409850/247‑1 (MQ=255)
                                                                                                                                                                                    aaTAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACgt                                                                 <  1:129227/251‑2 (MQ=255)
                                                                                                                                                                                         ccgcacgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATAGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTAGGTACGCATCAGCGGTGCCACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCAGGGGATTACGTTGTCg                                                            >  2:430795/7‑251 (MQ=255)
                                                                                                                                                                                           gcccgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTTCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCt                                                          <  1:301412/247‑1 (MQ=255)
                                                                                                                                                                                                 ttttttACCGTCGAGCGGCAGCCGGTCGCGCACTTCGCCGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGTATGCGATTTTCATCCGGCGAGTCGGGACGCCTCAGCGGTGAACCGTTCAGCGTTACCTGACGAGCCTAGGTGCGAGCAAGTTCTTCCAGCTGTGGTTTTCGTTGTCGCtggtggg                                                    >  2:263503/1‑250 (MQ=255)
                                                                                                                                                                                                      tACCGTCGGGCGCCAGCCGGGCGCGCACTTCGTCGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGAGCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGttttt                                                <  2:424130/250‑2 (MQ=255)
                                                                                                                                                                                                                 gcCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGGGCGATAAGATAGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCCGTTCTTCCCGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATg                                     >  1:375644/1‑250 (MQ=255)
                                                                                                                                                                                                                   cAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATAAGGCTTTCATACGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTAGATGGAGGTTTTGACTAAACCATGa                                    >  1:320636/1‑249 (MQ=255)
                                                                                                                                                                                                                       cGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTcc                              <  1:439504/251‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGGTTTGACTAAACCATGACGCTcc                              >  2:228040/1‑251 (MQ=255)
                                                                                                                                                                                                                             cgcACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCAt                        <  1:240000/251‑1 (MQ=255)
                                                                                                                                                                                                                                aCTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAg                     <  1:276184/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gcagCTTTCGCCAGAACGGGCGATAAGATCGCCGTGTTCGTCAGTCACGCCCGCGTCGAGGATACGGTTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCa                <  2:198195/248‑1 (MQ=255)
                                                                                                                                                                                                                                        gCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCAAGGGTGCGGGCCAGCTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGCCGCTACCTGCATCAGCGCCACt              >  1:454325/1‑250 (MQ=255)
                                                                                                                                                                                                                                          cTTTCGCCAGAACGGGCGCTAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTACACCATGACGCTCCCTGCATCAGCGCCACTac            >  2:314930/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    aaCGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGa  <  1:241672/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
CAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGA  >  NC_000913/4624241‑4624733

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