Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913718,6590TA44.1% 99.6 / 27.1 68D601V (GAT→GTT) speFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base T (12/26);  new base A (30/0);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 7.15e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGA  >  NC_000913/718512‑718805
                                                                                                                                                    |                                                                                                                                                   
ccTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCAt                                                                                                                                                  <  2:680674/151‑1 (MQ=255)
       cATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTc                                                                                                                                           <  2:799304/151‑1 (MQ=255)
             gCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGaca                                                                                                                                     <  2:520175/151‑1 (MQ=255)
                gACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGGTACAAAATATGCCttt                                                                                                                                            >  2:802465/1‑141 (MQ=255)
                  cGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGAAAAAAAAAAGCCTTTTCTTACACaa                                                                                                                                  >  1:1033918/1‑149 (MQ=255)
                  cGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGCTACAAAAAACATGCtttttttttcata                                                                                                                                   >  2:1105840/1‑139 (MQ=255)
                       cAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGGTACAAAACAAAGCTTTCCTGAAATAACAGaa                                                                                                                             >  2:287872/1‑148 (MQ=255)
                       cAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTGCACCACTGACg                                                                                                                            >  2:274813/1‑150 (MQ=255)
                               cTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAATCATGTATTTTCTTTCACCACTGACGCAGGGgat                                                                                                                    >  1:295012/1‑150 (MQ=255)
                                 tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGCAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGtatat                                                                                                                  <  1:680416/150‑2 (MQ=255)
                                    ccGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCgg                                                                                                              <  1:737720/151‑1 (MQ=255)
                                     cgcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAACATTCCTTTTCTTACACAACTGACGCAGGGGATaaaccg                                                                                                              >  2:886182/1‑146 (MQ=255)
                                      gcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTGACACAACAGACGCAGGGGATaaacag                                                                                                              >  2:894609/1‑145 (MQ=255)
                                      gcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTCCTGACACAACTGACGCAGGGGATAAACGGCg                                                                                                            >  2:1087669/1‑151 (MQ=255)
                                       cgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAACAATCATGTCTTTTCTTTAACCACTGACGCAGGGGATaaac                                                                                                                >  1:464020/1‑144 (MQ=255)
                                         tAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCATTTTCTTACACAACTGACGCAGGGGATaaac                                                                                                                >  2:112859/1‑142 (MQ=255)
                                           aaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAgcg                                                                                                       <  2:1059602/151‑1 (MQ=255)
                                                gCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCtt                                                                                                  <  2:1078053/151‑1 (MQ=255)
                                                gCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCtt                                                                                                  <  1:96653/151‑1 (MQ=255)
                                                gCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCATTTTCTTACACAACTGACGCACGGGATAACCgg                                                                                                              >  1:475387/1‑139 (MQ=255)
                                                  gTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGCCACCACTTACGCAGGGGATAACCGGCGCAGCGCTCTTc                                                                                                 >  1:392295/1‑150 (MQ=255)
                                                     gTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGc                                                                                             <  2:161164/151‑1 (MQ=255)
                                                            tcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTg                                                                                       <  2:720765/150‑1 (MQ=255)
                                                                gCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGa                                                                                   <  1:967691/150‑1 (MQ=255)
                                                                gCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGAt                                                                                  <  1:615974/151‑1 (MQ=255)
                                                                  ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAAAATGCATTTTCTTACACAACTGACGCAGGGGATAAACCGCGTAGCGCTCTTCATTCTTTTTGt                                                                                      >  2:225774/1‑145 (MQ=255)
                                                                      tCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAg                                                                             <  2:483562/150‑1 (MQ=255)
                                                                        aTGCTGATGCGTGGGAAGTGCTCCTTACGGACCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAgg                                                                            <  1:247908/149‑1 (MQ=255)
                                                                          gCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAg                                                                         <  2:537991/150‑1 (MQ=255)
                                                                            tGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATGCCTTTTCTTTCACAACTGACGCAGGGGATAAACCGCGTAGCGCTCTTCATTCTGTTTTTAGATGggagggagga                                                                       >  2:806527/1‑148 (MQ=255)
                                                                             gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCa                                                                       <  1:205421/149‑1 (MQ=255)
                                                                             gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCa                                                                       <  2:763472/149‑1 (MQ=255)
                                                                             gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAATCATGCATTTCCTGACACAACCGCCGCCGGGGATAACCCGCGGAGCGCCCTTCATTCTGTTTGTTGGTGGGAggga                                                                          >  2:799475/1‑144 (MQ=255)
                                                                               cgcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATGCATTTCCTTACCCAACAGACGCAGGGGATAAACCGCGTAGCGCTCTTTCTGCTGTTTGTAGGTGGAAGGGAGCAcct                                                                    >  1:1133195/1‑148 (MQ=255)
                                                                                  gTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATTTATTTTCTTTCCCAACCGTCGCAGGGGGTAAACGGCGGAGCGCCCTTTCTTCTGTTTTTAGGTGGAAGGGAGCAATTCt                                                                  >  1:613033/1‑149 (MQ=255)
                                                                                   tGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAATATGCAATTTCTTACACAACTAACGCAGGGGGTAAACCGCGGAGGGCTCTTCACGATGTTTGTTGAGGGaaagaaacaattc                                                                   >  2:698538/1‑136 (MQ=255)
                                                                                    gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCATTTTCTTACACAAATGACGCAGGGGATAAACGGCGGAGCGCTCTTCATGCTGTTTGTAGATGGAAGGGAGCACTTCt                                                                  >  2:578082/1‑147 (MQ=255)
                                                                                     ggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTAAAAAAAATGCATTTTCTTTCACAACTTACGCAGGGGGTAAACCGCGGGGCGCCCTTCTTTTCGTTTGTTGATGGGAGGGAGCACtttt                                                                  >  1:82327/1‑144 (MQ=255)
                                                                                        aGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAACAAACAGGTCTTTTCTTACACAAACAACGCAGGGGGTAAAAGGGGGAGGGCCCTTCATGTTGTTTTTGGGTGggagggagaaactttcgtaa                                                             >  1:536031/1‑136 (MQ=255)
                                                                                           gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTAACAAACATTCTTTTTCTTTCACCACTGTCGCCGGGGAGAACCGGCGGAGCGCCTCTCCTTCTCTTTTTTTATGGGAGGGAGCACCTCTCCTCAt                                                           >  1:718479/1‑146 (MQ=255)
                                                                                           gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCATTTTCTTACACAAAAGACGCAGGGGAAAAACCGCGGAGAGCGCTTCATGCTGTTTTTAGAGGGGAGGGAGCACTTCTTCTAAttggg                                                       >  2:188607/1‑146 (MQ=255)
                                                                                            cTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAAAATGCATTTTCTGAAACAAATGAAGCAGGGGATAAACGGCGGAGCGCCCTTTATTCTGTTTGTTGAGGGGAGGGAGGACTTCTTCTAATGGcgcg                                                     >  2:455986/1‑151 (MQ=255)
                                                                                             tCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAATATGCATTTTCTTGTCCAAAATAAGCAGGGGGTAAACGGCGTAGCGCTCTTCATTCTTTTTTTTGATGGGAGGGAGCAATTCTGCTTATgg                                                         >  1:1148871/1‑146 (MQ=255)
                                                                                              ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTc                                                   <  1:169862/151‑1 (MQ=255)
                                                                                                 tACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTAAAAAATATTTATTTTCTTTCACAACTTACCCAGGGGAGAAACGGCGTAGGGCTCTTTCTTTTGTTTTTTGATGGGAGGTAGAACCTTTTCTTATTGCGcggcggg                                                >  1:734277/1‑150 (MQ=255)
                                                                                                   cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTAAAAAAAATTTATTTTCTTTCACAACAAACGCGGGGGATAAACGGCGGGGCGCTCTTCTTTCTGTTTTTTGGTGGGAAGGAGGACTTTTGCTTAATgggcggcgg                                                  >  1:68193/1‑141 (MQ=255)
                                                                                                      aaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGAGCAATAAAACATTCATTTTTTTTCACCACTGACGCCGGGGGAAAACGGCGGAGCGCCCTTCATTGTTTTTTTTGGAGGAAGGGAGCACATCTGCTAAtaggtagccggga                                               >  1:971571/1‑135 (MQ=255)
                                                                                                        cATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATTTATTTTCTTATACAACTTACGCAGGGGATAAACGGCGGAGCGCTCTTTTTGTTGTTTTTTGGTGGGAGGGAGCACTTTTGCTAATGGGGcggcggggt                                              >  1:104460/1‑143 (MQ=255)
                                                                                                          tctcTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAATATTTTTTTTCTTGCACAACATACGCAGGGGAGAAACGGCGGAGGGCCCTTTATGTTGTTTTTTGATGGAAGGGAGCACTTCTTCTTATGGGGcggcggg                                                >  1:407950/1‑141 (MQ=255)
                                                                                                            tctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCATTTTCTTACACAAATGACGCAGGGGATAAACGGCGGAGCGCTCTTCCTGCTGTTTGTTGATGGAAAGGAGCACTTCTGCTAATGGCGCGTCGGGCTCAAGCAggt                                     >  2:1038067/1‑149 (MQ=255)
                                                                                                             ctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCTTTTTCTTTCACCACCGACGCAGGGGGTAACCGGCGGAGCGGCCTTCTTGTTTTTTGTTGATGGAAGGGAGGAACTCTGCCCATggcggggcgggcgccagacggtt                                    >  1:416097/1‑138 (MQ=255)
                                                                                                               ttttGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGttt                                    <  2:531925/149‑1 (MQ=255)
                                                                                                                tttGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGttttt                                  <  2:536205/150‑1 (MQ=255)
                                                                                                                     aGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGc                            <  2:1053686/151‑1 (MQ=255)
                                                                                                                      gttgttTCACGTTGTGGCGGGGCTACAAAAAATGCCTTTTCTTTCACAACCGACGCCGGGGGTAAACGGCGGAGCGCCCTTCCTGCTGTTTTTTGATGGAAGGGAGCACTTCTTCTTATGGCGCGGCGGGCTCAAACAGgttttttaa                              >  1:972639/1‑144 (MQ=255)
                                                                                                                            tCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACca                      <  1:353885/150‑1 (MQ=255)
                                                                                                                                   gtgGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGca              <  1:772123/151‑1 (MQ=255)
                                                                                                                                       cGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca           >  1:839211/1‑150 (MQ=255)
                                                                                                                                         ggCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGt        <  2:1001592/151‑1 (MQ=255)
                                                                                                                                           gCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCACCAAGtt       <  1:980387/151‑1 (MQ=255)
                                                                                                                                           gCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGtt       >  2:205961/1‑151 (MQ=255)
                                                                                                                                             aTACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTg      >  2:900578/1‑150 (MQ=255)
                                                                                                                                             aTACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGCCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTg      >  2:205279/1‑150 (MQ=255)
                                                                                                                                               aCAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg   >  2:407966/1‑151 (MQ=255)
                                                                                                                                                cAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCt    <  1:499710/149‑1 (MQ=255)
                                                                                                                                                cAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg   >  2:374115/1‑150 (MQ=255)
                                                                                                                                                 aaaTCATGTATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGa  >  1:76218/1‑150 (MQ=255)
                                                                                                                                                  aaTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGa  >  2:77200/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGA  >  NC_000913/718512‑718805

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