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

CCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGA  >  NC_000913/718512‑718805
                                                                                                                                                    |                                                                                                                                                   
ccTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCAt                                                                                                                                                  <  2:680674/151‑1 (MQ=255)
       cATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTc                                                                                                                                           <  2:799304/151‑1 (MQ=255)
             gCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGaca                                                                                                                                     <  2:520175/151‑1 (MQ=255)
                gACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GGGTACAAAATATGCCttt                                                                                                                                            >  2:802465/1‑141 (MQ=255)
                  cGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCGAAAAAAAAAAGCCTTTTCTTACACaa                                                                                                                                  >  1:1033918/1‑149 (MQ=255)
                  cGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑CTACAAAAAACATGCtttttttttcata                                                                                                                                   >  2:1105840/1‑139 (MQ=255)
                       cAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GGGTACAAAACAAAGCTTTCCTGAAATAACAGaa                                                                                                                             >  2:287872/1‑148 (MQ=255)
                       cAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAAACATGCATTTTCTTGCACCACTGACg                                                                                                                            >  2:274813/1‑150 (MQ=255)
                               cTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GGATAAAAATCATGTATTTTCTTTCACCACTGACGCAGGGgat                                                                                                                    >  1:295012/1‑150 (MQ=255)
                                 tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGCAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGtatat                                                                                                                  <  1:680416/150‑2 (MQ=255)
                                    ccGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCgg                                                                                                              <  1:737720/151‑1 (MQ=255)
                                     cgcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATAAAAAACATTCCTTTTCTTACACAACTGACGCAGGGGATaaaccg                                                                                                              >  2:886182/1‑146 (MQ=255)
                                      gcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAAACATGCATTTTCTGACACAACAGACGCAGGGGATaaacag                                                                                                              >  2:894609/1‑145 (MQ=255)
                                      gcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAAACATGCATTTCCTGACACAACTGACGCAGGGGATAAACGGCg                                                                                                            >  2:1087669/1‑151 (MQ=255)
                                       cgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATAACAATCATGTCTTTTCTTTAACCACTGACGCAGGGGATaaac                                                                                                                >  1:464020/1‑144 (MQ=255)
                                         tAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAAAAATGCATTTTCTTACACAACTGACGCAGGGGATaaac                                                                                                                >  2:112859/1‑142 (MQ=255)
                                           aaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAgcg                                                                                                       <  2:1059602/151‑1 (MQ=255)
                                                gCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCtt                                                                                                  <  2:1078053/151‑1 (MQ=255)
                                                gCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCtt                                                                                                  <  1:96653/151‑1 (MQ=255)
                                                gCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAAACATTCATTTTCTTACACAACTGACGCACGGGATAACCgg                                                                                                              >  1:475387/1‑139 (MQ=255)
                                                  gTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGCCACCACTTACGCAGGGGATAACCGGCGCAGCGCTCTTc                                                                                                 >  1:392295/1‑150 (MQ=255)
                                                     gTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGc                                                                                             <  2:161164/151‑1 (MQ=255)
                                                            tcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTg                                                                                       <  2:720765/150‑1 (MQ=255)
                                                                gCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGa                                                                                   <  1:967691/150‑1 (MQ=255)
                                                                gCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGAt                                                                                  <  1:615974/151‑1 (MQ=255)
                                                                  ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GGATACAAAAAATGCATTTTCTTACACAACTGACGCAGGGGATAAACCGCGTAGCGCTCTTCATTCTTTTTGt                                                                                      >  2:225774/1‑145 (MQ=255)
                                                                      tCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAg                                                                             <  2:483562/150‑1 (MQ=255)
                                                                        aTGCTGATGCGTGGGAAGTGCTCCTTACGGACCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAgg                                                                            <  1:247908/149‑1 (MQ=255)
                                                                          gCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAg                                                                         <  2:537991/150‑1 (MQ=255)
                                                                            tGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCCTACAAAACATGCCTTTTCTTTCACAACTGACGCAGGGGATAAACCGCGTAGCGCTCTTCATTCTGTTTTTAGATGggagggagga                                                                       >  2:806527/1‑148 (MQ=255)
                                                                             gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCa                                                                       <  1:205421/149‑1 (MQ=255)
                                                                             gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCa                                                                       <  2:763472/149‑1 (MQ=255)
                                                                             gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATAAAAATCATGCATTTCCTGACACAACCGCCGCCGGGGATAACCCGCGGAGCGCCCTTCATTCTGTTTGTTGGTGGGAggga                                                                          >  2:799475/1‑144 (MQ=255)
                                                                               cgcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCCTACAAAACATGCATTTCCTTACCCAACAGACGCAGGGGATAAACCGCGTAGCGCTCTTTCTGCTGTTTGTAGGTGGAAGGGAGCAcct                                                                    >  1:1133195/1‑148 (MQ=255)
                                                                                  gTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCCTACAAAACATTTATTTTCTTTCCCAACCGTCGCAGGGGGTAAACGGCGGAGCGCCCTTTCTTCTGTTTTTAGGTGGAAGGGAGCAATTCt                                                                  >  1:613033/1‑149 (MQ=255)
                                                                                   tGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAAATATGCAATTTCTTACACAACTAACGCAGGGGGTAAACCGCGGAGGGCTCTTCACGATGTTTGTTGAGGGaaagaaacaattc                                                                   >  2:698538/1‑136 (MQ=255)
                                                                                    gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAAAAATGCATTTTCTTACACAAATGACGCAGGGGATAAACGGCGGAGCGCTCTTCATGCTGTTTGTAGATGGAAGGGAGCACTTCt                                                                  >  2:578082/1‑147 (MQ=255)
                                                                                     ggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCCTAAAAAAAATGCATTTTCTTTCACAACTTACGCAGGGGGTAAACCGCGGGGCGCCCTTCTTTTCGTTTGTTGATGGGAGGGAGCACtttt                                                                  >  1:82327/1‑144 (MQ=255)
                                                                                        aGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATAACAAACAGGTCTTTTCTTACACAAACAACGCAGGGGGTAAAAGGGGGAGGGCCCTTCATGTTGTTTTTGGGTGggagggagaaactttcgtaa                                                             >  1:536031/1‑136 (MQ=255)
                                                                                           gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCGTAACAAACATTCTTTTTCTTTCACCACTGTCGCCGGGGAGAACCGGCGGAGCGCCTCTCCTTCTCTTTTTTTATGGGAGGGAGCACCTCTCC‑TCAt                                                           >  1:718479/1‑146 (MQ=255)
                                                                                           gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAAACATTCATTTTCTTACACAAAAGACGCAGGGGAAAAACCGCGGAGAGCGCTTCATGCTGTTTTTAGAGGGGAGGGAGCACTTCTTC‑TAAttggg                                                       >  2:188607/1‑146 (MQ=255)
                                                                                            cTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATAAAAAAAATGCATTTTCTGAAACAAATGAAGCAGGGGATAAACGGCGGAGCGCCCTTTATTCTGTTTGTTGAGGGGAGGGAGGACTTCTTC‑TAATGGcgcg                                                     >  2:455986/1‑151 (MQ=255)
                                                                                             tCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCCTACAAAATATGCATTTTCTTGTCCAAAATAAGCAGGGGGTAAACGGCGTAGCGCTCTTCATTCTTTTTTTTGATGGGAGGGAGCAATTCTGC‑TTATgg                                                         >  1:1148871/1‑146 (MQ=255)
                                                                                              ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTc                                                   <  1:169862/151‑1 (MQ=255)
                                                                                                 tACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCCTAAAAAATATTTATTTTCTTTCACAACTTACCCAGGGGAGAAACGGCGTAGGGCTCTTTCTTTTGTTTTTTGATGGGAGGTAGAACCTTTTC‑TTATTGCGcggcggg                                                >  1:734277/1‑150 (MQ=255)
                                                                                                   cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCCTAAAAAAAATTTATTTTCTTTCACAACAAACGCGGGGGATAAACGGCGGGGCGCTCTTCTTTCTGTTTTTTGGTGGGAAGGAGGACTTTTGCTTAATgggcggcgg                                                  >  1:68193/1‑141 (MQ=255)
                                                                                                      aaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGAGCA‑ATAAAACATTCATTTTTTTTCACCACTGACGCCGGGGGAAAACGGCGGAGCGCCCTTCATTGTTTTTTTTGGAGGAAGGGAGCACATCTGC‑TAAtaggtagccggga                                               >  1:971571/1‑135 (MQ=255)
                                                                                                        cATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GGATACAAAACATTTATTTTCTTATACAACTTACGCAGGGGATAAACGGCGGAGCGCTCTTTTTGTTGTTTTTTGGTGGGAGGGAGCACTTTTGC‑TAATGGGGcggcggggt                                              >  1:104460/1‑143 (MQ=255)
                                                                                                          tctcTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATAAAAAATATTTTTTTTCTTGCACAACATACGCAGGGGAGAAACGGCGGAGGGCCCTTTATGTTGTTTTTTGATGGAAGGGAGCACTTCTTC‑TTATGGGGcggcggg                                                >  1:407950/1‑141 (MQ=255)
                                                                                                            tctTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAAACATTCATTTTCTTACACAAATGACGCAGGGGATAAACGGCGGAGCGCTCTTCCTGCTGTTTGTTGATGGAAAGGAGCACTTCTGC‑TAATGGCGCGTCGGGCTCAAGCAggt                                     >  2:1038067/1‑149 (MQ=255)
                                                                                                             ctTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAAAAATGCTTTTTCTTTCACCACCGACGCAGGGGGTAACCGGCGGAGCGGCCTTCTTGTTTTTTGTTGATGGAAGGGAGGAACTCTGC‑CCATggcggggcgggcgccagacggtt                                    >  1:416097/1‑138 (MQ=255)
                                                                                                               ttttGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGttt                                    <  2:531925/149‑1 (MQ=255)
                                                                                                                tttGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGttttt                                  <  2:536205/150‑1 (MQ=255)
                                                                                                                     aGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGc                            <  2:1053686/151‑1 (MQ=255)
                                                                                                                      gttgttTCACGTTGTGGCGG‑GGCTACAAAAAATGCCTTTTCTTTCACAACCGACGCCGGGGGTAAACGGCGGAGCGCCCTTCCTGCTGTTTTTTGATGGAAGGGAGCACTTCTTC‑TTATGGCGCGGCGGGCTCAAACAGgttttttaa                              >  1:972639/1‑144 (MQ=255)
                                                                                                                            tCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACca                      <  1:353885/150‑1 (MQ=255)
                                                                                                                                   gtgGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGca              <  1:772123/151‑1 (MQ=255)
                                                                                                                                       cGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca           >  1:839211/1‑150 (MQ=255)
                                                                                                                                         g‑gCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGt        <  2:1001592/151‑1 (MQ=255)
                                                                                                                                           gCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCACCAAGtt       <  1:980387/151‑1 (MQ=255)
                                                                                                                                           gCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGtt       >  2:205961/1‑151 (MQ=255)
                                                                                                                                             aTACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTg      >  2:900578/1‑150 (MQ=255)
                                                                                                                                             aTACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGCCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTg      >  2:205279/1‑150 (MQ=255)
                                                                                                                                               aCAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg   >  2:407966/1‑151 (MQ=255)
                                                                                                                                                cAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCt    <  1:499710/149‑1 (MQ=255)
                                                                                                                                                cAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg   >  2:374115/1‑150 (MQ=255)
                                                                                                                                                 aaaTCATGTATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGa  >  1:76218/1‑150 (MQ=255)
                                                                                                                                                  aaTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGa  >  2:77200/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGG‑GCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGC‑TAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGA  >  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: ‑