Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,547,772	0	T	G	43.8%	119.0 / 15.0	73	V367G (GTA→GGA)	murP	N‑acetylmuramic acid permease, EIIBC component, PTS system
Reads supporting (aligned to +/- strand):  ref base T (18/23);  new base G (32/0);  total (50/23)
Fisher's exact test for biased strand distribution p-value = 5.46e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTG  >  NC_000913/2547629‑2547920                                                                                                                                                |                                                                                                                                                      tATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCgg                                                                                                                                                 <  1:621606/150‑1 (MQ=255)  aTTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCCCACAGGGCGCTACGCAGGCAGGGACGCgggg                                                                                                                                               >  2:177895/1‑151 (MQ=255)  aTTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCCCCCAATGCGCCCCGCCGGCCGGGACGCgggg                                                                                                                                               >  1:1052126/1‑151 (MQ=255)    tCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCCCCCAGTGCGCTACGCCGGCAGGGGCGCgggggg                                                                                                                                             >  1:1080769/1‑149 (MQ=255)    tCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGc                                                                                                                                              <  1:1129064/150‑1 (MQ=255)       tGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCCCAGAGCGCCACGCCGGCAGGGACGCgggggggtt                                                                                                                                          >  1:784399/1‑146 (MQ=255)             ccTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCCCACAATGCGCTACGCAGTCAGGGACGCGGGGCGG‑TTATTccc                                                                                                                                   >  2:518176/1‑151 (MQ=255)                cgcgTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCg                                                                                                                                  <  1:1044531/149‑1 (MQ=255)                 gcgTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGGGGGACCGCCCAGGGGGCTACGCCGCGCGGGAAGCGGGGGGA‑TTTTTCCCgggc                                                                                                                               >  1:372621/1‑149 (MQ=255)                  cgTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCGCACAGTGCGCCACGCAGGCCGGGACGCGGGGGGG‑TTATTTCCCGcc                                                                                                                               >  1:64866/1‑150 (MQ=255)                       aTGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTgg                                                                                                                         <  1:126538/151‑1 (MQ=255)                           aCAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGt                                                                                                                     <  2:796673/151‑1 (MQ=255)                                     aTTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCg                                                                                                           <  2:1135523/151‑1 (MQ=255)                                     aTTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCCCCCACGGCGGCACGCCGGCAGGGACGCGGGGGGG‑GTATTTCCCGCCTGCTGGGCGGTgg                                                                                                                  >  1:823490/1‑144 (MQ=255)                                      ttCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCCCACAGCGCGCTACGCAGTCAGGGACGCGGGGCGA‑TTATTCCCCGCCTGCTGGGCGTTGGTGGAcc                                                                                                            >  2:961948/1‑149 (MQ=255)                                      ttCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAAACCCACCGGGCGCTACGCTCTCTGGGAATCCGGGCGG‑TTTTTTCCCGCCCGCTCGGGGCTTGTGAAcc                                                                                                            >  2:1079846/1‑149 (MQ=255)                                         aaCAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCCCCCACTGCGCCACGCCGGCCGGGACGCGGGGCGG‑TTATTTCCCGCCCGCTggggg                                                                                                                      >  1:720271/1‑134 (MQ=255)                                         aaCAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGGGCGCTACGCAGGCAGGGACGGGGGGGGG‑TTTTTTCCGGCCTGCGGGGGGGTGGGGGACCGCTg                                                                                                        >  2:452165/1‑150 (MQ=255)                                         aaCAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACACTGCGCTACGCAGTCAGGGACGCGGGGGGA‑TTATTTCCGGCCTGCTGGGCGGTGGTGAAACGCTg                                                                                                        >  2:698282/1‑150 (MQ=255)                                             gCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAAACCCCCAGTGCGCCACGCCGTCAGGGACGCGGGGCGG‑TTTTTCCCCGCCTGCCGGGCGGTGGGGGACCGCTGGtttt                                                                                                   >  1:661289/1‑150 (MQ=255)                                              cTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGGCAGGTACGCGGGGCGA‑TTATTTCCCGCCTGCTGGGCGTTGGTGGACCGCTGAttt                                                                                                    >  2:863471/1‑149 (MQ=255)                                              cTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCCCACAGTGCGCTACGCCGTCAGGGACGAGGGGCGG‑TTATTTCCCGCCTGCTGGGCGGTGGTGAACCGCTGAttt                                                                                                    >  2:914605/1‑149 (MQ=255)                                              cTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAAACGCACAGGGCGCTACGCAGTCAGGGACGCGGGGCGA‑TTATTCCCCGCCTGCTGGGCGGTGGGGAACCGCTGAttt                                                                                                    >  2:1114124/1‑149 (MQ=255)                                              cTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGAGCAAACGCACAGTGCGCCACGCAGGCAGGGACGCGGGGGGG‑TTATTTCCGGCCCGCTGGGCGGTGGTGAACCGCTGGttt                                                                                                    >  2:828335/1‑149 (MQ=255)                                                tATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCAACAGACAGGGTGCTAAGCAGACAGGGACGCCGGGGGA‑TTATTTCCCGCCCGCTGGGCGGTTGTGAAACCCTGAttt                                                                                                    >  2:456715/1‑147 (MQ=255)                                                     cccATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGGGCGACCGCACAGTGCGCCACGCAGTCAGGGACGCGGGGGGA‑ATTTTTCCCGCCCGCTGGGCGTTGGGGGACCGCCGATTTACGgtg                                                                                              >  2:1158354/1‑148 (MQ=255)                                                      ccATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCGCCCACTGCGCTACGCAGGCAGGGACGCGGGGGGG‑TTATTTCCCGCCCGCCCGGGGTTGGGGGAACGCTGATTTTCGGGGTccc                                                                                          >  1:50454/1‑149 (MQ=255)                                                      ccATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCAc                                                                                          <  2:935118/151‑1 (MQ=255)                                                         tCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGGGCGCTACGCAGTCCGGGACGCGGGGGGG‑TTATTCCCCGCCCGCTGGGCGGTGGGGAACCGCTGGTTTTCGGTGTCAccc                                                                                        >  2:212319/1‑150 (MQ=255)                                                          ccTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGGGCAACCGCAACGTGCGCCACGCCGTCAGGTCCGGGGGGGGG‑TTTTTTCCGGCCCGCCGGGGGGTGGGGGACCGCTGGTTTTCGGGGTCAccc                                                                                        >  1:490346/1‑149 (MQ=255)                                                             ttCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCGCACACGGCGCTACGCAGTCAGGGACGCGGGGGGG‑TTATTTCCGGCCCGCTGGGCGGTGGGGGAACGCTGGTTTTCGGTGTCCCCCCGccc                                                                                   >  1:101213/1‑151 (MQ=255)                                                             ttCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGc                                                                                     <  2:661299/149‑1 (MQ=255)                                                               cAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGcccc                                                                                  <  1:919575/150‑1 (MQ=255)                                                                   ggCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATg                                                                             <  1:456712/151‑1 (MQ=255)                                                                    gCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGa                                                                            <  1:477337/151‑1 (MQ=255)                                                                     cAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCCCAGGGCGCTACGCAGGCAGGTAAACGGGGGGG‑TTATTCCCCGCCCGCTGGGCGTTGGGGAAACGCTGATTTACGGTGTCACCCTGCCCCGCCTGa                                                                            >  1:1064618/1‑150 (MQ=255)                                                                      aGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGGGCCACCCCACAGGGCGCTACGCCGGCAGGGACGCGGGGGGG‑GTATTTCCCGCCTGCTGGGCGGTGGGGGACCGCCGGTTTTTGGGGGCACCCCGCCCCCCCTGaa                                                                           >  1:900994/1‑150 (MQ=255)                                                                      aGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAg                                                                          <  2:815769/151‑1 (MQ=255)                                                                       ggCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAg                                                                          <  1:267935/150‑1 (MQ=255)                                                                         cgcgGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGGACGCGGGGCGA‑TTATTTCCCGCCCGCTGGGCGTTGGGGAACCGCTGATTTACGGTGTCACCCTGCCCCCCATGAAGcc                                                                        >  2:681894/1‑150 (MQ=255)                                                                          gcgGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCCACCCCACAGGGGGCTACGCAATGAGGGGCGCCGGGCGACTTTTTCCCGGCCT‑CTGGGGGGTGGTGCAGCGCCGATTTACGGGGTTCACCTCCCCCGccata                                                                            >  1:113122/1‑141 (MQ=255)                                                                           cgGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGtt                                                                     <  2:930079/151‑1 (MQ=255)                                                                            gggCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGGGCGCTACGCAGTCAGGGACGCGGGGCGG‑TTATTGCCGGCCTGCTGGGCGTTGGTGAACCGCCGATTTACGGGGTCACCCTGCCCCGCATGAAACCGttt                                                                    >  2:640205/1‑151 (MQ=255)                                                                               ccAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTTCCCGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGGGTCACCCTGCCCCGCATGAAGCCGTTTGtt                                                                 >  2:132887/1‑151 (MQ=255)                                                                               ccAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGtt                                                                 <  2:756127/151‑1 (MQ=255)                                                                                cAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTa                                                                <  2:297395/151‑1 (MQ=255)                                                                                 aGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGAGCGCTACGCAGGCAGGGGCGCGGGGGGG‑ATATTTCCGGCCCGCGGGGGGTTGGTGCAACGCCGATTTTCGGGGTCACCCCGCCCCGCATGAAGCCGTTTGtt                                                                 >  1:751308/1‑149 (MQ=255)                                                                                    tGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTTCCGGCCTGCTGGGCGGTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCCTGAAGCCGTTTGTTTcc                                                              >  2:622598/1‑149 (MQ=255)                                                                                     gggCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGc                                                            <  2:968319/150‑1 (MQ=255)                                                                                      ggcgcggcgcTGGCACTCTACTGGCGGGCGCAAACGCACAGGGCGCCACGCAGGCAGGGACGCGGGGCGG‑TTATTTCCGGCCTGCTGGGGGTTGGGGAACCGCTGATTTACGGGGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGcc                                                           >  2:455845/1‑150 (MQ=255)                                                                                           ggcgcTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCCACGCCGGCAGGGACGCGGGGGGG‑TTTTTTCCCGCCTGCTGGGGGTTGGGGGAACGCCGATTTACGGTGTCACCCCGCCCCCCCTGAAGCCGTTTGTTTCCGCCTGtttt                                                     >  1:968882/1‑150 (MQ=255)                                                                                             cgcTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAgg                                                   <  2:1148422/151‑1 (MQ=255)                                                                                              gcTGGCACTCTACTGGCGGGCGCCAACGGACAGTGCGCCACGCAGGCAGGGAAGCGGGGGGG‑GTATTCCCGGCCTGCTGGGGGGTGGGGAAACGCTGGTTTTCGGGGGCACCCTGCCCCGCATGAAGACGGTTTTTACCGCCTGGTTAggg                                                  >  1:328894/1‑150 (MQ=255)                                                                                                tGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCCGGGACGCGGGGCGG‑TTATTTCCGGCCTGCTGGGCGGTGGTGAACCGCCGGTTTTCGGGGTCACCCTCCCCCCCATGAAGCCGGTTGTTTCCGCCCGTTTTggggg                                                >  1:1083332/1‑151 (MQ=255)                                                                                                 ggCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAggcggc                                               <  2:540264/151‑1 (MQ=255)                                                                                                        tACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGGACGCGGGGGGG‑TTATTTCCCGCCTGCTGGGGGGTGGGGGACCGCTGATTTTCGGGGTCACCCCGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGGGGCGcggcgg                                        >  1:40006/1‑151 (MQ=255)                                                                                                         aCTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCggcggg                                       <  2:1185675/151‑1 (MQ=255)                                                                                                            ggcgggcgCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCggcgggc                                      <  1:1068452/149‑1 (MQ=255)                                                                                                            ggcgggcgCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGgcgggcgg                                    >  1:941242/1‑151 (MQ=255)                                                                                                             gcgggcgCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGgcgggcgg                                    >  1:977920/1‑150 (MQ=255)                                                                                                                ggcgCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGgttt                                 >  2:998101/1‑150 (MQ=255)                                                                                                                     aaCCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTAt                           >  1:1023325/1‑151 (MQ=255)                                                                                                                         gCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGc                       >  1:142644/1‑151 (MQ=255)                                                                                                                                     aCGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCt             <  2:311244/149‑1 (MQ=255)                                                                                                                                      cGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCtgg           >  1:734103/1‑150 (MQ=255)                                                                                                                                         aGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTggg       <  2:70849/151‑1 (MQ=255)                                                                                                                                         aGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTggg       >  2:817380/1‑151 (MQ=255)                                                                                                                                         aGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTggg       >  1:506892/1‑151 (MQ=255)                                                                                                                                         aGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTggg       >  1:356571/1‑151 (MQ=255)                                                                                                                                          gTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGgggggg      >  1:99639/1‑151 (MQ=255)                                                                                                                                            cAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTc    >  2:992917/1‑151 (MQ=255)                                                                                                                                             aGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTgggggc    >  1:617345/1‑148 (MQ=255)                                                                                                                                               gTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTg  >  1:398755/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                      TATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGA‑TTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTG  >  NC_000913/2547629‑2547920

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