Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,941	0	A	C	43.3%	103.0 / 16.5	67	E392D (GAA→GAC)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (30/8);  new base C (1/28);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.23e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATA  >  NC_000913/3026794‑3027087                                                                                                                                                       |                                                                                                                                                   cTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTGT‑ACC‑TGACAAAGAGGCTGATTTCGTGGTGATGGAAc                                                                                                                                                   <  2:627211/149‑1 (MQ=255)  tACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAAccc                                                                                                                                                 >  1:448267/1‑150 (MQ=255)  tACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAAcc                                                                                                                                                  >  2:742510/1‑149 (MQ=255)      aaGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGc                                                                                                                                            <  2:871294/151‑1 (MQ=255)       aGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGcc                                                                                                                                           <  1:304735/151‑1 (MQ=255)             gCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCg                                                                                                                                     >  1:906102/1‑151 (MQ=255)                aTTACAGGGCTACCGCCTCTCGGCATAGGAAGGGTTTTACCTGGCCACGCTCGGCGGAGCAAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑CCC‑GGGAAAAGAGGCTGATTTCGGGGTGAGGGACCCCACCGCCACTCCGCTa                                                                                                                                  <  2:1041513/151‑1 (MQ=255)                 ttACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAc                                                                                                                                 >  2:1309933/1‑151 (MQ=255)                    cAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAcagc                                                                                                                              >  1:443509/1‑151 (MQ=255)                           aTCGCCTCCCGGCATATGAAGCGTTTTCCCGGGCCCCGCTGGGCGGAGCGAAACCTCGGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑CCC‑GGGAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCt                                                                                                                          <  1:553828/148‑1 (MQ=255)                            tCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTgcgc                                                                                                                      >  1:1165720/1‑151 (MQ=255)                            tCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTgcgc                                                                                                                      >  2:424475/1‑151 (MQ=255)                             ccccTCTGGGCATATGAAGGTTTTTACCTGGCCACGCTCGGCGGGGCAAAACCTCTGGCCCTTGACGATTTG‑TTGGCAACTTTTTT‑CCC‑TGGAAAAGAGGTTGATTTCGTGGTGATGAACCCCGCCGCCACTCCGCTACAGCAGCTgc                                                                                                                        <  1:378133/146‑1 (MQ=255)                             cGCCTCCCGGAATATGAGGTTTTTTCCCGGGCCCCCCGCGGGGGAGGGAAATCTCGGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑CCC‑GGCAAAAGAGGCTGATTTGGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTgcg                                                                                                                       <  2:425370/149‑1 (MQ=255)                                  cgcGGCATATGAAGCTTTTTACCGGGCCGCGCTCGGCGGAGCGAAATCTCTGGGCCTTAACGATTTGATTGGCAAC‑TTTTT‑GCC‑TGGCAAAGAGGTTGATTTCGTGGTGGTGGACCCAACCGCCACTCCGCTACAGCAGCTGCGCTATg                                                                                                                  <  2:1327053/147‑1 (MQ=255)                                     ggCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAc                                                                                                             >  1:472237/1‑151 (MQ=255)                                        atatGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActct                                                                                                          >  2:295425/1‑151 (MQ=255)                                         tatGAAGCGTTTTACCTGGCCACGCTCGGCGGCGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCGGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActct                                                                                                          >  1:1261112/1‑150 (MQ=255)                                             aaGCGTTTTACCTGGCCCCGCTCGGCGGAGAAAAATCTCGGGGCCTTGACGATTTGATTGGAAAT‑TTTTT‑CCC‑GGGAAAAGAGGCTGATTTGGTGGTGAGGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTg                                                                                                         <  1:556816/147‑1 (MQ=255)                                             aaGCGTTTTACCGGGCCACGCGCGGCGGAGCGAAATCTCTGGGCTTTGACGATTTGATTGGCACC‑TTTTT‑ACC‑TGGAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGtt                                                                                                       <  1:280012/149‑1 (MQ=255)                                              aGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCt                                                                                                    >  2:965956/1‑151 (MQ=255)                                              aGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCt                                                                                                    >  2:342290/1‑151 (MQ=255)                                                    ttACCTGGCCCCGCTGGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCACC‑TTTTT‑CCC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt                                                                                               <  1:209734/150‑1 (MQ=255)                                                    ttACCTCGCCAGGCTCGGCGGAGCGAAATCTCTGGGCTTTGACGATTTGATTCGCACC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGCGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt                                                                                               >  2:909681/1‑150 (MQ=255)                                                       ccgggCCACGCTCGGCGGAGCAAAATCTCTGGGCCTTGACGTTTTGATTGGCAAT‑TTTTT‑CCC‑GGGAAAAGGGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGAc                                                                                           <  1:477971/148‑1 (MQ=255)                                                         tGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACa                                                                                          >  2:378825/1‑150 (MQ=255)                                                          ggCCCCGCTGGGCGGAGCAAAATCTCGGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑CCC‑TGGAAAAGAGGCTGATTTCGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACaaa                                                                                        <  1:708402/151‑1 (MQ=255)                                                             cccccgcggcggAGCGAAATCTCGGCGCTTTGACGTTTTGATTGGCAAC‑TTTTT‑CCC‑TGACAAAGAGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACaaa                                                                                        <  2:217428/142‑1 (MQ=255)                                                              aCGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGt                                                                                    >  2:206087/1‑151 (MQ=255)                                                              aCGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGt                                                                                    >  2:460114/1‑151 (MQ=255)                                                               cGCTCGGCGGAGCGAAATCTCTGGGCTTTGACGATTTGATTGGCAAT‑TTTTT‑CCC‑TGGCAAAGAGGCTGTTTTCGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGtt                                                                                   <  1:203358/151‑1 (MQ=255)                                                                 cgcggcggAGCAAAATCTCTGGGCTTTGACGTTTTGATTGCCAAT‑TTTTT‑CCC‑GGGAAAAGGGGCTGATTTCGGGGTGATGAAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCg                                                                                 <  2:1179765/149‑1 (MQ=255)                                                                  tCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCg                                                                                 >  1:336416/1‑150 (MQ=255)                                                                  tCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGt                                                                                >  1:297839/1‑151 (MQ=255)                                                                   ggggggcgggAAATCATTGGGCCTTAACG‑TTTAAATGGAAAA‑TTTTTCGCC‑GGGCAAAAAGCCTTTTTGCGGGGTGTGAGACCCGCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc                                                                                  <  1:700131/139‑1 (MQ=255)                                                                    ggcgggggaaaaTCTCTGGGCTTTGCCGTTTTGATTGCCATC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTGGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCg                                                                                 <  2:24532/139‑1 (MQ=255)                                                                    ggcggAGCGAAATCTCGGGCCCTTGACGTTTTGATTGCCACC‑TTTTT‑ACC‑TGGAAAAGGGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCg                                                                                 <  2:336410/148‑1 (MQ=255)                                                                    ggcggAGAGAAATCTCGGGGCCTTGACGTTTTGGTTGGCAAC‑TTTTT‑CCC‑TGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtg                                                                               <  1:639729/150‑1 (MQ=255)                                                                     ggggAGCGAATTCTCGGGGCTTTGACGATTTGATTGGCACC‑TTTTT‑CCC‑TGCCAAAGAGTCTAATTTGGTGGTGTTGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc                                                                                  <  1:253881/144‑1 (MQ=255)                                                                      cggAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatg                                                                            >  1:715419/1‑151 (MQ=255)                                                                      cggAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatg                                                                            >  1:195784/1‑151 (MQ=255)                                                                        gAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgat                                                                          >  1:1341601/1‑151 (MQ=255)                                                                             aaaTCTGTGGCCCTTGGCGTTTTGATTGGCACT‑TTTTC‑CCC‑GGGAAAAGGGGTTTTTTGCGTGGTAAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatga                                                                        <  2:698393/148‑1 (MQ=255)                                                                                    tGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCg                                                               >  1:143167/1‑150 (MQ=255)                                                                                     gggCCTTGACGATTTGATTGGCAAT‑TTTTT‑ACC‑TGGAAAAGAGGCTGATTTGGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTgg                                                                  <  2:448259/146‑1 (MQ=255)                                                                                     gccccTTGACTATTTAATGGGCACT‑TTTTT‑CCT‑TGGAAAAGGGTCTGTTTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGAt                                                             <  2:629370/148‑1 (MQ=255)                                                                                        ccTGGACGTTTTGATTGGCAAC‑TTTTT‑ACC‑GGGAAAAGAGGTTTTTTTCGTGGTGATGGACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATg                                                            <  2:660794/149‑1 (MQ=255)                                                                                          ttGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCg                                                        >  2:724655/1‑151 (MQ=255)                                                                                            gACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGtt                                                      >  2:855322/1‑151 (MQ=255)                                                                                              cGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCg                                                    >  2:436631/1‑151 (MQ=255)                                                                                                aTTTGTTTGGCAAT‑TTTTT‑ACC‑GGGAAAAGGGGGTGATTGCGGGGTGATGGAACCAACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGa                                                   <  1:1007239/150‑1 (MQ=255)                                                                                                 tttGATTGGCATC‑TTTTT‑CCC‑GGAAAAAGAGGCTGATTTGGTGGAGGTGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGt                                                       <  2:500490/145‑1 (MQ=255)                                                                                                    ttgatgcgcagttttt‑ACCATGGAAAAGAGGTTGATTTCGTTGTAAGGCACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGAt                                                  <  1:1008433/138‑1 (MQ=255)                                                                                                        ggCAAC‑TTTTT‑CCC‑GGGCAAAGAGGCTGATTTCGGGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCg                                            <  2:472229/149‑1 (MQ=255)                                                                                                         ggcattttttc‑ACC‑TGAAAAAGAGGCTGTTTGCGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCg                                            <  2:462066/143‑1 (MQ=255)                                                                                                            aatttttc‑ACG‑TGGAAAAGGGGGTTATTCCGTGGTGATGCACCCCGCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTAcc                                             <  1:339216/142‑1 (MQ=255)                                                                                                                         tGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATgg                             >  2:338605/1‑150 (MQ=255)                                                                                                                             aaaGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgtc                        <  1:424483/151‑1 (MQ=255)                                                                                                                              aaGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCt                       >  1:686874/1‑151 (MQ=255)                                                                                                                                  ggttgtttgcGTGGTGTGGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgt                    <  2:1016373/141‑1 (MQ=255)                                                                                                                                  ggCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtg                   <  1:537212/151‑1 (MQ=255)                                                                                                                                   gCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtgt                  >  2:431109/1‑151 (MQ=255)                                                                                                                                        tttCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGa              <  1:561164/150‑1 (MQ=255)                                                                                                                                        tttCGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTa                 <  1:965974/147‑1 (MQ=255)                                                                                                                                            gtggtgATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGc          >  1:644111/1‑150 (MQ=255)                                                                                                                                               gtgATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACt      >  2:26779/1‑151 (MQ=255)                                                                                                                                                    ggAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaata  >  2:525209/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAAC‑TTTTT‑ACC‑TGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATA  >  NC_000913/3026794‑3027087

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