Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,159	0	G	A	42.9%	73.3 / 41.3	70	intergenic (+139/‑75)	thrC/yaaX	L‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (31/9);  new base A (0/30);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 3.89e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCT  >  NC_000913/5012‑5302                                                                                                                                                        |                                                                                                                                                catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGaa                                                                                                                                               <  1:425965/150‑1 (MQ=255) catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAAt                                                                                                                                              <  1:727710/151‑1 (MQ=255) catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAAt                                                                                                                                              <  2:590362/151‑1 (MQ=255)    caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑TTG‑AAAAATTTTTGGGGAAAAATGCAGGGGAAAAGGGAGAAATTCTAATTAAATGGGGAAACTT‑AGGGATTAGAATGGGGGAGAATa                                                                                                                                             >  1:224649/1‑149 (MQ=255)    caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATaa                                                                                                                                            <  1:449802/150‑1 (MQ=255)      gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATaacaa                                                                                                                                         <  2:310374/151‑1 (MQ=255)      gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATaaca                                                                                                                                          >  2:274275/1‑150 (MQ=255)        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATGAGGATTGCGGAGAATaacaaca                                                                                                                                       <  2:1040141/151‑2 (MQ=255)          aaTCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                     >  1:833032/1‑151 (MQ=255)               aTTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                    <  2:1123044/147‑1 (MQ=255)                 tCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                              >  1:418776/1‑151 (MQ=255)                      aTCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                          >  2:658688/1‑150 (MQ=255)                          cAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTa                                                                                                                      >  1:889638/1‑150 (MQ=255)                           aaccAGGCGGGGTTTCT‑TTTTAGGCAGCCGGTCTTTTTA‑TTA‑AAGAAATTAGGGAAAAAAAGGACAGGAAAAAGGGAGAAATTCCAAATAAAAGGGGTAACTT‑AGAGATAAGGATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAAt                                                                                                                    <  2:960093/148‑1 (MQ=255)                           aaacAGGGCGGGTTTTCTTTTTATGGA‑CCCCGCTTTTTT‑ATGAAAAAAATAAGGGAAAAAAAGGACA‑GGAAAAAGGAAAAAATCCAAAAAAAAGGGGAAAATTAAGAATTTGGGTTGGGGGAAAATAACAACCGCCGTTCTCATCGAGTaa                                                                                                                     <  1:1023039/148‑1 (MQ=255)                               aGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                >  1:945793/1‑151 (MQ=255)                                ggCGGGTTTTGT‑TTTTATGCAGCCCGGTTTTTTT‑ATG‑AAAAAATTAGGGAAAAAAAGGACAGGAAAAAAGGAAAAATTCAAAAAAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAAAATACCAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                               <  2:27701/151‑1 (MQ=255)                                 ggcgggttttttttttATGCAGCCCGGTTTTTTT‑AAA‑AAAAAATAAGGAAAAAAAAGGACAGGAAAAAAGGAAAAATTTTAAAAAAATGCGGTAATTT‑AGAGATTGGGATGGGGAAAAATAACAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                 <  2:197664/137‑1 (MQ=255)                                   cGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                            >  1:873747/1‑151 (MQ=255)                                             ttttATGCA‑CCCGGTTTTTTTAAAG‑AAAAATTTAGGGAAAAAAAGGACAGGAAAAAAGGAAAAATTCAAAAAAAAGGCGTTACCTT‑AAAGATTAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                    <  2:670881/150‑1 (MQ=255)                                              tttATGCAGCCCGTTTTTTTT‑ATA‑AAAAAATTAGGGAAAAAAAGGACAGGAAAAAAGGAAAAATTCTCAATAAAGGCGGTAACTA‑AGAGATTAGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                   <  1:316017/150‑1 (MQ=255)                                              tttAGGCACCCCGGTTTTTTT‑AGA‑AAAAAATTAGGGAAAAAAAGGAAAGGAAAAAAGGAAAATTTCTAAAAAAAGGGGATAATTA‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                   <  1:477686/150‑1 (MQ=255)                                                 aTGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                               >  1:297457/1‑151 (MQ=255)                                                  tGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                              >  2:1101227/1‑151 (MQ=255)                                                       cccgtttttttt‑ATA‑AAAAAATTAGGGAAAAAAATGAAGGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAATTA‑AGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                           <  2:652301/143‑1 (MQ=255)                                                       cccGGTTTTTTT‑ATG‑AAAAATTTAGGGAAAAAAAGGAGAGGAAAAAAGGAAAATTTCTAAATAAATGGGGTAACTA‑AGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                           <  2:516312/149‑1 (MQ=255)                                                               tttt‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                  >  1:554659/1‑150 (MQ=255)                                                               tttt‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                  >  1:170639/1‑150 (MQ=255)                                                                 tt‑AAG‑AAAAAATTAGGAAAAAAAAGGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGGGGTAATTA‑AGAGATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAg                                                                                <  2:24048/150‑1 (MQ=255)                                                                     tG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                            >  1:192402/1‑151 (MQ=255)                                                                        aaaaatttAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGGGTTAATTG‑AGAGATGGGGATTGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                            <  1:1039495/143‑1 (MQ=255)                                                                        aaaaaaTTATGGAAAAAAATGAAAGGAAAAAAGAAAAAATTCTAAAAAAATGGGGTAACTT‑AGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                          <  2:690821/148‑1 (MQ=255)                                                                        aaaaaaTTATGGAAAAAAAGGAAGGGAAAAAAGAAAAAATCCAAAAAAAATGGGATAATTA‑AGAAATTAGGTTTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                              <  1:1166822/144‑1 (MQ=255)                                                                         agaaATTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTCAAAAAAGGGGGTAACTT‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc                                                                         <  2:1164929/151‑1 (MQ=255)                                                                          gaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                        >  2:617452/1‑151 (MQ=255)                                                                            aaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                      >  2:844572/1‑151 (MQ=255)                                                                            aaTTAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCAAAAAAAGTGCGGTAACTT‑AGAGATTAGGTTGGCGGAAAAAAAAATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                      <  1:204634/151‑1 (MQ=255)                                                                              ttATGGAAAAAAAGAAGAGAAAAAAAGAGAAAATTCTAAATAATTGGGTTATCTA‑AGAGATTAGGATGGCGGAAAATAAAAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt                                                                     <  2:62350/150‑1 (MQ=255)                                                                               tATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                    >  2:910742/1‑150 (MQ=255)                                                                                aTGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                   >  2:92785/1‑150 (MQ=255)                                                                                aTGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                   >  2:1053201/1‑150 (MQ=255)                                                                                aTGGAAAAAAATGACAGGAAAAAAGAAAAAATCCTAAAAAAAGGCGGTAATTA‑AGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                   <  2:1099973/150‑1 (MQ=255)                                                                                      aaaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                >  2:249899/1‑147 (MQ=255)                                                                                        aaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                          >  2:1073979/1‑151 (MQ=255)                                                                                         aaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                          >  1:675579/1‑150 (MQ=255)                                                                                          aggaaagaaaaaaaGGAAAAATTCTAAATAAATGGGGTAATTA‑AGATATAAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa                                                        <  2:140114/142‑1 (MQ=255)                                                                                               tgggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                   >  1:772294/2‑151 (MQ=255)                                                                                                gggAAAAGTGATAATTTCAAAATAAATGGGTTATGTG‑AGGGATTGGGTTTGCTAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                  <  2:677652/151‑1 (MQ=255)                                                                                                ggaaaaaaGGAAAAATTCTAAAAAAATGCGGTAACTT‑AGAGATTAGGTTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                    <  2:833224/146‑1 (MQ=255)                                                                                                ggaaaaaaGAAGAAATTCTAAAAAAATGGGTTAACTA‑AGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                    <  2:1009791/146‑1 (MQ=255)                                                                                                ggaaaaaaGAAAAAATTCTAAAAAATCGCGGTAAATA‑AAAAATTAGAAGGAGGAAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                   <  1:919039/147‑1 (MQ=255)                                                                                                      aaGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                            >  2:1121993/1‑151 (MQ=255)                                                                                                       aGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                           >  1:117547/1‑151 (MQ=255)                                                                                                          aaaaaTTCTCAAAAAATGCGTTAACTT‑AGAGATTAGGATGGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                        <  2:1180619/149‑1 (MQ=255)                                                                                                            aaaTTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                      >  1:309365/1‑151 (MQ=255)                                                                                                            aaaTTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                      >  2:171423/1‑151 (MQ=255)                                                                                                            aaaTTCTCAAAAATTGCGGTAACTT‑AAAAATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                      <  2:708914/151‑1 (MQ=255)                                                                                                             aaTTCTCAAAAAATGCGTTAATTT‑AGAGATAAGGTTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                     <  1:546183/151‑1 (MQ=255)                                                                                                                   cAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                               >  2:767453/1‑151 (MQ=255)                                                                                                                        aaggCGATAAATA‑AGAGATTAGGTTTGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                          <  1:409914/148‑1 (MQ=255)                                                                                                                        aaTGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                          >  2:1176029/1‑151 (MQ=255)                                                                                                                          tGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                        >  1:1155613/1‑151 (MQ=255)                                                                                                                           gCGGTAACTT‑AGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                        <  1:274210/150‑1 (MQ=255)                                                                                                                                aaCTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                  >  2:240134/1‑151 (MQ=255)                                                                                                                                  cTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                 >  2:391570/1‑150 (MQ=255)                                                                                                                                           ttAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag          <  2:1046500/149‑1 (MQ=255)                                                                                                                                            aaGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                 <  1:762929/140‑1 (MQ=255)                                                                                                                                              ggtttGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAg     <  2:534678/148‑1 (MQ=255)                                                                                                                                                  gggggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCt  <  2:946015/147‑1 (MQ=255)                                                                                                                                                   ggggAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGc   <  2:929615/146‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC‑TTTTATGCAGCCCGGCTTTTTT‑ATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCT  >  NC_000913/5012‑5302

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