Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,457	0	T	G	41.9%	109.3 / 12.8	63	G232G (GGT→GGG)	dapF	diaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base T (10/26);  new base G (26/0);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.30e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGA  >  NC_000913/3995309‑3995599                                                                                                                                                     |                                                                                                                                                 aGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTTTGCAAGGgtt                                                                                                                                                >  2:633107/1‑151 (MQ=255)    gTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTGCAGAAGGGAGTCAGCAAGGGTTgt                                                                                                                                              >  2:402485/1‑149 (MQ=255)      tttCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTgg                                                                                                                                           >  2:1197533/1‑151 (MQ=255)      tttCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTg                                                                                                                                            <  2:922941/150‑1 (MQ=255)       ttCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTgg                                                                                                                                           <  2:877252/150‑1 (MQ=255)        tCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTTCTTGcc                                                                                                                                         >  1:449755/1‑151 (MQ=255)          cGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCga                                                                                                                                       <  2:1223672/151‑1 (MQ=255)            gAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGGAAGGAATCAGCAAGGGTTGCTgg                                                                                                                                           >  1:998551/1‑145 (MQ=255)                 cgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTAGGGATTCAGCAAGGTTTGCTGGCCgaagaa                                                                                                                                   >  2:560184/1‑148 (MQ=255)                      atatCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTCAGCAAGGGTTGCTGGCCgaagaa                                                                                                                                   >  2:1041358/1‑143 (MQ=255)                        atCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTgg                                                                                                                         <  1:140608/151‑1 (MQ=255)                          cGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCGGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGaa                                                                                                                       <  1:819872/151‑1 (MQ=255)                             ttttATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTTCTGGCCGAAGAAGGACGCGTGGAACTc                                                                                                                    >  2:67608/1‑151 (MQ=255)                                tATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGCGTAGGGATTCAGCAAGGGTTTGTGGCCGAAGAAGGACGCGTGGGACTcccc                                                                                                                 >  1:1042497/1‑151 (MQ=255)                                 aTGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCcg                                                                                                                <  2:38008/151‑1 (MQ=255)                                  tGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTGGAGGAGGGGTTTAACAAGGGTTTCTGGCCCGAGAAGGACGCGCGGGACTTCCcgg                                                                                                               >  2:676746/1‑151 (MQ=255)                                   gCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGTAGGGGGGGGTCTGCAAGGGTTGGTTGGCGGAGAAATACGCGGGGGAATcccccg                                                                                                               >  1:549819/1‑148 (MQ=255)                                   gCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCcgg                                                                                                               <  2:464725/150‑1 (MQ=255)                                    cAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGGAGGGAATAAGCCAGGGTTGCTGGCCGAAAGAGGACACGGGGAACTcccc                                                                                                                 >  1:43023/1‑147 (MQ=255)                                    cAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCATGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGAAGAAAGACGCGTGGAACTcccc                                                                                                                 >  1:781533/1‑147 (MQ=255)                                     aaGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAAATCCCccgcgg                                                                                                            >  2:258166/1‑151 (MQ=255)                                         ggTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGgtcg                                                                                                         <  2:171726/150‑1 (MQ=255)                                         ggTTAAGCGCGAGCAGATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGgtcg                                                                                                         <  1:218153/150‑1 (MQ=255)                                          gTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTGCTGGCCGAAGAAGTAAGCGGGGGACTCCCCCGCGggcg                                                                                                         >  2:644306/1‑149 (MQ=255)                                           ttAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAAGGTTGCTGGCCGAAGAAGGACGCGGGGGACTCCCCGGGGGGCGTCt                                                                                                      >  2:822643/1‑151 (MQ=255)                                            tAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCt                                                                                                      <  1:505036/150‑1 (MQ=255)                                             aaGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCtt                                                                                                     <  2:1016636/150‑1 (MQ=255)                                             aaGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTg                                                                                                    <  2:1161861/151‑1 (MQ=255)                                             aaGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTTAGCAAGGGTTTGTGGGCGGAGAAGGACGCGGGGGAATCCCCGGCGGGCGTCtt                                                                                                     >  2:177512/1‑150 (MQ=255)                                               gcgcGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTTAGCAAGGGTTTCTGGCCGGAGAAAAACGCGGGGGACTCCCCGGCGGGCGTCt                                                                                                      >  2:724283/1‑147 (MQ=255)                                                cgcgAGGATATTCGTTTACGCGTTTGTTGGCGTGGGGGAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCTGTCGTCTTGat                                                                                                  <  1:903102/150‑1 (MQ=255)                                                    aGCATATTCGTTTACGCGTTTATGAGCGTGGGGCCGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGCTGCCGGCGAGAGTCATCACGGAATGCTGGCCGGAGACGCAAGCGCGGTAGTTCCCGGAGGCAGCCTTGATATTGc                                                                                             >  1:686396/1‑151 (MQ=255)                                                        tatTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAACAAGGGTTGCTGGCCCAAGAAGTACGCGGGGAGCTCCCCCGCGGGCGTCCTGatattg                                                                                              >  1:681754/1‑144 (MQ=255)                                                                 aCGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTGCTGGCCGGAAAAGGACGCGGGGAACTCCCCGGCGGGCGGCTTGATATTGCCTGGGAAAggcc                                                                                 >  2:808423/1‑148 (MQ=255)                                                                  cgcgTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGTGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCgg                                                                              <  1:1172968/151‑1 (MQ=255)                                                                      tttATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAAGGTTGGTCGGCGGAGAAGTAAGCGTGGAACTCCCCCGCGGGCGGCTTTATATCGCCC‑GGGAAAGGCCggg                                                                             >  2:445678/1‑148 (MQ=255)                                                                         aTGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTAGGGGTTCCGCAAGGGTTGCTGGCCGAAGAAATAAGCGTGGAAATCCCCGGCGGGCGTCTTGATATCGCCT‑GGGAAGGGCCGGGGGA‑ccc                                                                      >  1:600838/1‑151 (MQ=255)                                                                         aTGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑cc                                                                       <  2:1225410/150‑1 (MQ=255)                                                                             gCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTGCTGGCCGGAAAAGTAAGCGTGGGACTCCCCGGCGGGCGGCTTTGTATCGCCC‑GGGAAAGGCCGGGGCA‑ccccgt                                                                   >  2:584396/1‑147 (MQ=255)                                                                              cGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGTATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAGCTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCC‑CCCGtt                                                                   <  1:1158322/149‑1 (MQ=255)                                                                                   ggCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATg                                                            <  1:364671/151‑1 (MQ=255)                                                                                       ggAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTg                                                        <  2:743534/151‑1 (MQ=255)                                                                                        gagaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGCAAGGGGTCAGCCAGGGGTGCTGGCCGGAGACGGAAGCGGGGGAACCCCCGGCGGGCGTCCTGATATCACCT‑GCAAAGGTGCGGGGCC‑CCCGTtatattttaatg                                                        >  1:89736/1‑143 (MQ=255)                                                                                         agaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTTCTGGCCGAAGAAGGACGCGGGGAAATCCCCGGCGGGCGGCTTGGTATCGCCC‑GGGAAGGGCCGGGGGA‑CCCGGtatat                                                               >  1:927395/1‑142 (MQ=255)                                                                                          gaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGcc                                                     <  1:1197178/151‑1 (MQ=255)                                                                                          gaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGcc                                                     <  1:812942/151‑1 (MQ=255)                                                                                               ccAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCcggcg                                                <  1:787256/151‑1 (MQ=255)                                                                                                   gCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAAGGTTTCTGGGCGAAGAAATACGCGTGGAACTCCCCGGCGGGCGTCTTGATATTGCCT‑GGGAAGGTTCGGGGCA‑CCCCTTATATATGACTTGCCCCGCGGGAc                                            >  2:774912/1‑151 (MQ=255)                                                                                                     cTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTTCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGGCCGGGGCA‑CCCGTTATATATTACTGGCCCGGCGGGACa                                           >  2:788424/1‑150 (MQ=255)                                                                                                        cGGCAGCGGCGCGTGTGCGGCGGGTGCAGTAGGGGTTCAGCAAAGGTTTCGCGCCCAGGAAGGACGCGCGGGGATCCCCCGCGGGCGGCTTGGTATCGCCC‑GGGCAGGGCCGGGGGACCCCGTTA‑GTATGACTGGgccccgggggaa                                           >  1:179108/1‑138 (MQ=255)                                                                                                             gcggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTacga                                  <  2:628048/151‑1 (MQ=255)                                                                                                              cggcgCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGGAGAAGGACGCGTGGAACTCCCCGGCGGGCGTCTTGGTATTGGCC‑GGGAAGGGCCGGGGCA‑CCCGCTGTATATGACTGGCCCGGGGGGGCATGTCTac                                    >  1:383678/1‑148 (MQ=255)                                                                                                                  gcgTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGAt                             <  1:86598/151‑1 (MQ=255)                                                                                                                          ggcggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGGTTTATTCAt                     >  2:763036/1‑151 (MQ=255)                                                                                                                            cggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCt                   >  1:165601/1‑151 (MQ=255)                                                                                                                             ggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGGCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCt                   >  2:780243/1‑150 (MQ=255)                                                                                                                              gTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTAt                 <  1:232460/151‑1 (MQ=255)                                                                                                                                   aGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGc            >  2:793873/1‑151 (MQ=255)                                                                                                                                      aGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAAc         >  2:769982/1‑151 (MQ=255)                                                                                                                                      aGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAAc         <  1:67596/151‑1 (MQ=255)                                                                                                                                        ggATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAAcc        >  2:52194/1‑150 (MQ=255)                                                                                                                                         gATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAg      <  2:512471/151‑1 (MQ=255)                                                                                                                                               gCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGga  >  2:406990/1‑149 (MQ=255)                                                                                                                                                cAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGCCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAg      >  1:887655/1‑144 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCT‑GGAAAGGTCCGGGTCA‑CCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGA  >  NC_000913/3995309‑3995599

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