Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	41.4%	107.3 / 17.6	59	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (32/2);  new base C (0/24);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 2.53e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAA  >  NC_000913/993770‑994061                                                                                                                                                     |                                                                                                                                                 aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCccacca                                                                                                                                               <  1:526903/151‑1 (MQ=255)   cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCaccacc                                                                                                                                             <  1:599568/151‑1 (MQ=255)     gACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGAGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGc                                                                                                                                           >  2:477162/1‑151 (MQ=255)      aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGcc                                                                                                                                          >  2:694898/1‑151 (MQ=255)       cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGcc                                                                                                                                          >  1:548492/1‑150 (MQ=255)       cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGcc                                                                                                                                          >  2:39089/1‑150 (MQ=255)         cgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCcac                                                                                                                                        >  2:1087475/1‑150 (MQ=255)         cgcATCCTGAAACATCATGCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTCCTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCca                                                                                                                                         <  2:555478/149‑1 (MQ=255)            aTCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAAc                                                                                                                                    >  2:72241/1‑151 (MQ=255)                  aaaCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                              >  2:180536/1‑151 (MQ=255)                    aCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGc                                                                                                                            >  1:619104/1‑151 (MQ=255)                     catcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGcc                                                                                                                           >  1:638377/1‑151 (MQ=255)                        catTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                         >  1:312666/1‑150 (MQ=255)                         atTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGa                                                                                                                        >  1:364267/1‑150 (MQ=255)                             gcgTGTCTTCTTGATTTTCAGCCGGCGGTGTGGGCCCCGCAAACACTCCGCCGGCGGTTGGTTTTTCCAGCCCGGCCAGCAGGCGCAGCAGGGTACTTTTCCCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGa                                                                                                                        <  1:1112085/146‑1 (MQ=255)                              cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                  >  1:218007/1‑151 (MQ=255)                              cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                  >  1:47109/1‑151 (MQ=255)                               gtgtCTTCCTGGATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                  >  2:758616/1‑150 (MQ=255)                               gtgtCTCCCTGATTTTCAGCCAGGGGGGGGGGGCCCGCAAACACTCCGCCGGGGGTGGGCTTTCCCAGCCCGGCCAGCAGGCGCAGCAGGGTATTTTTCCCCCCC‑CCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAat                                                                                                                      <  1:935011/146‑1 (MQ=255)                                tgtCTTCTTAAATTTCAGCCGGGGGGGTGGCCCCCGCCACCCCACCCCCGGCGTTGGGCGTTCCCACCCCTGCCAGCGGGCACAGCGGGTTTTTTTTCCCACCC‑CCGGGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGaa                                                                                                                       <  2:129949/144‑1 (MQ=255)                                  tCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                              >  1:148072/1‑151 (MQ=255)                                   cTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                             >  1:413412/1‑151 (MQ=255)                                      ccaaaaTTCCAGCCAGGGGTGGGGGCCCCCCTAACCCCCCCCGGGGGGTTGGTTTTCCCACCCCGGCGGGCGGGCGCAGCAGGGTTTTTTGCCCCCCA‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg                                                                                                          <  2:839762/147‑1 (MQ=255)                                        tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                        >  1:283916/1‑151 (MQ=255)                                           aTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                     >  1:428382/1‑151 (MQ=255)                                           aTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                     >  1:52848/1‑151 (MQ=255)                                                  ccgggggTGTGGTGCCCGCAACCCCATCGCGGGGGGTTGGCTTTCCCAGCCCTGCAGGCAGGCGCAGGGGGGTACTTTTCCCCCCC‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                   <  1:358838/141‑1 (MQ=255)                                                   cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                             >  1:211290/1‑151 (MQ=255)                                                     ggggTGGGGTGCCCGCTACCACACCGCCTGTGGTGGGCTTTCCCAGCCCTGCCAGGAGGCGCAGCGGGGTTTTTTTGCCCCCA‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                   <  1:121239/141‑1 (MQ=255)                                                      cGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                           >  2:274784/1‑150 (MQ=255)                                                          gtggggCCCGATAACCCACCGCGTGCGTTTGGCGTTCCCAGCCCTCCAAACAGGCGAAGCGGGTTTCTTTTCCCACAA‑CCGCCGCGGCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                    <  2:813542/137‑1 (MQ=255)                                                            ggtcCCCACTAACACACCGCCTGGGGTTGGCTTTTCCAGCCCTGCCAGCGGGCGCAGCAGGTTTTTTTTCCCCCCC‑CCGTTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                              <  1:1084994/141‑1 (MQ=255)                                                                      accccaccGCCTGGGGTGGGCTTTTCCAGCCCTCCCAGCGGGCGCAGGGGGTTTTTTTCCCCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                              <  2:490183/140‑1 (MQ=255)                                                                        cccaccGCCGGGGGTTGGCTTTCCCAGCCCTGCCAGCGGCCGCGGAGGGGTCTTTTTCCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                          <  1:133249/144‑1 (MQ=255)                                                                        cccaTCGCCGGGGTTTGGCTTTTCCAGCCCTGCCAGCGGGCGCGGGAGGTTCCTTTTCCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                        <  1:207061/149‑1 (MQ=255)                                                                             cGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                   >  1:60879/1‑151 (MQ=255)                                                                               ccTGGGGTTGGCTTTCCCAGCCCGGCCAGCAGGCGCGGCAGGTTACTTTTGCCACCC‑CCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                                  <  1:8170/150‑1 (MQ=255)                                                                               ccTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                 >  1:770078/1‑151 (MQ=255)                                                                                 tGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                >  1:349359/1‑150 (MQ=255)                                                                                  gCGGTTGGCTTTTCCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTTCCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                <  1:681025/149‑1 (MQ=255)                                                                                    ggTTGGCTTTCCCAGCCCCGCCAGCAGGCCCAGCAGGTTCTTTTTCCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                              <  1:558645/149‑1 (MQ=255)                                                                                      ttGGCTTTCCCAGACCTGCCGGCGGGCGCAGGAGGTTCTTTTCGCCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                           <  1:970976/150‑1 (MQ=255)                                                                                        ggCGTTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCCCCC‑CCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                        <  1:351863/151‑1 (MQ=255)                                                                                               ccAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                 >  1:105446/1‑151 (MQ=255)                                                                                                 aGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                               >  2:699438/1‑151 (MQ=255)                                                                                                   cccTGCCGGCAGGCGAGGAAGGTTTCTTTCGCCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                  <  2:817885/145‑1 (MQ=255)                                                                                                   aCCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                             >  1:450400/1‑151 (MQ=255)                                                                                                   aCCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                             >  1:232109/1‑151 (MQ=255)                                                                                                    ccTGCCAGCAGGCCCAGCGGGGTCCTTTTCCCCCCC‑CCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                              <  1:646045/149‑1 (MQ=255)                                                                                                              ggCGCAGCGGGTTCTTTTTCCCACCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                    <  1:390601/149‑1 (MQ=255)                                                                                                                       ggttatttttTCCCCCC‑CCTCAGCGCCCCGCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                         <  1:545921/145‑1 (MQ=255)                                                                                                                       gggtcttttcGCCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAg                          <  1:134624/144‑1 (MQ=255)                                                                                                                         gttattttGCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                        <  1:922636/147‑1 (MQ=255)                                                                                                                              ttccccaccc‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                   <  2:494047/146‑1 (MQ=255)                                                                                                                                  ccacca‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt              >  2:165276/1‑151 (MQ=255)                                                                                                                                   ccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt              <  1:925331/145‑1 (MQ=255)                                                                                                                                      ca‑ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc          >  1:122043/1‑151 (MQ=255)                                                                                                                                        ccgccTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct         <  2:647925/147‑1 (MQ=255)                                                                                                                                              gcgcCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct      <  2:759996/148‑1 (MQ=255)                                                                                                                                               cgcccccccccCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa   <  2:975149/147‑1 (MQ=255)                                                                                                                                                ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa  >  1:814522/1‑150 (MQ=255)                                                                                                                                                ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa  >  2:1100730/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAA  >  NC_000913/993770‑994061

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