Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,270,455	0	G	T	42.1%	69.0 / 26.5	57	intergenic (+853/+170)	yhaC/garK	pentapetide repeats‑related protein/glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base G (8/25);  new base T (24/0);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.55e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTAT  >  NC_000913/3270305‑3270604                                                                                                                                                       |                                                                                                                                                        gTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTg                                                                                                                                                         <  2:141292/151‑1 (MQ=255)   tCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGct                                                                                                                                                       >  1:885070/1‑151 (MQ=255)          gggCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGtt                                                                                                                                                <  1:760364/151‑1 (MQ=255)            gCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGCTTgtcct                                                                                                                                                     >  2:1022554/1‑140 (MQ=255)               gTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCCCTTTTGCAc                                                                                                                                            >  1:689147/1‑150 (MQ=255)                   cTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTCtttccct                                                                                                                                           >  1:820229/1‑142 (MQ=255)                    ttATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTCACttg                                                                                                                                         >  2:97818/1‑146 (MQ=255)                     tATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACTGGTc                                                                                                                                       >  2:1011311/1‑149 (MQ=255)                     tATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTca                                                                                                                                      <  2:790821/150‑2 (MQ=255)                     tATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTc                                                                                                                                       <  1:252789/149‑1 (MQ=255)                       tGAACCACTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTGTCTTTGTTGCACTGGTCG‑Tgg                                                                                                                                  >  1:273802/1‑151 (MQ=255)                        gAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑Tggg                                                                                                                                 <  1:1118111/151‑1 (MQ=255)                         aaCCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGt                                                                                                                                <  1:238838/151‑1 (MQ=255)                             ccTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTTCACTGGTCG‑TGGGgtt                                                                                                                              >  2:360947/1‑149 (MQ=255)                                  ttGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGTCG‑TGGGGTTtcccccc                                                                                                                       >  1:440253/1‑151 (MQ=255)                                       cTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCCGTTGCCCTGCTCG‑TGGGTTTcccccccccgcg                                                                                                                  >  2:282014/1‑146 (MQ=255)                                         tGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑Gt                                                                                                                <  1:246009/150‑1 (MQ=255)                                          gCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTTCACCTTTTGCACTGGGCGCTTGGTTTTCCCCCCggggtt                                                                                                                 >  1:304667/1‑147 (MQ=255)                                            tCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTCTTTCGCTGGTCG‑TGGGGTTCCCCCCCCCGC‑Gttttc                                                                                                            >  2:193073/1‑148 (MQ=255)                                                   ggAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCAc                                                                                                      <  1:1052803/150‑1 (MQ=255)                                                   ggAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCAc                                                                                                      <  2:461555/150‑1 (MQ=255)                                                   ggAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGGTTTCTCTCTGTTTCACTGGTCG‑TGGGGGTTCCCCCCCCGC‑GTTACCTgggccc                                                                                                     >  2:157683/1‑147 (MQ=255)                                                     aGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCt                                                                                                   <  2:325534/151‑1 (MQ=255)                                                     aGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCt                                                                                                   <  2:73350/151‑1 (MQ=255)                                                        ttACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGGTTGCTCCCTTTTTCACCCGTCG‑CGGGGTTCCCCCCCCGGC‑GTTACCTGGCCCCCTGcc                                                                                                >  2:344296/1‑151 (MQ=255)                                                                 cACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACCGGGCG‑TGGGGTTCCCCCCCCGGC‑GGTACCTGGCACCCCGCCCCATGGAg                                                                                        >  1:1087135/1‑150 (MQ=255)                                                                  aCGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGGACAGGTCG‑GGGGGTTCCCCCCCCCGCGGTTACC‑GGCACCCCCCCCtct                                                                                            >  1:36934/1‑143 (MQ=255)                                                                  aCGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTTCACCGGTCC‑TGGGGTTCCCCCCCCCGC‑GGTACCCGGCcccccccccccttg                                                                                          >  2:613707/1‑137 (MQ=255)                                                                   cGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGccc                                                                                     <  2:885098/151‑1 (MQ=255)                                                                        tGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCCCTGTTGCCCTCGGCG‑TGGGGTTTCCCCCCAGGC‑CGTACCCCGCAACCCGCCCTATGGAGCCCGGAc                                                                                 >  1:174775/1‑150 (MQ=255)                                                                         gTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGAc                                                                                 <  2:982795/149‑1 (MQ=255)                                                                         gTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACt                                                                                <  2:15860/150‑1 (MQ=255)                                                                          ttACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTTTTTCACTGGTCG‑TGGGGTTCCCCCCCCGGC‑GGTTCCTGGCACCCCGCCCCATGGAGCCCGGGCttt                                                                              >  1:814971/1‑151 (MQ=255)                                                                          ttAACAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCGCTTGTCG‑TGGGGTTCCCCCCCCGGC‑GTTACCTGGCACCCCGCCCCATTGAGCCCGGACtt                                                                               >  2:238265/1‑150 (MQ=255)                                                                                gCCGCGCGGTGCGCTCTTACCGCGCCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTtc                                                                             <  2:81201/146‑1 (MQ=255)                                                                                gCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTcctccc                                                                        <  2:820254/151‑1 (MQ=255)                                                                                  cgcgcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACTGGTCG‑GGGGGTTTCCCCCCCCGC‑GTTACCTGGCCCCCCGCCCCATGGAGGCCCGGCTTTtcccccc                                                                       >  1:139943/1‑150 (MQ=255)                                                                                    cgcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTCCCCGGTCG‑TGGGGTTCCCCCCCCGGC‑GCTACCCCGCCCCCCCCCCCCTGGAGCCCGGACTTTcctccccccc                                                                    >  2:654129/1‑151 (MQ=255)                                                                                     gcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGTCG‑TGGGGTTTCCCCCCCCGC‑GGTTCCCGGCACCCTGCCCCATGGAGCCCCGACTTTcctccccccc                                                                    >  2:367696/1‑150 (MQ=255)                                                                                     gcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGCTTTCTCTCCTTTGCCCTTGTCG‑TGGGGTTCCCCCCCCCGC‑GCTCCCTTGCACCCCGCCCCCTGGCGCCCCGGCCTTccccccctccc                                                                   >  2:447961/1‑150 (MQ=255)                                                                                        gTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTCACCGGGCG‑TGGGTTTTCCCCCCCCGC‑GGTACCCGGCCCCCCGCCCTATGGAGCCCGGGCCTTTCTCcccccccccc                                                                >  1:966712/1‑151 (MQ=255)                                                                                          gcgcTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCg                                                               <  1:6391/150‑1 (MQ=255)                                                                                               ctTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGGCG‑TGGGGTTCCCCCCCCCGC‑GGTTCCTCGCCCCCCTCCCCCTGGAGCCCGGGCTTTTCTCCCCCCCCCCCCCCTccc                                                         >  2:748358/1‑151 (MQ=255)                                                                                                 tACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTcccc                                                        <  1:360941/150‑1 (MQ=255)                                                                                                  aCCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCg                                                      <  1:136960/151‑1 (MQ=255)                                                                                                   ccGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTccccc                                                       <  1:268946/149‑1 (MQ=255)                                                                                                     gCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGTCG‑TGGGGTTTCCCCCCCCGC‑GGTACCTGGCACCCCGCCCCATGGAGCCCGGACTTTTCTCCCCCCCCCCCCGCTccccc                                                       >  1:1033011/1‑147 (MQ=255)                                                                                                           tttCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTTCACTGGTCG‑TGGGGTTCCCCCCCCGGG‑GGTACCTGGCACCCCGCCCCATGGAGCCCCGGCTTTCCTCCCCCCCCCCCGGCTCCCCCCCAGAgg                                                >  1:1038738/1‑148 (MQ=255)                                                                                                            ttCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGgac                                              <  2:394387/149‑1 (MQ=255)                                                                                                                  cTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGGCA‑TGGGGTTCCCCCCCCGGC‑GGTTCCCGGCCCCCCGCCCCCTGGAGCCCGGGCCTTCCCCCCCCCCGCCCGTCTCCCCCCCAGAGgacgacgac                                        >  2:982236/1‑149 (MQ=255)                                                                                                                     aCCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAg                                    <  1:196672/150‑1 (MQ=255)                                                                                                                      ccTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAgcg                                  <  2:1087177/151‑1 (MQ=255)                                                                                                                          aTCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTTTTGCACTGGTCG‑TGGGGTTTCCCCCCCGGC‑GTTACCCGGCCCCCCTCCCCCTGGGGCCCCGGCTTTCCTCCCCTCCCCCCGCCTCCCCCCAAGAGGACGACGACGAAgcggc                                >  2:92548/1‑149 (MQ=255)                                                                                                                                     gcgGGCCATCGGCGGTTTTCTCTCTGTTGGACTGGTCG‑TGGGGTTTCCCCCCCGGC‑GGTACCCGGCACCCCGCCCCATGGAGCCCGGGCTTTCCTCCCCCCCGCCCGGCTCCCCCCAAGAGGACGACGACGAAGCGGCGGCCGTCttg                      >  2:264841/1‑146 (MQ=255)                                                                                                                                      cgGGCCATCGGCGGTTTTCTCTCTTTTGCACTGGTCG‑TGGGGTTCCCCCCCCGGC‑GGTACCTGGCACCCTGCCCTATGGAGCCCCGGCTTTCCTCCCCCCCCCCCCCCTCCCCCCCAGAGGGCGACGACGAAACGGCGACCGTCtcg                      >  2:155389/1‑145 (MQ=255)                                                                                                                                         gCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCtt               <  2:1085842/151‑1 (MQ=255)                                                                                                                                         gCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCtt               <  1:254588/151‑1 (MQ=255)                                                                                                                                          ccATCGGCGGTTTTCTCTCTTTTGCACTGGTCG‑TGGGGTTCCCCCCCCGGC‑GGTACCTTGCACCCCGCCCTATGGAGGCCGGGCTTTCCTCCCCCCCGCCCCTCTCCCCCCAAGAGGAGGACGACGAAGCGGCGACTGTCTGGGCAGc                 >  2:686452/1‑148 (MQ=255)                                                                                                                                                      tGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGta   <  2:814997/150‑1 (MQ=255)                                                                                                                                                       gCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGtat  >  1:1030543/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                        GTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCG‑TGGGTTTCCCCCCCAGGC‑GTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTAT  >  NC_000913/3270305‑3270604

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