Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,516	0	C	A	50.0%	52.9 / 48.9	66	Q295K (CAA→AAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (30/3);  new base A (0/33);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.98e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436367‑436666                                                                                                                                                       |                                                                                                                                                       tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                        >  2:427846/1‑150 (MQ=255)             ttGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgccg                                                                                                                                            >  1:361427/1‑150 (MQ=255)                 aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATc                                                                                                                                       >  1:27144/1‑151 (MQ=255)                    ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGa                                                                                                                                     >  1:10497/1‑150 (MQ=255)                       aCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGCGTTGGGTTTCACTGTGCCGG‑GACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGGACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGa                                                                                                                                     >  2:21182/1‑147 (MQ=255)                       aCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAggg                                                                                                                                 >  2:496209/1‑151 (MQ=255)                       aCAGGCGCTGCGCGGGGCGCTCTCTGGCGGTGAAGATTACGATTTGTTTTTCACTGGGCGGG‑AACTGACCCGGGGCGCGGGGGTTGTGGCTCTCGGCCCCTGGGGCGTACCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAggg                                                                                                                                 <  1:962049/151‑1 (MQ=255)                            cgctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTg                                                                                                                            >  2:140361/1‑151 (MQ=255)                             gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCttt                                                                                                                             >  2:663461/1‑149 (MQ=255)                             gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                           >  1:32596/1‑151 (MQ=255)                              ctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                          >  2:534649/1‑151 (MQ=255)                                 cgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                        >  1:775375/1‑150 (MQ=255)                                    gggcccctctcgGGGGGGAAAGTTTAAGGGTTGTGTTTCACTGTGCCGG‑AACTGACCGGGGTCGCGTGGGATGTGGTTCCCGGCCCCTGGGGCGTACTTTTTCCTTTTTCCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                     <  2:745572/144‑1 (MQ=255)                                     gggccttctcgGGGGGTAAAGATTCCGTTTTGTTTTTCAGTG‑GCCGGTAACGGACCCGTGGCCCGTGGGTTGTGGTTCTGGGACCCCGGGGCGTCCTTTTTCTTTTTATGGGGAAAATTTCCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                      <  1:1138635/142‑1 (MQ=255)                                     gggCGTTCTCGGGGGGTGAAGATTACGAGTTGTTTTTCACTGTGCGGG‑AACTGACCCGTGGCGCGGTGGTTGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                                    <  1:1115374/150‑1 (MQ=255)                                     gggCGCTTTCGGGGGGTGAAGATTACGAGTTGTGTTTCACTGTCCCGG‑AACTGAACGGTGGCGCGCTGGATGTGGCTCTCGGCCCCCGGGGCGTCCCTTTTCCCTTTTTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                                    <  1:504838/150‑1 (MQ=255)                                     gggCGCTCTGTGGCGGTGAAGATGCCGTGTTGTTTTTCACTGTGCCGA‑AACGGACCCGGGGCGCGTGGGTTGTGGCTCTCGGACCCGTGGGCGTCCCGTTTACCTGTGTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                     <  1:1147718/149‑1 (MQ=255)                                         gcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc                                                                                                               >  1:1100757/1‑151 (MQ=255)                                           tctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                             >  2:329601/1‑151 (MQ=255)                                                ggCGGTGAAGATTACGGGTTGTTTTTCCCTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTCCCGTTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                          <  1:833337/149‑1 (MQ=255)                                                ggCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAc                                                                                                        >  2:493291/1‑151 (MQ=255)                                                   ggggaggatgaggggttgttTTTCGGTGTGGCGG‑AGCTGCGGCGTGGCGGGGTGGTTGTGGTTCGCGGACCCGGGGGCTTGCTTTTTTCCTGGCTGGGAAAAGTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                     <  2:434922/137‑1 (MQ=255)                                                     tGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                                   >  1:27346/1‑151 (MQ=255)                                                      gAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                                   >  1:763558/1‑150 (MQ=255)                                                       aaGTTTACGAGTTGTGTTTCCCTGTGCCGG‑ACCTGAACCGTGGCGCGCGGGATGTGGCTCTCGGCCACCGGGGCGTACCTTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                                 <  2:662054/151‑1 (MQ=255)                                                       aaGATTACGTTTTGTGTTTCACTGTGCGGG‑AACTGACCCGTGGCGCGCTGGATGTGCCTCCCGGCCCCCGGGGCGTACTTTTTCCTTTTACCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                                 <  2:265156/151‑1 (MQ=255)                                                        aGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                                 >  1:907685/1‑150 (MQ=255)                                                               gtgttgttTTTCACTGGGCGGG‑ACCGGACCCGGGGCGCGTTGGAGGTGGCCCGCGCCCACGTGGGCGTCCTGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                                             <  1:979721/145‑1 (MQ=255)                                                               gggttgttTTTCCCTGTGCGGG‑AACTGAACGGTGGCGCGGTGGATGGGGCTCCCGGACACCTGGGCGTCCCTTTTCCTTGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                                             <  1:71441/145‑1 (MQ=255)                                                                  ttttgtTTCATTGGCCGGG‑ACCTGACCGGTGGCGCGGAGATGGTGTCTCCCGGCCCCTGGGGTGTACTTTTTCCCTGTTTCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                        <  1:1125050/146‑1 (MQ=255)                                                                  ttGTTTTTCACTGTGCCGA‑AACTGAACCGTGGCGCGCTGGATGTGTCTCTCGGACCCCGGGGCGTCCCGTTTCCTTGTAGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                         <  1:846163/148‑1 (MQ=255)                                                                  ttGTGTTTCACTGTGCGGG‑AACTGAACCGGGGCGCGCTGGATGGGGCTCTCGGACACCTGGGCTTCCCGTTTCCCTGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                        <  1:145722/149‑1 (MQ=255)                                                                   tgtgTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                      >  2:562168/1‑150 (MQ=255)                                                                   tgtgTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                      >  2:1090075/1‑150 (MQ=255)                                                                    gtgtTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGCTGTGGCTCTCGGACACCTCGGCGCACCGTTTACCTGTATCGGGCACATGCCCGCCGACATCGACGGGCTTTGTTTTATTCGTGACGGCGAACCTTTTACATTAGACTgg                                                                                       >  2:46286/1‑148 (MQ=255)                                                                     tgtTTCCCTGGGCGGA‑AACGGCCCGGTGGCGCGCTGGTTGGGGCTCGCGAACACCTGGCTTTGCCGTTTACCTTTAGCGGGAAAAAGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                                                    <  2:327057/150‑1 (MQ=255)                                                                     tgtTTCACTGTGCCGG‑AACTGAACGGGGGCGCGCTGGGTGGGGCTCTGGGACCGGGGGGCGTCCCTTTTACCTGTACCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                     <  1:583919/149‑1 (MQ=255)                                                                       tttCACTGTGCCGA‑AACTAACCCGTGGCGCGTGGGATGTGGCTCCCGGCCCCCGGGGCTTCCCTTTTCCCTTTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                                  <  1:739628/150‑1 (MQ=255)                                                                          ccttgtgCCGA‑AACTGACCCGTGGCGCGCTGGATGTGGCTCCCGGCCCCCGGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                               <  1:417951/147‑1 (MQ=255)                                                                                 ccGG‑AACTGACCCGTGGCGCGTTGGTTGTGGCTCTGGGACCCTTGGGCTTCCCTTTTCCCTGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                                          <  1:211032/148‑1 (MQ=255)                                                                                   gg‑AACTGACCCGTGGCGCGCTGGGTGTGCCCCGCGGCCCCCGGGGGGTACTTTTTACCTGTCGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                     <  2:135611/151‑1 (MQ=255)                                                                                    gaaaCTGAACGGTGGCGCGGTGGTTGTGTCTCTGGGCCCCCGGGGGGTCCATTTTTCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                                <  1:862186/138‑1 (MQ=255)                                                                                          gacccGGGGCCCGTTGGATGTGACTCTCGCTCCCCGGGGGGACCCGTTTCCTTTTTTGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCa                                                                         <  2:969409/138‑1 (MQ=255)                                                                                             ccGGGGCGCGCTGGATGTGGCTCTCGCCCACCGGGGCGTACCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                             <  1:671200/150‑1 (MQ=255)                                                                                                        tGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                                  >  2:665389/1‑150 (MQ=255)                                                                                                         ggtggtgGCTCTCGAACACGGGGGCTTCCCTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                                  <  2:629088/145‑1 (MQ=255)                                                                                                         ggATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                                 >  1:761210/1‑150 (MQ=255)                                                                                                          gttggggCTCGCGGACCCCGGGGCGTCCCGTTTCCTTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                                  <  1:842585/143‑1 (MQ=255)                                                                                                               gcctctCGGCCCCTGGGGCGTCCCTTTTCCCTTTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaaa                                           <  1:1060714/148‑1 (MQ=255)                                                                                                                 ctctCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                        >  2:416644/1‑151 (MQ=255)                                                                                                                 ctccgggACCCCGGGGCGTCCGGTTTCCTTTTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAg                                          <  1:19130/144‑1 (MQ=255)                                                                                                                   ctgggacccgggggCGTCCTTTTTCCTTGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaaa                                           <  1:150012/135‑1 (MQ=255)                                                                                                                   ctcGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc                                      >  2:11655/1‑151 (MQ=255)                                                                                                                   cccgcacacgggggCGTACCTTTTCCCTGTATCGGAAAAGTCACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                       <  2:100340/139‑1 (MQ=255)                                                                                                                            cTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                             >  2:590043/1‑151 (MQ=255)                                                                                                                             ggggcttccggTTTCCCTGTGTGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAGCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGAt                              <  1:58787/139‑1 (MQ=255)                                                                                                                              gggCGTCCCTTTTCCCTGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGAt                              <  1:847692/148‑1 (MQ=255)                                                                                                                              gggCGTACCTTTTCCTTGTAGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGt                                                   <  2:495453/127‑1 (MQ=255)                                                                                                                              gggCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                             >  1:614098/1‑149 (MQ=255)                                                                                                                                     ctgtttccttttATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                     <  1:985022/139‑1 (MQ=255)                                                                                                                                       gTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg                  >  1:911987/1‑151 (MQ=255)                                                                                                                                        tttACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTgg                 >  1:1076184/1‑151 (MQ=255)                                                                                                                                          tACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc                >  2:610988/1‑150 (MQ=255)                                                                                                                                            ccTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATc             >  1:477206/1‑151 (MQ=255)                                                                                                                                              tGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt            <  1:427955/150‑1 (MQ=255)                                                                                                                                                tATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt         >  1:756486/1‑151 (MQ=255)                                                                                                                                                    gggCAATTGCCCGCCGATATCGAATGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACAATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc                <  2:922642/140‑1 (MQ=255)                                                                                                                                                       cAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTc  >  1:430839/1‑151 (MQ=255)                                                                                                                                                       cAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTc  >  2:935917/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436367‑436666

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