Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,517 A→G G91G (GGT→GGC 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,5170AG96.2% 147.8 / ‑4.3 52G91G (GGT→GGCytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (29/21);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

ATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCA  >  NC_000913/4624287‑4624761
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aTCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCt                                                                                                                                                                                                                                  <  1:390533/251‑1 (MQ=255)
              gTTAGTGGAGTTATCCACCTGTTGGGCTTCGGTGGCGCGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGCGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTCAAATAGACGCCCGCTTTTTTACAGTCGCGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCCTACACGCCCGCGTCGTGGATACGGCTTTCCTCCggttact                                                                                                                                                                                                                    >  2:340132/1‑245 (MQ=255)
                 aGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCAAGGATACGGCTTTCATCCGTTAACTGAt                                                                                                                                                                                                                 >  2:512957/1‑251 (MQ=255)
                      tgtTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGCCGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATGCGGCTTTCATCCGTTAACTGCTCGAgg                                                                                                                                                                                                            >  2:463936/1‑251 (MQ=255)
                         tATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGGGGATACGGCTTTCATCCGTTacttgatgtaggatc                                                                                                                                                                                                         >  2:328791/1‑237 (MQ=255)
                         tATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATc                                                                                                                                                                                                         >  1:98986/1‑251 (MQ=255)
                           tCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGc                                                                                                                                                                                                       >  1:31257/1‑251 (MQ=255)
                                  gTTGGTCTTCGGGGGCGAGCGTATGGGTGTCGGGGGTCAGGCGCAGATAGCCGAGGGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGAGCGCCTGaa                                                                                                                                                                                                  <  2:42438/249‑1 (MQ=255)
                                     gggcttgggggtgtgacGTAGGGGGGTCGAGGGTCAGGCCCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCGTCTGGTAAAAATAGCCGCCCGCTTTTTTCCCGGCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGGCGAGGAGACGGGTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGt                                                                                                                                                                                             <  2:98428/235‑1 (MQ=255)
                                     gggCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGt                                                                                                                                                                                             >  1:98428/1‑251 (MQ=255)
                                       gCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGttt                                                                                                                                                                                           <  1:512957/251‑1 (MQ=255)
                                          tCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGttttt                                                                                                                                                                                         <  1:201613/250‑1 (MQ=255)
                                           cGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTc                                                                                                                                                                                        <  1:218311/250‑1 (MQ=255)
                                                          gggTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGtc                                                                                                                                                                        >  1:226254/1‑250 (MQ=255)
                                                           ggTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTACGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTAc                                                                                                                                                                       <  1:328791/251‑1 (MQ=255)
                                                             tgtCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGc                                                                                                                                                                     >  1:148078/1‑251 (MQ=255)
                                                                     gtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCgg                                                                                                                                                             >  1:58007/1‑251 (MQ=255)
                                                                        aGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGc                                                                                                                                                          >  1:1217/1‑251 (MQ=255)
                                                                                        actatcagcgttcagtatcAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGaa                                                                                                                                                                                                  <  1:487136/177‑1 (MQ=255)
                                                                                        actatcagcgttcagtatcAATTGGTTAGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGaa                                                                                                                                                                                                  >  2:487136/19‑195 (MQ=255)
                                                                                          ggTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGCGGCTCGCCTGACTGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGGGCAAGGTTCAGCGtttcctg                                                                                                                                        >  2:329231/1‑245 (MQ=255)
                                                                                                tttttCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTGCCGGCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCAGCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTGCCTGAcgcg                                                                                                                                   <  2:31257/250‑1 (MQ=255)
                                                                                                     cccggCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAggg                                                                                                                             >  1:384792/5‑251 (MQ=255)
                                                                                                     cccggCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAgg                                                                                                                              <  1:567033/246‑1 (MQ=255)
                                                                                                      ccggCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCCgggg                                                                                                                            >  2:308691/4‑250 (MQ=255)
                                                                                                        ggCAATTGTGTCGACATTCTGCTGGTTAAAATAGCCGCCAGGTTTTTTACCGCGGAGCGCCAGCCGGTCGCGCACGTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTGCGTCATACACGCCCGCGTCGAGGATACGGCTTTCAGCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGc                                                                                                                          <  2:98986/250‑1 (MQ=255)
                                                                                                           aaTTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCggg                                                                                                                       <  1:330109/251‑1 (MQ=255)
                                                                                                           aaTTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCggg                                                                                                                       >  1:425255/1‑251 (MQ=255)
                                                                                                                   cGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACCTTCATGGTTACCTGACGCGCCCGGGTGCGGGCCAGttct                                                                                                               >  2:419664/1‑251 (MQ=255)
                                                                                                                                   aaaaTAGCCGCCCGCTTTTTAACCGCGGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTAc                                                                                               <  2:150537/251‑1 (MQ=255)
                                                                                                                                         cccgcccgcTTTTTTACCGGCGAGCGCCAGCCGGTCGCGCACTGCGACGTCTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCAGACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTc                                                                                         <  2:410702/250‑1 (MQ=255)
                                                                                                                                             cccgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGCTACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTCCGTTGGCGCtg                                                                                     >  2:664578/4‑251 (MQ=255)
                                                                                                                                             cccgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGAGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCCGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGATGTCGCtg                                                                                     >  2:621571/4‑251 (MQ=255)
                                                                                                                                               cgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCTCTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGGCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGGTACCTGACGCGCCAGGGTGCGGGTCAGTTCTTCCAGCTGTGGATTACGTTGTCGCtggt                                                                                   <  2:348039/250‑1 (MQ=255)
                                                                                                                                                gtttgtttACCGGCGAGTGCCCGCTGATCGTGCACTTCGACGCTTTCGCCAGAACCGGCGATAAGATCGCCTGTGTCGTCAGACCCGCCCGCGTCGAGGATACGGCTTTCATCCTGTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGTCGAGTCGGTACTCATCAGCGGTGCAACGTGCATCTTTACCTGTCGCGCCAGGGTTCGTGCCAGTTCTTCCAGCATTGGATTActttgtctctggga                                                                                  <  2:114783/246‑14 (MQ=255)
                                                                                                                                                  ttttttCCCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCATAGCGGGCGTTTAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGTCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGTtggtgg                                                                                 <  2:432368/250‑1 (MQ=255)
                                                                                                                                                    ttttACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGtt                                                                               <  1:463936/250‑1 (MQ=255)
                                                                                                                                                             cGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGGGGTTTTGAATaatc                                                                     >  2:305276/1‑249 (MQ=255)
                                                                                                                                                                 cgcCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATg                                                                 >  1:522389/1‑251 (MQ=255)
                                                                                                                                                                    cAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGCATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCCACGTTCAGCGTTACCTCACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACg                                                              >  2:225354/1‑251 (MQ=255)
                                                                                                                                                                       ccGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTc                                                           >  1:491351/1‑251 (MQ=255)
                                                                                                                                                                       ccGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCAACAGCGGAGCAACGTTCAGCGTTACCTGTCGCGCCAGGGAGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTTACTAAACCCTGACGCTc                                                           >  1:454807/1‑251 (MQ=255)
                                                                                                                                                                          gTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCt                                                        <  2:491351/251‑1 (MQ=255)
                                                                                                                                                                                   ttCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGGTTACGTTGTCGCTGGTGGTTTTGACTCAACCATGACGCTCCCTGCATCAgcg                                               >  2:317253/1‑251 (MQ=255)
                                                                                                                                                                                      gACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAgc                                                >  1:194562/1‑247 (MQ=255)
                                                                                                                                                                                      gACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGcc                                             >  1:62972/1‑250 (MQ=255)
                                                                                                                                                                                      gACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTGTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGGACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAgc                                                <  2:194562/247‑1 (MQ=255)
                                                                                                                                                                                                             aaaGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAgg                     <  2:58007/250‑1 (MQ=255)
                                                                                                                                                                                                              aaGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAgg                     >  1:155489/1‑250 (MQ=255)
                                                                                                                                                                                                                  tCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAAt                >  1:483642/1‑251 (MQ=255)
                                                                                                                                                                                                                        tGTTAGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGCTCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGGGCGGGCCAGTTCTTCCAGCTGTGGATTAAGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGAATCAGCGCCACTAACAAGGCAAGACAGAACAGGCCAATCACTGc          >  2:171024/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTACACACCCGCGACGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCAGAATACGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCCGGGAGCGGGCCAGTTCTTCCAGCTGGGGATTACGTTGTCGCTTGTGGTTTTGACTAAACCATGACGCTCCCTGCATTAGCGCCACTAACAAGGCACGACAGAACAGGACAATCACTGCCAGATGcc  >  2:50341/1‑250 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
ATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCA  >  NC_000913/4624287‑4624761

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: