Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,080,837	0	G	T	43.7%	62.7 / 37.2	71	intergenic (+24/‑519)	putP/efeU	proline:sodium symporter/ferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base G (12/28);  new base T (31/0);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 1.24e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACAC  >  NC_000913/1080691‑1080987                                                                                                                                                   |                                                                                                                                                       aGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTtttgtttg                                                                                                                                                    <  1:625418/151‑1 (MQ=255)     ttGGGATTGTAGTGTTTAGTTTGCTTGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCtt                                                                                                                                                >  1:102093/1‑151 (MQ=255)     ttGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTgg                                                                                                                                                   >  1:899252/1‑148 (MQ=255)     ttGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTGTTTGGCtt                                                                                                                                                >  2:1301574/1‑151 (MQ=255)       gggATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTc                                                                                                                                               <  2:443239/150‑1 (MQ=255)       gggATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTc                                                                                                                                               <  1:609530/150‑1 (MQ=255)        ggATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGtt                                                                                                                                                >  2:1001809/1‑148 (MQ=255)            tGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTcagca                                                                                                                                           >  2:533473/1‑149 (MQ=255)             gTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCagcagc                                                                                                                                         <  1:891452/150‑1 (MQ=255)               aGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCCGCAGCggg                                                                                                                                      >  2:275124/1‑151 (MQ=255)                gtgtTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCggg                                                                                                                                      >  2:150047/1‑150 (MQ=255)                  ggttAGTTTGCGGGGTAAAGCGCCGTCAGCCGCGATGCAAAAACGCTTTGCCGACGCCGATGCGCACTAACATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTg                                                                                                                                   <  1:1052245/149‑1 (MQ=255)                   tttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCCTCCCCAACggggt                                                                                                                                    >  1:971944/1‑147 (MQ=255)                       gTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTGTTTTGCTTTCGCCGCggggt                                                                                                                                    >  1:213676/1‑143 (MQ=255)                            cTGGGTAAAGCGCCGTCAGGGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTCGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCtt                                                                                                                          <  1:1015964/150‑1 (MQ=255)                                gTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAAtgt                                                                                                                     <  1:1071012/151‑1 (MQ=255)                                  aaaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAGCAGCGGGGTGCGCTCCCtt                                                                                                                          >  1:763489/1‑144 (MQ=255)                                  aaaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAAtgtg                                                                                                                    <  1:957663/150‑1 (MQ=255)                                  aaaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAAtgtg                                                                                                                    <  2:197874/150‑1 (MQ=255)                                    aGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCCTCAACAGCGGGGTGCGCTCCCCTAATGTGcc                                                                                                                  >  2:1205112/1‑150 (MQ=255)                                    aGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCt                                                                                                                 <  2:521281/151‑1 (MQ=255)                                     gcgcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGTTTCAGCAGCGGGGTGCGCCCCCTTTAtgt                                                                                                                     >  1:85986/1‑146 (MQ=255)                                       gcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGc                                                                                                              <  2:1290143/151‑1 (MQ=255)                                         cGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGTTTTTGCAGCCGGGGGCGCTCCCCTTAtgtggccccccc                                                                                                            >  2:763103/1‑143 (MQ=255)                                             aGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCCTTCGCAACGGGGTTCGGCCCCTTAATTGGCCTccccc                                                                                                            >  1:16378/1‑146 (MQ=255)                                               cggcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCGCTCCCTTAATGTGCCTCGCCCTATaa                                                                                                       >  2:1225936/1‑150 (MQ=255)                                                   gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTTAGCAGCGGGGTGCGCTCCCCTTATTTGCCTccccctcttaa                                                                                                      >  2:1065422/1‑140 (MQ=255)                                                    aTGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCCGCCACGGGGTGCGCCCCCTTAATGTGCCTCCCCCTATaaaat                                                                                                    >  1:133416/1‑146 (MQ=255)                                                       cAAAAACGCTTTGCCGCGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGACTTGTTTGGCGTTAGCACCCGGTTGCTGTCCCTTTTTGTGCCTCGCCATATAAAtattattg                                                                                              >  2:865518/1‑144 (MQ=255)                                                        aaaaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCCGCAGCGGGGTGGGCTCCCCTTATTTGGCTCGCCATATAAAATTAATgg                                                                                              >  1:1227806/1‑150 (MQ=255)                                                        aaaaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCCTCAGCAGCGGGGTGCGCCCCCTTAATGTGCCCCCCCCTATAAATTGAATgg                                                                                              >  1:328059/1‑150 (MQ=255)                                                         aaaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTg                                                                                            <  1:1050554/151‑1 (MQ=255)                                                            aCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAGAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTTTGCAGCGGGGTTCGGTCCCCTAATTTGGCTCGCCCTATAAATTGaa                                                                                                 >  1:567523/1‑143 (MQ=255)                                                            aCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCTGCAGCGGGGTGTGCTCCCTTAATGTGCCTCGCCCTATAAATTTGATgggggcg                                                                                         >  2:1200701/1‑146 (MQ=255)                                                            aCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAg                                                                                         <  2:419534/151‑1 (MQ=255)                                                                tttGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGATTTGTTTTGCTTCAACCACGGGTTGCGCTCCCCTTATGTGACCCGCCCTATAAAATTAATAGGGCAagg                                                                                       >  2:186195/1‑149 (MQ=255)                                                                      gAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGGTTCGCTCCCCTTATGTGCCTCCCCCTATAAATTGAATGGGGCCGGGAgcg                                                                                   >  2:384111/1‑147 (MQ=255)                                                                      gAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGGTTCAGCAGCGGGGTGCGCTCCCTTAATTTGCCTCCCCCTATAAATTTAATGGTGCAGGGAgcgc                                                                                  >  2:457401/1‑148 (MQ=255)                                                                      gAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTGCTTCCGCAGCGGGGTGGGCTCCCCTTATGTGGCTCCCCCTATAAATTTGATGGGGGAGGGGgcgcgc                                                                                >  1:1148918/1‑150 (MQ=255)                                                                         gCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAggg                                                                            <  2:280259/151‑1 (MQ=255)                                                                            gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGGTTCAGCAGCGGGGTGCGCCCCCTTTATGTGCCTCGCCCTATAAATTTAATGGGGGAGGGGGGGCGCAgggggg                                                                         >  1:673258/1‑150 (MQ=255)                                                                            gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGggc                                                                         <  2:1091566/151‑1 (MQ=255)                                                                             aTGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCCTCACGAGCGGGGTTGGCTCCCCTAATGTGCCTCGCCCTATAAATTTAATGGTGCAGGGAGCGCGGCgggggg                                                                         >  1:217930/1‑149 (MQ=255)                                                                               gcgcACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGggcggc                                                                      <  1:543916/151‑1 (MQ=255)                                                                                cgcACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGcc                                                                     <  2:645324/151‑1 (MQ=255)                                                                                   aCTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCATCAGCGGGTTTCTGACCCCTAATGTGTCTCTCCACATAAAATGAATGGGGGGGGGAGCGGCCAGGGGGGGGccca                                                                   >  2:1025452/1‑148 (MQ=255)                                                                                   aCTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAACGGGGTGCGCTCCCTTAATGTGGCTCGCCATATAAATTTAATGGTGCAGGGAGCGCGCAGGGGGGGGCCaa                                                                   >  2:1287104/1‑150 (MQ=255)                                                                                    cTAGCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAAt                                                                  <  1:1229867/150‑1 (MQ=255)                                                                                     tATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGTTTCAGCAGGGGGTTTCGGTCCCCTAATTTGGCTTGCCATTTAAAAAGAATGGGGTAGGGGGCGGGCCGGGGGGGGCCa                                                                    >  2:270530/1‑147 (MQ=255)                                                                                     tATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTGGTTCCGCAACGGGGTGGGCCCCCTTAAAATGGCCCCCCCTCTAAATTGAATGGTGCAGGGGGCCGGCAGGGGGGGGGCCATc                                                                 >  1:1234770/1‑150 (MQ=255)                                                                                           tCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGCTCCGCAGCGGGGTTCGCCCCCCTTATGTGGCTCCCCCTTTAAATTTAATGGGGCAGGGGGGGCGCCGGGGGCGGCCAAAcgccccc                                                           >  1:984733/1‑150 (MQ=255)                                                                                                 ccgccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCct                                                    <  1:255153/151‑1 (MQ=255)                                                                                                  cgccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAACAGCCGGGTGCGGCCCCCTTATGTGGCTCGCCCTATAAAATGAATTGGGCAGGGGGCGCGCCGGGGGGGGGCCAACGCCGcccccccc                                                     >  1:1098823/1‑149 (MQ=255)                                                                                                    ccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCAGCAGCGGGGTGCGCTCCCCTTATGTGCCTCGCCATATAAATTGTAAAGGGCCGGGGGGGCGCCGGGGGGGGCCCAACGCccccccccccc                                                    >  1:591482/1‑147 (MQ=255)                                                                                                    ccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCctgct                                                 <  2:298207/151‑1 (MQ=255)                                                                                                        cACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGGTGCGCCCCCCTAATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCCCCCCCGCCCCCCGCTGTcc                                             >  2:337630/1‑151 (MQ=255)                                                                                                           ggTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTTAGCAGCGGGTTGCGCTCCCCTAATGTGCCTCGCCCTATAAAATTAATTGGGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCCCCCCCTGCTGTccc                                            >  1:56127/1‑149 (MQ=37)                                                                                                           ggTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCgg                                          <  2:673281/151‑1 (MQ=255)                                                                                                           ggTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTGGTTCAGCAGCGGGGTGCGCTCCCCTTATGTGGCTCCCCCTATAAATTTAATGGTGGAGGGGGGCCGCCGGGGGGGGCCAATCCCCGccccccccccccgttcccg                                          >  1:135243/1‑140 (MQ=255)                                                                                                            gTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCGCTCCCTTTATGTGCCTCGCCATATAAATTTAATTGTGGAGGGAGCGCGCAGGGGGCGGGCCATCGCCGCCCCCCCCCGCTGTccccg                                          >  2:991434/1‑148 (MQ=255)                                                                                                              tGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTTGGTTCAGCAGCGGGGTGCGCTCCCCTAATGTGCCTCGCCATATAAAATGAATTGGGCCGGGGGCGCGCCGGGGGGGGCCAATCGCCGCCCCCCCCCGCTGTCCCGGcc                                        >  2:481215/1‑150 (MQ=255)                                                                                                               gCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCtt                                      <  1:149612/151‑1 (MQ=255)                                                                                                                  ggAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCgg                                   <  2:861519/151‑1 (MQ=255)                                                                                                                    aaaGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTTATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCCTCgggg                                 >  2:694828/1‑151 (MQ=255)                                                                                                                         tAAGGGACTTAGCCTGCGGCGGGTTTTTTTGGGTTCCGCAGCGGGGTGGGCTCCCCTAATTTGGCTCGCCCTATTAATTTAATGGTGGCGGGAGCGCGCCGGGGGGGGCCAAACCCCGCCGCCCCCCGCTGTCCCCGCCCTCGGGGGACg                             >  1:600352/1‑150 (MQ=255)                                                                                                                            gggACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCtt                          <  1:1017037/150‑1 (MQ=255)                                                                                                                              gACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAg                       <  1:1029977/151‑1 (MQ=255)                                                                                                                                  tAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCg                     <  2:247876/149‑1 (MQ=255)                                                                                                                                        gcggcgGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTg               <  2:1227881/149‑1 (MQ=255)                                                                                                                                        gcggcgGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTg               <  2:961255/149‑1 (MQ=255)                                                                                                                                          ggcggTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGc           <  2:1049425/151‑1 (MQ=255)                                                                                                                                             ggTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccac        >  2:792648/1‑151 (MQ=255)                                                                                                                                                   gtttgGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAacac  >  1:60892/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                       AGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCGCTCCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACAC  >  NC_000913/1080691‑1080987

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