Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	266,349	0	C	A	39.0%	126.7 / 20.4	59	V69F (GTT→TTT)	yafY	lipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (7/29);  new base A (23/0);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.74e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATG  >  NC_000913/266205‑266495                                                                                                                                                   |                                                                                                                                                    ccGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCaa                                                                                                                                               <  1:346119/151‑1 (MQ=255) ccGTCACCTCCTGGTATTCCAGTTTGATGTGATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCaa                                                                                                                                               <  1:649791/151‑1 (MQ=255)  cGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAg                                                                                                                                              <  2:605924/151‑1 (MQ=255)   gTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAA‑GT‑AAAAACTGCaa                                                                                                                                               >  1:703768/1‑149 (MQ=255)       cctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCCTCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTggg                                                                                                                                         <  1:95034/151‑1 (MQ=255)       cctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GT‑AAAAACCGCAAGGTggg                                                                                                                                         >  1:1161246/1‑151 (MQ=255)        ctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GG‑AAAACACGCAAgggg                                                                                                                                           >  2:205451/1‑146 (MQ=255)        ctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GT‑AAAAACCGCAAgggggg                                                                                                                                         >  1:322884/1‑150 (MQ=255)        ctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑AT‑AAAAACACCAAgggggg                                                                                                                                         >  2:117139/1‑150 (MQ=255)           ctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GG‑AAAAACTGCCAAGGGGGAcc                                                                                                                                      >  1:1080199/1‑150 (MQ=255)                 ttCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGA‑AAAACCTCCAAGG‑GGGGCCGGGTaa                                                                                                                                >  1:937657/1‑150 (MQ=255)                   ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAATGGCCGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAa‑gta                                                                                                                            <  2:302469/151‑1 (MQ=255)                       tttGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAAACTGCCAGGGGGGGCCGGaaaaag                                                                                                                              >  2:559177/1‑144 (MQ=255)                           aTGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCGGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCt                                                                                                                    <  2:1094801/151‑1 (MQ=255)                            tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAAACCGCCAGGGGGGGCCGGGTAA‑ATAGTTaa                                                                                                                       >  1:307040/1‑147 (MQ=255)                                aTCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GG‑AAAAACTGCAAAGTGGGGGCCGGTAA‑GGAGGCAACCttt                                                                                                                  >  2:713403/1‑146 (MQ=255)                                         acacAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCAttcttc                                                                                                       <  1:441017/150‑1 (MQ=255)                                         acacAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCg                                                                                                      <  2:358994/151‑1 (MQ=255)                                          cacaGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GG‑AAAAACTGCAAGGGGGGGCCGGGTAA‑AGAGTCAACCTGTGACCAttctt                                                                                                        >  1:558688/1‑148 (MQ=255)                                                        caGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGGC‑GT‑GAAAACCGCGAGGGGGGAGCCGCAAA‑GGAGTCAAACTGGGAACATTTTTTGCTGATAACATCCttg                                                                                       >  1:221857/1‑149 (MQ=255)                                                            ggCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GT‑AAAAACCCCGAGGGGGGAGCGGGTAA‑AGAGTACAAATGGGACCAATCTTTCCTGGCATCAACCcggggg                                                                                    >  1:542947/1‑150 (MQ=255)                                                               tGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAA‑GT‑AAAAACCGCCAAGCGGGGCCGGAGAA‑GTAGTGGAACTGTGGACAATCTTCGCTTACATCATCCttgtgg                                                                                    >  1:94259/1‑147 (MQ=255)                                                                ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTg                                                                                 <  2:331873/149‑1 (MQ=255)                                                                  aCGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACca                                                                             <  2:846462/151‑1 (MQ=255)                                                                  aCGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GT‑AAAAACCCCAAGGTGGGGGCGGGTAA‑GAAGGGAAACTGTGACCATTTTTCCCTGCCCTCATTCTGGGGGGTGGcc                                                                              >  1:10524/1‑150 (MQ=255)                                                                    ggTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATTGGGCGACTGGCGGGGAA‑GG‑AAAAAATGCCAGGGGGGGCCGGGAAA‑GAAGGTAACCTGTGACCAATCTTTGCCGACATCCTCCCGGTGGGTGACca                                                                             >  1:67888/1‑149 (MQ=255)                                                                      ttttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GG‑AAAAACTGCAAGGGGGGGCCGGATAA‑ATAGTCAAACCGTGGCCCATTTTTGCTGACATCATCCTGGGGGGTGagc                                                                              >  2:957334/1‑144 (MQ=255)                                                                       tttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GT‑AAAAACCGCAAAGGGGGGCCGGATAA‑ATAGGTAAACCGTGACCATTTTTCGCTGACATCATCCTGGGGGGTGAcc                                                                              >  1:927007/1‑145 (MQ=255)                                                                        ttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GTGAAAACCGCCAAGG‑GGGGCCGGGGAA‑AAAGTCAAACCGTGACCATTCCTCTCTGACAGCATCCTTGTGGGTGGCcaca                                                                           >  2:668570/1‑147 (MQ=255)                                                                         tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAc                                                                       <  1:404845/150‑1 (MQ=255)                                                                          ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACt                                                                      <  1:1022057/150‑1 (MQ=255)                                                                           tCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTg                                                                     <  1:211918/150‑1 (MQ=255)                                                                                cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgc                                                               <  2:693339/151‑1 (MQ=255)                                                                                cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GT‑AAAAACCGCAAGGGGGGGCCGGGTAA‑GTAGGCAAACCGTGACCATTTTTTGCTGACATCATCCTGGGGGGTGACCCCATCACTGaaa                                                                  >  1:1054556/1‑148 (MQ=255)                                                                                cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GT‑AAAAACCGCAAAGGGGGGTCGGGTAA‑GTAGGTCACCCGTTGCCATTCTTCTCTGTCATTATCCTGGGGGCGGACCCAAACACTGa                                                                    >  2:312325/1‑146 (MQ=255)                                                                                   atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGt                                                            <  2:109766/151‑1 (MQ=255)                                                                                     ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAA‑GG‑AAAAACCGCAAAGGGGGGCCGGATAA‑ATAGTTCACCCGTGACCATTCTTCCCTGGCATCCTCCTGGTGGCGGGCCACCTCACTGAAAAGCGgggg                                                          >  1:1169151/1‑151 (MQ=255)                                                                                      cACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTggg                                                         <  1:821490/151‑1 (MQ=255)                                                                                        cGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt                                                       <  1:836702/151‑1 (MQ=255)                                                                                        cGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGAAC‑GG‑TAAACACGCCAGGTGGGGGCGGGTAAGGAAGCCAACCT‑TGAACCCTCTTCCCTGAACACATCCTGGTGGGCGGCCACAACACCGAAAAGCGGAgggg                                                        >  1:634124/1‑150 (MQ=255)                                                                                            gtagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAgg                                                   <  1:776870/151‑1 (MQ=255)                                                                                                 agtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCtt                                              <  1:156955/151‑1 (MQ=255)                                                                                                 agtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCt                                               <  1:606052/150‑1 (MQ=255)                                                                                                      cTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTTC‑GT‑AAAAACCGCAAGGGGGGACCGGATAA‑ATAGGGAACCCGTGACCCTTTTTCGCTGACATCATCCTGGTGGCTGACCACATCACCGGAAAGCGgggggggggggggccttaa                                           >  1:366042/1‑136 (MQ=255)                                                                                                      cTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTAGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAAt                                          <  2:593098/150‑1 (MQ=255)                                                                                                       tGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAC‑GG‑AAAAACCGCAAGGGGGGGCCGGGTAA‑AGAGTCAACCCGTGACCATTCTTCGCTGACCTCATCCTTGTGGCTGACCCCAACACTGAAACGCGGGGGGGGGGGGGCCTTAAAttt                                        >  2:483331/1‑151 (MQ=255)                                                                                                        gCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAAtt                                         <  1:78281/149‑1 (MQ=255)                                                                                                         ccATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTAc                                      <  2:937609/151‑1 (MQ=255)                                                                                                              cTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGa                                   <  2:806749/149‑1 (MQ=255)                                                                                                               gcGCAGCCATCGGGCGACTGGCGGGGAC‑GT‑AAAAACTTCCAGGGGGGGCCGGATAA‑GAAGTTAACCCGTGACCCTTCTTCTCTTACATCATCCTTGTGGGTGACCACATCACTGAAAAGCGGGGGGGGGGGGGGCTTTAAttt                                        >  2:217462/2‑143 (MQ=255)                                                                                                                         cGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGt                       <  1:1096403/150‑1 (MQ=255)                                                                                                                          gggCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtc                      <  2:51352/150‑1 (MQ=255)                                                                                                                          gggCGACTGGCGGGGAC‑GG‑AAAAACTGCAAGGGGGGGCCGGGTAA‑GGAGTCAACCTGTGGCCATTTTTCGCTGACATCATCCTGGTGGGTGACCACATCACTGAAACGCGGGGGGGGGGGGTCCTTAAATTTACTGAAGGCCTTCAgggc                      >  2:184364/1‑148 (MQ=255)                                                                                                                           ggCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcat                    <  1:50228/151‑1 (MQ=255)                                                                                                                                 tGGCGGGGAC‑GT‑AAAACCTGCAAGGGGGGAACGGATAA‑GTAGTCAACCTGTGACCATTTTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGGGGGGTGGGGTCCTTAAATTTACTGATGGCCGTCAGGTCATCACtt               >  2:162696/1‑150 (MQ=255)                                                                                                                                  ggCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc              <  1:396050/150‑1 (MQ=255)                                                                                                                                            gT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGc     >  2:977276/1‑150 (MQ=255)                                                                                                                                               taaaCCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATg  >  2:1169872/2‑151 (MQ=255)                                                                                                                                               aaaaCCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAt   <  1:939013/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAC‑GT‑AAAACCTGCAAGGTGGGACCGGATAA‑GTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATG  >  NC_000913/266205‑266495

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: ‑