Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,148,330	0	T	G	44.6%	108.4 / 14.5	65	G219G (GGT→GGG)	yegD	Hsp70 chaperone family protein
Reads supporting (aligned to +/- strand):  ref base T (6/30);  new base G (29/0);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 6.53e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCAAGGGATTCTGGAACGCGCGGCGAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGTTTAAAAACCTGATGCCACTGCTG  >  NC_000913/2148186‑2148477                                                                                                                                                 |                                                                                                                                                    ggCGCAAGGGATTCTGGAACGCGCGGCGAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACgact                                                                                                                                               <  1:516548/151‑1 (MQ=255)             tCTGGAACGCGCGGCGAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGttttct                                                                                                                                    >  1:285121/1‑146 (MQ=255)               tGGAACGCGCGGCGAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTAAGACTGACTGTTTAtt                                                                                                                                    >  2:922785/1‑148 (MQ=255)                ggAACGCGCGGCGAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATtg                                                                                                                                   <  1:753149/148‑1 (MQ=255)                 gAACGCGCGGCGAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATtgctgc                                                                                                                               <  2:46800/151‑1 (MQ=255)                        cggcgAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGAAGGCTGACTGGTTATTTCTGCTGATggg                                                                                                                        >  2:874705/1‑151 (MQ=255)                           cgAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGcc                                                                                                                     <  2:878962/151‑1 (MQ=255)                           cgAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGGCTGACTGTTTATTGCTGCTGATgggg                                                                                                                       >  2:951579/1‑149 (MQ=255)                            gAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCg                                                                                                                    <  2:641664/151‑1 (MQ=255)                            gAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGGGGACAGGCGGTTCATTGCTGCtggtgg                                                                                                                         >  1:656954/1‑146 (MQ=255)                            gAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGAGGACTGGCCGTTTCTTTGTGCTGGTgggggc                                                                                                                     >  1:880739/1‑148 (MQ=255)                             aaGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGc                                                                                                                   <  2:252040/151‑1 (MQ=255)                               gCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCa                                                                                                                  <  1:808591/150‑1 (MQ=255)                                 gTGCGGGATTCAGGGACGTGGTATTCCAGTGCGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCa                                                                                                                  <  1:574929/148‑1 (MQ=255)                                    ccGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCg                                                                                                            <  2:658935/151‑1 (MQ=255)                                      ggATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGtt                                                                                                          <  2:288071/151‑1 (MQ=255)                                      ggATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGCGCTGGTGGTGGATATCGGCGGTGGGATGGCCGGCATTTTATTGCTGCTCGTgggggggccggggg                                                                                                             >  1:875106/1‑138 (MQ=255)                                       gATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTc                                                                                                         <  2:921534/151‑1 (MQ=255)                                        aTTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGGCTGTTTATTGCTGCTGATGGGGGCGCAggggg                                                                                                             >  2:129012/1‑145 (MQ=255)                                        aTTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGGGGGACGGCCGACTGGTCATTTCTGCTGGTGGGGGCGCAgggggggt                                                                                                          >  1:772320/1‑145 (MQ=255)                                           cAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGAAGACTGACTGTTTATTGGTGCTGGTGGGGgcgccgggggggt                                                                                                          >  2:780855/1‑137 (MQ=255)                                            aGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTTATTGCTGCTGGTGGGGgcgccgggggggt                                                                                                          >  2:656356/1‑136 (MQ=255)                                             gggACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGtct                                                                                                   <  1:942208/151‑1 (MQ=255)                                                aCGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGGCTGTTTATTGCTGCTGGTGGGGGCGCAGTgggggt                                                                                                          >  2:659154/1‑137 (MQ=255)                                                   tGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGAGGACTGACTGTTTATTGCTGCTGATGGGGGCGGCGGGGCGTTCGCGTCTCGCTc                                                                                              >  1:20316/1‑150 (MQ=255)                                                   tGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTTACTGTTTATTTCTGCCGGTGGGGgcgccgggggggt                                                                                                          >  2:2779/1‑129 (MQ=255)                                                        ttCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGGAGACTGACTGTTCTTTGGTTCTGGTGGGGGCGCAGCGGGGGTCGCGTCTCGGTCg                                                                                             >  1:453070/1‑146 (MQ=255)                                                          ccAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAagcc                                                                                      <  2:578135/151‑1 (MQ=255)                                                               aCGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTTATTGCTTCTGGTGGGGGCGCCGGGGGGGTTGCGCCTCCATCGTGGagccagc                                                                                   >  1:813689/1‑149 (MQ=255)                                                                cGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCctg                                                                                <  2:115264/151‑1 (MQ=255)                                                                  aGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCctgct                                                                              <  1:905145/151‑1 (MQ=255)                                                                  aGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGGTCATTTCTGCTGGTGGGGGCGCAGTGGCGTTTGCGGCTCCGTCGTGGaaccagcc                                                                                  >  2:789765/1‑147 (MQ=255)                                                                              ctgggctggATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACAGAATGTTTCTTTCTGCTGGTGGGGGCGGAGGGGGGGTCGGGGTTCGCTCGGGGAGCCAGCCCGCTGGGGCCCAGTgg                                                                  >  2:58429/1‑151 (MQ=255)                                                                               tgggctggATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTCATTGCTGCTGGTGGGGGCGCAGTGGGGGTCGCGTCTCGATCGTGAAGCCAGCCCGCTGGGGCACAGTgg                                                                  >  2:522420/1‑150 (MQ=255)                                                                                gggctggATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTCATTGCTGCTGATGGGGGCGCCGGGGGGGTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGGCACAGTgggt                                                                >  2:541183/1‑149 (MQ=255)                                                                                    tggATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTTATTTCTGGTGGTGGGGGGGCCGTGGGGGTCGCGTCTCGATCGTGAAGCCAGCCCGCTTGGGGACAggggg                                                                 >  1:562744/1‑145 (MQ=255)                                                                                     ggATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCg                                                            <  1:498645/150‑1 (MQ=255)                                                                                      gATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGGCTGACTGTTTTTTTCTGCTGGTGGGGGCGCCGGGGGGTTTGCGGCCCGGTCTCGGAGGCAGCCCGCTGGGGCCCCGGGGTTg                                                               >  1:376282/1‑146 (MQ=255)                                                                                      gATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTa                                                          <  1:15105/151‑1 (MQ=255)                                                                                         tACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTTATTGCTGCTGGTGGGGCCGCCGGGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGGCACAGGGGTTGCCGTAtt                                                        >  2:632501/1‑150 (MQ=255)                                                                                           cGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTCATTGCTGCTGGTGGGGGCGCAGGGGGGTTTGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGGGGGTGCCGTAtt                                                        >  2:403432/1‑148 (MQ=255)                                                                                            gAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCg                                                    <  2:685545/151‑1 (MQ=255)                                                                                               gCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTCATTTCTGCTGGTGGGGGCGCCGTGGGGGTCGCGGCTCCATCCTGGAGCCAGCCTGCTGGGGCACAGGGGTTGGCGCTTTGgggg                                                   >  1:636774/1‑149 (MQ=255)                                                                                                ccACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACCGACTGTTCATTGCTTCTGGTGGGGGCGCAGGGGGGTTCGCGGCTCCATCGGGGAGCCAGCCCGCTGGGGCACAGGGGGTGGCGGATTGgggg                                                   >  1:768306/1‑148 (MQ=255)                                                                                                 cACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAAc                                               <  2:633666/151‑1 (MQ=255)                                                                                                 cACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTCATTGCTGCTGATGGGGGCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCCGCTGGGGCACAGTGGGTGCCGTAtttgggggaa                                                >  2:245906/1‑142 (MQ=255)                                                                                                   ccTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGa                                             <  1:541125/151‑1 (MQ=255)                                                                                                   ccTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTCATTGCTGCTGATGGGGGCGCAGTGGGGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACCGTGGGTGCCGTATTGGGGGGaa                                                >  2:178912/1‑148 (MQ=255)                                                                                                    cTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGAGGACTGAAAGGTCATTGGTGGTGATGGGGGCGCAGTGGCGGTCGCGCCTCGATCGTGAAGCCAGCCTGCTGGGGGACAGCGGGTGCCGGGTTGGCggg                                                  >  1:334569/1‑144 (MQ=255)                                                                                                     ttGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTCATTTCTGCTGATGGGGGCGCAGGGGGGTTCGCGTCTCGATCGTGGAGCCAGCCCGCTGGGTCACAGTGGGTGCCGCATTGGCggggaacgac                                           >  2:766645/1‑143 (MQ=255)                                                                                                       gCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTg                                         <  2:880842/151‑1 (MQ=255)                                                                                                             agaaAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATc                                   <  2:636845/151‑1 (MQ=255)                                                                                                             agaaAAACGGGTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTCATTGCTGCTGATGGGGGCGCCGGGGCGTTCGCGTCTCGATCGTGAAACCAGCCTGCTGGGGCACAGGGGGTGCCGTATTGgggggtaagaacttg                                        >  2:851551/1‑135 (MQ=255)                                                                                                                aaaaCGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATcgcg                                <  2:847185/151‑1 (MQ=255)                                                                                                                     gggTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGc                           <  1:764786/151‑1 (MQ=255)                                                                                                                     gggTGCTGGTGGTGGATATCGGCGGTGGGACGACTGACTGTTTATTGGTGCTGATGGGGGCGCCGTGGCGTTCGCGTCTCGCTCGTGGAGCCAGCCCGCTGGGGCACAGTGGGTGCCGTATTTGCGGGAAAGaatttggtattgcg                                >  2:145740/1‑133 (MQ=255)                                                                                                                      ggTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCg                          <  1:452115/151‑1 (MQ=255)                                                                                                                        tGCTGGTGGTGGATATCGGCGGTGGGACGGCTGACTGTTCATTGCTGCTGATGGGGGCGCCGGGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGGCACAGTGGGTGCCGGATTGGCGGGGAAGATCTggg                                       >  2:588576/1‑135 (MQ=255)                                                                                                                         gCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGttt                       <  2:908517/151‑1 (MQ=255)                                                                                                                            ggtggtggATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGTTTaaa                    <  1:713124/151‑1 (MQ=255)                                                                                                                               ggtggATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGTTTAAAAAc                 <  1:753116/151‑1 (MQ=255)                                                                                                                                      atCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTGACGATCTGGATATCGCGCTGGCGTTTAAAAACCTGATg           <  1:758684/150‑1 (MQ=255)                                                                                                                                      atCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGCGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGTTTAAAAACCTGATg           <  2:770264/150‑1 (MQ=255)                                                                                                                                          gcggTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGTTTAAAAACCTGATGCCAct      >  1:346316/1‑151 (MQ=255)                                                                                                                                          gcggTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGTTTAAAAACCTGATGCCAc       >  1:623850/1‑150 (MQ=255)                                                                                                                                             gtggtACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGTTTAAAAACCTGATGCCActgct   <  1:941001/151‑1 (MQ=255)                                                                                                                                             gtggtACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGTTTAAAAACCTGATGCCAc       >  2:626676/1‑147 (MQ=255)                                                                                                                                             gtggtACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGGAATGATCTGGATATCGCGCTGGCGGTTTAAAATCTGATGCCActgct   >  1:797209/1‑151 (MQ=255)                                                                                                                                               ggtACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGTTTAAAAACCTGATGCCActgctg  <  2:871919/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GGCGCAAGGGATTCTGGAACGCGCGGCGAAGCGTGCCGGATTCAGGGACGTGGTATTCCAGTACGAGCCGGTCGCGGCTGGGCTGGATTACGAAGCCACCTTGCAGGAAGAAAAACGGGTGCTGGTGGTGGATATCGGCGGTGGTACGACTGACTGTTCATTGCTGCTGATGGGGCCGCAGTGGCGTTCGCGTCTCGATCGTGAAGCCAGCCTGCTGGGTCACAGTGGTTGCCGTATTGGCGGTAACGATCTGGATATCGCGCTGGCGTTTAAAAACCTGATGCCACTGCTG  >  NC_000913/2148186‑2148477

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