Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,096,090	0	A	G	42.4%	94.3 / 19.3	67	I83V (ATT→GTT)	hisF	imidazole glycerol phosphate synthase, catalytic subunit with HisH
Reads supporting (aligned to +/- strand):  ref base A (8/30);  new base G (28/0);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.40e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT  >  NC_000913/2095942‑2096240                                                                                                                                                     |                                                                                                                                                        gCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGAt                                                                                                                                                        <  2:281638/150‑1 (MQ=255)    gCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGtt                                                                                                                                                       >  1:266030/1‑148 (MQ=255)          aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTg                                                                                                                                             <  1:650699/151‑1 (MQ=255)              cgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAg                                                                                                                                         <  1:61465/151‑1 (MQ=255)                  acgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAGTCTTTGGAAGa                                                                                                                                        >  2:706843/1‑148 (MQ=255)                    gctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAGTCTCTGGAAGATGcc                                                                                                                                    >  1:864914/1‑150 (MQ=255)                        aagaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGa                                                                                                                                <  2:298976/150‑1 (MQ=255)                         agaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAGTCTCTGGGAGATGCCGCGaa                                                                                                                               >  2:384122/1‑150 (MQ=255)                             ggCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTc                                                                                                                          <  1:255880/151‑1 (MQ=255)                             ggCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTc                                                                                                                          <  2:512906/151‑1 (MQ=255)                                   gACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTcc                                                                                                                     <  2:119723/150‑1 (MQ=255)                                    aCGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTT‑AAGTCTCTGGAAGAAGCCGCGAAAAATctttcttt                                                                                                                   >  1:83820/1‑151 (MQ=255)                                      gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAAGTTCTTGGAAG‑TGCCGCGGAAATtcttttcttt                                                                                                                  >  2:653525/1‑150 (MQ=255)                                      gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTg                                                                                                                 <  1:457131/151‑1 (MQ=255)                                      gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTg                                                                                                                 <  1:457150/151‑1 (MQ=255)                                       aaCTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGCTT‑AAGTCTCTGGAAGATGCCGCGACAATTCTTCCCTTTg                                                                                                                 <  1:75313/150‑1 (MQ=255)                                       aaCTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAAGCTCTGGGAAAAGCCCCGGAAAATCTTTCCttt                                                                                                                  >  1:161174/1‑149 (MQ=255)                                         cTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGAAT‑AAGTCTCTGCGAAAGACCCGGAAAATTC‑TTCCTTTTgcg                                                                                                              >  1:243894/1‑150 (MQ=255)                                           ggTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGgcgc                                                                                                             <  2:362934/150‑1 (MQ=255)                                            gtgtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTT‑AAGTTTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGCgg                                                                                                           >  2:623501/1‑151 (MQ=255)                                              gtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGAt                                                                                                         <  2:456922/151‑1 (MQ=255)                                                   aCGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAGTCTCTGGAAGATGCCGCGAAAAATTTTTTCTTTTGCGCGGGTaaaa                                                                                                     >  2:578662/1‑150 (MQ=255)                                                         tCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTGGCGCGGGTAAAATTTTCa                                                                                               >  2:911955/1‑150 (MQ=255)                                                         tCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTT‑AAGTCTCTGGAAGAAGCCGCGAAAATTCTTTTCTTTGGGGCGGATAAAATTTCCa                                                                                               >  1:426830/1‑150 (MQ=255)                                                          cACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCAt                                                                                              <  1:195320/150‑1 (MQ=255)                                                          cACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAGAGTTTGGGAGAAGCCGCGAAAATTCTTTTCTTTTGCGCGGGTAAAATTTcc                                                                                                >  1:898711/1‑148 (MQ=255)                                                            ccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGCGGGGTT‑AAGTATCTGGGAGATGCCGCGAAAATTTTTTTCTTTGGGGCGGACAAAATTTCCaata                                                                                            >  1:763187/1‑147 (MQ=255)                                                              gCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGCGGGGTT‑TAGTGTTTGGAAGATGGCGCGAAAATTTTTTTCTTTTGGGTGGGGAAAATTTTCCACAAct                                                                                         >  1:336584/1‑151 (MQ=255)                                                                ttCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAGTATCTGGAAGATGCCCCCTAAAATTTTTCCTTTGGAGCGGATAAAAttttc                                                                                                >  1:659723/1‑140 (MQ=255)                                                                  ccAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTcc                                                                                     <  1:697582/151‑1 (MQ=255)                                                                       gATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTT‑AAGTCTCTGGAAGATGCCCCGGAAATTTTTTTCCTTTGCGGGGGTAAAATTTTCATCAACTCTcc                                                                                     >  1:493539/1‑146 (MQ=255)                                                                            ccGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCa                                                                           <  1:879207/151‑1 (MQ=255)                                                                                 tgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTT‑TAGTCTCTGGGAGAAGCCGCGAAAATTTTTTTCTTTTGCGCGGGTAAAATTTTCATAAACACTCCCGCGCTGGCAGAccc                                                                      >  1:11385/1‑151 (MQ=255)                                                                                   gTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑TAGTCTCTGGAAAATGCCGCGAAAATTCTTTCCTTTGGCGCGGGTAAAATTTCCATCAAACCTCCTGCGCTGGGCGACCCa                                                                     >  1:605985/1‑150 (MQ=255)                                                                                   gTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAGTCTCTGGGAGATGCCGCGAAAAATCTTTTCTTTTGGGCGGGTAAAAATTTCATCAACTCTCCTTCGCTGGCAGAccc                                                                      >  1:674674/1‑149 (MQ=255)                                                                                   gTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTTGCGCGGGTAAAATTTTCATCAACTCTCCTGCGCTGGCAGAccc                                                                      >  2:387463/1‑149 (MQ=255)                                                                                        tAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAGTCTCTGGAAGATGCCGCGGAAATTTTTTTCTTTTGCGCGGGTAAAAATTTCATCAACTCTCCTGCGCTGGCAGAACCAACatt                                                                >  1:302950/1‑150 (MQ=255)                                                                                         aaaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAATCTCTGGAAGATGCCGCGAAAAATCTTTTCTTTTGCGCGGGTAAAATTTCCCTCAACTCTCCTGCGCCGGCAGACCCAACatt                                                                >  2:629072/1‑149 (MQ=255)                                                                                          aaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTT‑TAGTCTCTGGGAGATGCCCCGGAAATTTTTTTCTTTTGCGCGGGTAAAATTTCCATTAAACCTCCCGCGCCGGCAGACCCCACatt                                                                >  1:723739/1‑148 (MQ=255)                                                                                           aaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTT‑AAGTCTCTGGAAAATGCCGCGAAAATTCTTTTCTTTGGCGCGGATAAAAATTCCATCAACTCTCCTGCGCTGGCAGACCCAACattaaat                                                            >  1:509785/1‑149 (MQ=255)                                                                                           aaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATT‑AAGTCTCTGGGAGATGGCGCGAAAAATCTTTCCTTTTGGGCGGGTAAAAATTTCATCAAATCTCCTGCGCTGGCAGACCCAACattaatt                                                            >  2:698718/1‑151 (MQ=255)                                                                                                     tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCCgg                                                  >  1:7848/1‑151 (MQ=255)                                                                                                        cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTgg                                                  <  2:167845/148‑1 (MQ=255)                                                                                                        cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGAGGGGTT‑GAATGACTGGTAGATGGCGCGGAAAGTCTTTCCTTTTGGGGGGGTCAAACATCCATTAACACCCCCGCGCAGGCCGACCCCACCTTAATTACTCGCCCGGCCg                                               >  1:695600/1‑151 (MQ=255)                                                                                                             ggCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTTCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCg                                           >  2:1064867/1‑150 (MQ=255)                                                                                                                 gAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTg                                      <  1:859149/151‑1 (MQ=255)                                                                                                                   aGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCCCGAAAAATTTTTCCTTTGGGGCGGGTAAAAATTCCATCAACTCTCCTGCGCTGGCAGACCCCACATTAATTACTCGCCTGGCCGATCGCTTTgg                                     >  1:293800/1‑150 (MQ=255)                                                                                                                      gATCGACATTCCGTTTTGTGTGGCGGGGGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTTGCGTg                                 >  2:34112/1‑151 (MQ=255)                                                                                                                       aTCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTg                                 >  2:859388/1‑150 (MQ=255)                                                                                                                           aCATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAgtg                            <  1:669733/151‑1 (MQ=255)                                                                                                                           aCATTCCGTTTTGTGTGGCGGGGGGGGTT‑TAGTCTCTGGGAGATGCCGCGAAAATTTTTTTCTTTTGCGCGGATAAAATTTTCATCAACTCTCCTGCGCTGGCAGACCCCACATTAATTAATCGCCTGGCCGATCGCTTTTGCGTg                                 >  2:513181/1‑146 (MQ=255)                                                                                                                           aCATTCCGTTTTGTGTGGCGGGGGGGGTT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTTGCGCGGGTAAAATTTTCATCAACTCTCCTGCGCTGGCAGACCCCACATTAATTACTCGCCTGGCCGATCGCTTTGGCGGGCAgtg                            >  1:767132/1‑151 (MQ=255)                                                                                                                            cATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAgtgt                           <  2:761665/151‑1 (MQ=255)                                                                                                                                  gTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgtg                     <  1:998990/151‑1 (MQ=255)                                                                                                                                     ttGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGt                   <  2:307064/150‑1 (MQ=255)                                                                                                                                         gtgGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGt               <  1:599837/150‑1 (MQ=255)                                                                                                                                         gtgGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGt               <  1:384105/150‑1 (MQ=255)                                                                                                                                           ggCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTAt             <  1:387302/150‑1 (MQ=255)                                                                                                                                             cGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTg           >  1:435551/1‑150 (MQ=255)                                                                                                                                               ggTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATa        <  2:573640/151‑1 (MQ=255)                                                                                                                                                  gggATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATAcc      <  2:492109/150‑1 (MQ=255)                                                                                                                                                  gggATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATAcc      <  1:1063026/150‑1 (MQ=255)                                                                                                                                                  gggATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATAcc      <  1:979820/150‑1 (MQ=255)                                                                                                                                                  gggATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCt     >  1:773278/1‑151 (MQ=255)                                                                                                                                                  gggATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCt     <  2:234492/151‑1 (MQ=255)                                                                                                                                                  gggATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCt     <  2:305828/151‑1 (MQ=255)                                                                                                                                                   ggATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTg    <  2:945494/151‑1 (MQ=255)                                                                                                                                                     aTT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGt  >  2:614183/1‑151 (MQ=255)                                                                                                                                                     aTT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGt  >  2:67431/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                        GCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATT‑AAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT  >  NC_000913/2095942‑2096240

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