Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,228,222	0	T	G	53.8%	55.1 / 78.4	54	pseudogene (388/1287 nt)	mdtQ	pseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base T (3/21);  new base G (28/0);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 1.05e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.79e-10

GAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCG‑AGATG‑CCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCC  >  NC_000913/2228074‑2228367                                                                                                                                                     |                                                                                                                                                       gaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGgtt                                                                                                                                                      >  2:662709/1‑151 (MQ=255) gaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTAAACTGCCAGCCAGCAATTGGCGGGttt                                                                                                                                                      >  2:1180542/1‑151 (MQ=255)  aagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTg                                                                                                                                                     <  2:680918/151‑1 (MQ=255)   agTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGgtt                                                                                                                                                      >  2:366451/1‑149 (MQ=255)   agTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTg                                                                                                                                                     <  2:570644/150‑1 (MQ=255)         cccGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTcgcg                                                                                                                                               <  2:223987/150‑1 (MQ=255)          ccGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGgtt                                                                                                                                                      >  1:626501/1‑142 (MQ=255)           cGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGtt                                                                                                                                             >  2:517070/1‑150 (MQ=255)            gTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGtt                                                                                                                                             >  2:255063/1‑149 (MQ=255)                tGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCACCACTTTGCGGGATTggtttcgcgccagc                                                                                                                                        >  2:14577/1‑137 (MQ=255)                 gATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGTTcccgc                                                                                                                                        >  2:85122/1‑146 (MQ=255)                  atatAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGCTCCgccg‑cc                                                                                                                                    >  1:633826/1‑151 (MQ=255)                     taGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTtctctgggttcccccccc                                                                                                                                    >  1:602799/1‑135 (MQ=255)                      aGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGTCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTg                                                                                                                               <  1:485992/151‑1 (MQ=255)                      aGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATCGGCGGGGTTTCTCGCCTTTCGccc‑ccc                                                                                                                                   >  2:712626/1‑148 (MQ=255)                         tGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCGCCG‑CCC‑CTGG‑tt                                                                                                                           >  2:325917/1‑151 (MQ=255)                          ggCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑tt                                                                                                                           <  2:923156/150‑1 (MQ=255)                             gATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAAc                                                                                                                       <  2:366569/151‑1 (MQ=255)                                cGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGt                                                                                                                    <  2:1127488/151‑1 (MQ=255)                                    cacgATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCGCTCCCAGTACAGGCGGGCTACGCTGCCTGCCAGCAATGGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTAc                                                                                                                  <  1:13337/147‑1 (MQ=255)                                     gcgATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCG‑GCC‑GCGC‑TTTaaacctaacc                                                                                                                >  2:327116/1‑142 (MQ=255)                                           aTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCggg                                                                                                         <  1:581969/151‑1 (MQ=255)                                           aTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCCCCG‑CCC‑CTGG‑TTTTACCCTACCCCGGCggg                                                                                                         >  2:995212/1‑151 (MQ=255)                                                gtTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGt                                                                                                     <  1:899557/150‑1 (MQ=255)                                                     gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCAGCCAGCCAACAATTGGCGGGGCTTATCCCGGTGCTCCC‑GCC‑CTGC‑TTTAAACATTCCCCGGCggggtggacac                                                                                                 >  1:798523/1‑141 (MQ=255)                                                           tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTACCAGCCAGCAATTGGCGGGGTTTCTCGCGCTCCGCCT‑CCC‑GTGC‑TTTTAACTTTCCCCGGCGGGCCGGAATCTcccc                                                                                            >  2:757850/1‑146 (MQ=255)                                                           tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGATAAGCTGCCCGCCAGCAATTGGCGGGGTTACGCGCGGTTCCCCG‑CCC‑CTGG‑TTTTAACCTTCCCAGGCGGGGCGTTATCCCCGccc                                                                                          >  1:1304619/1‑150 (MQ=255)                                                            ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGTAGTAACCTCCGCCCg                                                                                         <  2:1153777/150‑1 (MQ=255)                                                            ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTTCGCCG‑CCC‑GGGG‑TTTAAACGGACCCAAGGCGGCCGTAAACTCCGccc                                                                                          >  2:1185075/1‑149 (MQ=255)                                                               tATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCG‑CCC‑GCGC‑TTTAACCCGACCCAGGCGGGGCGTAACCTCCGCCCGCtt                                                                                      >  2:702403/1‑150 (MQ=255)                                                                 tttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCTGCAATTGGCGGGGTTGGTTCCGGTCCGCCC‑CCC‑CAGC‑TTTATCCCCACCCCGGGGGGCGGGAACCTcccccccaat                                                                                      >  1:483326/1‑144 (MQ=255)                                                                 tttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGGTCCGCCG‑CCC‑GTGC‑TTTTACCCTACCCAGGCGGGGAGGAAACCCCCCCCCATtttt                                                                                   >  2:650691/1‑151 (MQ=255)                                                                  ttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCtt                                                                                   <  1:755876/150‑1 (MQ=255)                                                                      tgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAAATGGCGGGGTTGCTTGCGTTTCCCCG‑CCC‑GAGA‑TTTTACCGCAACCACGGGGGTCTTAAACCCCGCCCCATTTCTTcccc                                                                              >  2:1132434/1‑151 (MQ=255)                                                                        cAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAg                                                                            <  1:933116/151‑1 (MQ=255)                                                                           gACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGata                                                                         <  2:1003102/151‑1 (MQ=255)                                                                            aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGatat                                                                        <  1:1221049/151‑1 (MQ=255)                                                                               gtgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTTAACGTACCCAGGCGGGCCGTAACCCCCGCCCGATTTTTTtccccg                                                                            >  2:503494/1‑142 (MQ=255)                                                                                tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCCCGCGGTTCCCCG‑CCCGGGGC‑TTAAACCG‑ACCCAGGCGGGCAGCAACCTCCCCCCCATTCTTTCCCCCGATATc                                                                       >  2:654127/1‑148 (MQ=255)                                                                                    tAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCG‑AGATg                                                                <  2:458198/150‑1 (MQ=255)                                                                                     aaCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCG‑AGATG‑cc                                                             <  2:1201879/151‑1 (MQ=255)                                                                                     aaCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCCCCGCCCC‑CGGC‑TTTAACCG‑ACCCCGGCGGGGAGGAACCCCCCCCCGCGTTTTTTCCCCCATATCG‑AGATggc                                                              >  1:1131907/1‑149 (MQ=255)                                                                                      aCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTCTTTGCGGACCGCCG‑CCC‑CTGC‑TTTAAACCCAACCAGGCGGGGCGGAAACTCCCCCCGATTCTTTCCCCAGaaatt                                                                       >  2:1218601/1‑141 (MQ=255)                                                                                      aCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTCCCAGCCAGCAATTGGAGGGGTTCCCCGCGGTTCCCCC‑CCC‑CTGC‑TTTAACCCTTTCCCGGCGAGCAGGAAACTCCGCCCCCTTTTTTTCCCCGata                                                                         >  2:883665/1‑140 (MQ=255)                                                                                       cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCG‑AGATG‑CCAg                                                           <  1:1215368/151‑1 (MQ=255)                                                                                        gccgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCG‑CCC‑CGGCTTTTACCCGTCCCCCGG‑GGGCAGTAAACCCCCCCCCCTTCTTTTCCCCGAAATCGAAGATG‑c                                                              >  1:1198314/1‑148 (MQ=255)                                                                                                      ccACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAAACGTACCCAGGCGGGCAGTAAACCCCGCCCGATTCTTTTCCCAGATATTG‑AGATG‑CCAGCCCCCCGtttaacc                                             >  2:892702/1‑145 (MQ=255)                                                                                                              aGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGGTTCGCCGCCCC‑GGCC‑TTTAACCCCCCCCAGCGGGGACGAAACCCCCCCCCGTTTTTTTCCCCA‑AAATCG‑CGAAGACCGCCCCCCCCTTAACCCAAAAgacc                                     >  1:961771/1‑148 (MQ=255)                                                                                                                tACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGGTTCCCCC‑CCC‑GGGC‑CTTTACCCTACCCAGGCGGGCAGTAAACTCCCCCCCGTTTTTTTCCCCGATATCC‑AGATG‑CCCGCCCCCCCTTAAACCaaaaaaaccc                                   >  2:955798/1‑144 (MQ=255)                                                                                                                       gggCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCG‑GGATG‑CCAGCTCGCCGTTAAGCCAAAAGTACCGTTGgtg                             >  2:731780/1‑149 (MQ=255)                                                                                                                          cTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCCCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAAACTCCGCCCGATTCTTTCCCCAGATATCC‑AGATG‑TCAGCCCGCCCTTAAGCCAAAAGTACCCTTGGTGTTcc                         >  1:1036134/1‑150 (MQ=255)                                                                                                                              gCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCCCCG‑CCC‑GTGC‑TTTAAACGCACCCAGGCCGGCGGGAAACTCCCCCCCATTCTTTCCCCAAGTTTCA‑AGGTG‑CCCGCCCGCCGTTTAGCCAAAAATACCATTGGGATACCACcgg                    >  1:1148454/1‑151 (MQ=255)                                                                                                                                  ccagccagCAATTGGCGGGTTTCCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTCCCCAGGGGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCG‑AGATG‑CCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                  <  2:490209/149‑1 (MQ=255)                                                                                                                                    agccagcAATTGGCGGGGTTTCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCCATTCTTTCCCCAGATATCG‑AGATG‑CCAGCCCGCCGGTAAGCCAAAAGTACCGGTGGTGTACCACGGGCCgg                >  1:1207138/1‑149 (MQ=255)                                                                                                                                    agccagcAATTGGCGGGCTTGCTCGCGTTGCGCCG‑CCC‑GGGG‑TTTTACCCGACCCCGGGGGGCAGGAAACCCCCCCCCATTTTTTCCCCAGATATCG‑AGATT‑CCAGCCCGCCGTTAAGCCAAAAAGACCGATGGCGGACCACggggc                  >  1:136596/1‑145 (MQ=255)                                                                                                                                      ccagcAATTGGCGGGGTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTTACCGGACCCACGCGGGGAGTAACCTCCGCCCGATTCTTTCCCCAGAAATCG‑AGATG‑CCAGCCCGCCCTTAAGCCCAAAGTACCGTTGGTGTACCCCCGGCCGGGCg             >  1:149491/1‑150 (MQ=255)                                                                                                                                           aaTGGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCG‑AGATG‑CCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTg       <  2:829563/151‑1 (MQ=255)                                                                                                                                                gCGGTTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTACCCGTCCCCAGGCGGGCAGTAACCCCCGCCCGATTCTTTCCCCAGATATCG‑AGATG‑CCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGcc  <  2:316800/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCG‑CCC‑GTGC‑TTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCG‑AGATG‑CCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCC  >  NC_000913/2228074‑2228367

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