Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,999,673	1	.	G	57.8%	130.6 / 110.7	64	coding (219/756 nt)	ygeR	LysM domain‑containing M23 family putative peptidase; septation lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (16/11);  new base G (22/15);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

TGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGC  >  NC_000913/2999524‑2999818                                                                                                                                                       |                                                                                                                                                  tGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTaa                                                                                                                                                    <  1:347176/150‑1 (MQ=255)    ccTACTCGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGTCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGCCCAACTTCAATGGTGTAA‑gggg                                                                                                                                               <  1:498339/151‑1 (MQ=255)     cTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑ggggg                                                                                                                                              >  2:1135212/1‑151 (MQ=255)     cTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑ggggg                                                                                                                                              >  2:30150/1‑151 (MQ=255)          ggcggcCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCGCCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGa                                                                                                                                           <  1:465748/149‑1 (MQ=255)          ggcggcCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑gggggggga                                                                                                                                          <  1:1047868/150‑3 (MQ=255)                cAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGt                                                                                                                                   >  1:401961/1‑151 (MQ=255)                         ttCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTcg                                                                                                                            <  1:680918/150‑1 (MQ=255)                                ccGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGtt                                                                                                                    <  2:629175/151‑1 (MQ=255)                                   cTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGGCGCGCCAGTTCtt                                                                                                                 >  1:816793/1‑151 (MQ=255)                                     gatgaCGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCtttt                                                                                                               <  1:114450/150‑1 (MQ=255)                                       tgaCGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTAc                                                                                                             >  1:712210/1‑151 (MQ=255)                                        gaCGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTAct                                                                                                            <  2:132303/151‑2 (MQ=255)                                        gaCGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACg                                                                                                            <  1:749372/151‑1 (MQ=255)                                         aCGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGc                                                                                                           >  1:443891/1‑151 (MQ=255)                                          cGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGc                                                                                                           >  2:898514/1‑149 (MQ=255)                                               tAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTTCGGGTAATTTTGCTACTTTTCGCCCCACCCAGTTTTAGCTTCTGTTCAACTTTCATGGTGGAA‑GTGGTGGAAATTTCCTTCACTCGCCCCAGTTCTTTTACGCTGGTTcc                                                                                                    >  1:1111673/1‑151 (MQ=255)                                                aaCCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTcc                                                                                                    >  1:993455/1‑151 (MQ=255)                                                aaCCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTcc                                                                                                    >  2:941509/1‑151 (MQ=255)                                                 aCCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTccc                                                                                                   >  1:306432/1‑151 (MQ=255)                                                 aCCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCg                                                                                                  <  1:811041/151‑1 (MQ=255)                                                   cGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCCGTCGCGCCAGTTCTTTTACGCTGGTTccc                                                                                                   >  1:321118/1‑149 (MQ=255)                                                    gATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtg                                                                                                <  2:389102/151‑1 (MQ=255)                                                    gATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtg                                                                                                >  2:1219393/1‑151 (MQ=255)                                                    gATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtgg                                                                                               >  2:637083/1‑151 (MQ=255)                                                      tGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtgg                                                                                               >  1:844519/1‑150 (MQ=255)                                                         ggTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtggtggg                                                                                           >  2:41126/1‑148 (MQ=255)                                                                gTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATa                                                                                    <  1:200104/151‑1 (MQ=255)                                                                     tGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGata                                                                               >  2:606213/1‑150 (MQ=255)                                                                      gATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGata                                                                               >  1:1228860/1‑150 (MQ=255)                                                                          tGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGc                                                                           <  1:204960/150‑1 (MQ=255)                                                                           ggCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCg                                                                          >  1:1274753/1‑150 (MQ=255)                                                                                     ttACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGtt                                                              >  1:538041/1‑151 (MQ=255)                                                                                      tACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGttt                                                             >  2:809665/1‑151 (MQ=255)                                                                                       aCGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGttt                                                             <  2:1082143/151‑1 (MQ=255)                                                                                       aCGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTc                                                            >  1:102001/1‑151 (MQ=255)                                                                                               tACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgtgt                                                     >  1:664251/1‑151 (MQ=255)                                                                                                aCTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGgtgt                                                     >  2:209189/1‑150 (MQ=255)                                                                                                   gCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAAc                                                 >  2:815992/1‑151 (MQ=255)                                                                                                   gCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACg                                                <  1:1142215/151‑1 (MQ=255)                                                                                                      aCTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACgg                                               >  2:1114991/1‑150 (MQ=255)                                                                                                        ttttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGcc                                           >  2:247275/1‑151 (MQ=255)                                                                                                         tttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCg                                          >  2:726600/1‑151 (MQ=255)                                                                                                          tttggcccacccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAg                                             <  1:782236/142‑1 (MQ=255)                                                                                                          ttCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCg                                          >  2:876331/1‑151 (MQ=255)                                                                                                               cccacccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATAc                                    <  1:1201045/151‑1 (MQ=255)                                                                                                                ccacccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATAc                                    >  2:338175/1‑151 (MQ=255)                                                                                                                  acccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGtt                                 >  1:166426/1‑151 (MQ=255)                                                                                                                   cccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGtt                                 <  2:625908/151‑1 (MQ=255)                                                                                                                     caGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCt                              >  2:518086/1‑151 (MQ=255)                                                                                                                      aGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCt                              >  2:653606/1‑151 (MQ=255)                                                                                                                      aGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCt                              <  2:844556/151‑1 (MQ=255)                                                                                                                            aGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATc                         <  1:1064392/150‑1 (MQ=255)                                                                                                                             gTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCgg                       >  1:862893/1‑150 (MQ=255)                                                                                                                               ttCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGgatg                    >  2:655536/1‑151 (MQ=255)                                                                                                                                gatcggcaaagcgaaTGGTGTGG‑GAAGGGGGAAGGCCGTTGAGTCGAGCCAGGTGGTTGACGCTTGTTTCCGAGATGCGCGAAATACGATATAGAGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGATCCTGTATCGga                      >  1:58729/14‑148 (MQ=255)                                                                                                                                   gACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGAttt                <  1:834345/151‑1 (MQ=255)                                                                                                                                   gACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGAtt                 <  2:1182651/150‑1 (MQ=255)                                                                                                                                     ccAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTg               >  2:1063121/1‑151 (MQ=255)                                                                                                                                      cAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGgat                     <  1:97230/144‑1 (MQ=255)                                                                                                                                       aaCTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTg               <  1:642794/148‑1 (MQ=255)                                                                                                                                          ttCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTAcc          >  2:497209/1‑151 (MQ=255)                                                                                                                                            cAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGa       <  2:522021/151‑1 (MQ=255)                                                                                                                                             aaTGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCg        <  1:30158/150‑1 (MQ=255)                                                                                                                                                ggTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAAc     <  2:175862/150‑1 (MQ=255)                                                                                                                                                ggTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACa    <  1:747579/151‑1 (MQ=255)                                                                                                                                                ggTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAg   <  2:286845/151‑1 (MQ=255)                                                                                                                                                  tgtAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTGACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGc  >  1:544419/1‑151 (MQ=255)                                                                                                                                                  tgtAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGc  >  2:1170040/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                  TGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGC  >  NC_000913/2999524‑2999818

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