Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,098	0	A	C	44.7%	99.5 / 10.5	76	Y66S (TAC→TCC)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (23/19);  new base C (0/34);  total (23/53)
Fisher's exact test for biased strand distribution p-value = 2.76e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGG  >  NC_000913/1490951‑1491248                                                                                                                                                    |                                                                                                                                                       ttttACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATa                                                                                                                                                        <  1:383941/148‑1 (MQ=255)       tGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAAcccc                                                                                                                                               >  2:1099864/1‑151 (MQ=255)        ggTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAAccc                                                                                                                                                <  1:139031/149‑1 (MQ=255)         gTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCccgc                                                                                                                                             >  2:1080870/1‑151 (MQ=255)            accGCAGCGTATTGCGCAATGGATTTTCGTGTTTTCTTCCCCCGTAGCGAGCGCCCCCCCATCAACATGATTCATTTTTCCTGTGCCATCCCAATTCTCTTGCTGATGGTCTTTCGTCTTTTTTTAAGGCTAAAATCCCCAACCccgccg                                                                                                                                           <  1:771271/148‑1 (MQ=255)               gCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGAtt                                                                                                                                        >  2:631802/1‑150 (MQ=255)                  gCGTATTGCGCAATGGAGTTTCTTGGTTTCTTCCCCCGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATAc                                                                                                                                    <  2:944225/151‑1 (MQ=255)                    gTATTGCGCAATGGAGTTTCGTGTTTTCTTCCCACGTAGCGATCGCCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAGGGCTGAAATACCCAACCCCGCCGATta                                                                                                                                       <  1:1196874/146‑1 (MQ=255)                      aTTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGCCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTa                                                                                                                                 <  1:780525/150‑1 (MQ=255)                         gcgcAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGcc                                                                                                                             >  2:204729/1‑151 (MQ=255)                          cgcAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGcc                                                                                                                             >  1:1108391/1‑150 (MQ=255)                            cAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCaa                                                                                                                           >  1:337288/1‑150 (MQ=255)                            cAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTAAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCaa                                                                                                                           <  1:1129382/150‑1 (MQ=255)                                  aGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCgat                                                                                                                    >  1:165167/1‑151 (MQ=255)                                   gTTTGGTGGTTTCCTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTCCCGGGGGCATCTCAATTCTCTTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTAAAACACCCACCCCCGCCGATTATACCTAAGCCAAAACCgat                                                                                                                    <  2:787264/150‑1 (MQ=255)                                      tCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCTTTCTTCTTTTGTAAAGGCTAAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCgatgat                                                                                                                 <  1:699831/150‑1 (MQ=255)                                            tttCTCCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTAAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACgg                                                                                                            <  1:1173570/149‑1 (MQ=255)                                              tcttcCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGAc                                                                                                         <  1:792487/150‑1 (MQ=255)                                                  cccACGTAGCGATCGGCCCCTCATCAACATGTTTCATGTTCCCTGTGGCATCTCAATCCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTg                                                                                                       <  1:56680/148‑1 (MQ=255)                                                  cccACGTAGCGATCGGCCACTCATCAACATGATTCATTTTTCCTGTGGCATCTAAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACt                                                                                                        <  1:1080816/147‑1 (MQ=255)                                                  cccACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACt                                                                                                        <  2:697831/147‑1 (MQ=255)                                                        tAGCGTTCGGCCCCTCTTCAACATGTTTCATGTTTCTTGTGGCAGCTAAATTCTCGTGTTGATGGTCGTTCGTCTTTTGTTAAGGCTAAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCAttt                                                                                               <  1:563480/150‑1 (MQ=255)                                                         aGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTgg                                                                                             >  1:1158930/1‑151 (MQ=255)                                                         aGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTgg                                                                                             >  1:1033755/1‑151 (MQ=255)                                                          gCGATCGGCCACTCATCAACATGTTTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTggg                                                                                            <  1:705017/151‑1 (MQ=255)                                                            gATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTAAATTCTCGTGCTGTGGGTCTTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTggg                                                                                            <  1:1153741/149‑1 (MQ=255)                                                             aTCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGAc                                                                                          >  1:398668/1‑150 (MQ=255)                                                                ggCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACAttt                                                                                      >  2:85473/1‑151 (MQ=255)                                                                  ccACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTgg                                                                                    <  2:295328/151‑1 (MQ=255)                                                                  ccACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTgg                                                                                    >  1:918524/1‑151 (MQ=255)                                                                   cACTCATCAACATGTTCAATTTTTCCTGGGGCATCTAAATTCTCTTGCTGATGGTCGTTCGTTTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGt                                                                                   <  1:581363/151‑1 (MQ=255)                                                                   cACTCATCAACATGATTCATGTTTACTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGt                                                                                   >  2:1183044/1‑151 (MQ=255)                                                                        atcaACATGATTCTTGTTTCCTGTGGCATCTAAATTCTCGTGCTGATGGTCGTTCGTCTTTGGTAAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGAtt                                                                               <  1:351933/150‑1 (MQ=255)                                                                             cATGATTCATGTTTCCTGTGGCCTCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCAAACCCCGCCGTTTATACCTAAGCCAAAACCGATGATGACGGGACTGACGCATTTGGGACATTTGGGGATTTAtctc                                                                         >  2:927843/1‑151 (MQ=255)                                                                             cATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAtctc                                                                         <  2:430453/151‑1 (MQ=255)                                                                                gATTCATGTTTCTTGTGCCTTCTAAATTCCCTTGCTGAGGGTCGTTCGTCTTGTTTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAtct                                                                          <  1:758406/147‑1 (MQ=255)                                                                                gATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTg                                                                      <  2:294741/151‑1 (MQ=255)                                                                                  ttCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCt                                                                       <  2:454389/148‑1 (MQ=255)                                                                                    cATGTTTCCTGTGACTTCTCAATCCTCGTGCTGATGGTCGTTCGTCTTTTTTAAAGGTTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGtt                                                                    <  2:25332/149‑1 (MQ=255)                                                                                        tttCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTAtt                                                                <  1:1139138/149‑1 (MQ=255)                                                                                         ttCCTGTGGCATCTAAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTgc                                                              <  1:583691/150‑1 (MQ=255)                                                                                           ccTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCt                                                           <  1:370396/151‑1 (MQ=255)                                                                                              gtgGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGcc                                                        <  2:1176342/151‑1 (MQ=255)                                                                                              gtgGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGcc                                                        >  2:1084461/1‑151 (MQ=255)                                                                                               tgGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                                         <  1:133168/149‑1 (MQ=255)                                                                                                ggCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCa                                                       <  1:1099810/150‑1 (MQ=255)                                                                                                ggCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAg                                                      <  2:964599/151‑1 (MQ=255)                                                                                                ggCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAg                                                      <  1:612699/151‑1 (MQ=255)                                                                                                ggCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                                         <  2:333392/148‑1 (MQ=255)                                                                                                   aTCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGt                                                     <  1:1042366/149‑1 (MQ=255)                                                                                                        aaTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggt                                              <  2:885470/151‑1 (MQ=255)                                                                                                        aaTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                                                  <  2:150416/147‑1 (MQ=255)                                                                                                        aaTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttg                                                <  2:402571/149‑1 (MQ=255)                                                                                                        aaTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttg                                                <  1:305712/149‑1 (MQ=255)                                                                                                                   tGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgatg                                     <  1:238215/149‑1 (MQ=255)                                                                                                                   tGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTgatgat                                   >  2:648353/1‑151 (MQ=255)                                                                                                                    gATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGatgatg                                  <  2:986749/151‑1 (MQ=255)                                                                                                                     aTGGTCGTTCGTCTTTTGTTAAGGCTAAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgatga                                    <  1:643993/148‑1 (MQ=255)                                                                                                                       ggTCGTTCGTCTTTTGTTAGGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGatgatg                                  <  1:647669/148‑1 (MQ=255)                                                                                                                       ggTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGta                                <  1:86871/150‑1 (MQ=255)                                                                                                                            ttcttcTTTTTTAAAGGCTAAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAAcc                           <  2:1041059/150‑1 (MQ=255)                                                                                                                               gtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCgggg                       >  1:611327/1‑151 (MQ=255)                                                                                                                                tCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGc                      >  1:932304/1‑151 (MQ=255)                                                                                                                                  ttttGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCa                     >  2:219438/1‑150 (MQ=255)                                                                                                                                      gTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCg                >  2:837709/1‑151 (MQ=255)                                                                                                                                         aaGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTgg             <  2:189735/151‑1 (MQ=255)                                                                                                                                         aaGGCTAAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTg              <  2:843204/150‑1 (MQ=255)                                                                                                                                           ggCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGtt           <  2:312370/151‑1 (MQ=255)                                                                                                                                           ggCTAAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGtt           <  1:173309/151‑1 (MQ=255)                                                                                                                                             cTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTg         >  2:828242/1‑151 (MQ=255)                                                                                                                                              tGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTg         >  1:151586/1‑150 (MQ=255)                                                                                                                                                aaaTCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGc        <  2:398689/149‑1 (MQ=255)                                                                                                                                                aaaTACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTg         <  2:526066/148‑1 (MQ=255)                                                                                                                                                 aaTCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgtt     <  2:729774/151‑1 (MQ=255)                                                                                                                                                 aaTCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgtt     <  2:723963/151‑1 (MQ=255)                                                                                                                                                 aaTACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgtt     >  1:554036/1‑151 (MQ=255)                                                                                                                                                 aaTACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgtt     <  2:165178/151‑1 (MQ=255)                                                                                                                                                    aCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgtttgg  >  1:600177/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                       TTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGG  >  NC_000913/1490951‑1491248

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