Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,807,369	0	A	C	40.6%	101.4 / 14.8	64	T80P (ACC→CCC)	proX	glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein
Reads supporting (aligned to +/- strand):  ref base A (28/10);  new base C (0/26);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 4.44e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACA  >  NC_000913/2807221‑2807513                                                                                                                                                     |                                                                                                                                                 tGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAAcc                                                                                                                                                >  1:753718/1‑151 (MQ=255) tGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAAcc                                                                                                                                                >  2:724704/1‑151 (MQ=255)    tAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTc                                                                                                                                             >  1:314190/1‑151 (MQ=255)          aGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCAccgcc                                                                                                                                       >  2:151443/1‑151 (MQ=255)            ttCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCAccgccg                                                                                                                                      >  1:178159/1‑150 (MQ=255)             tCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTg                                                                                                                                    >  1:501930/1‑151 (MQ=255)             tCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTCGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTg                                                                                                                                    >  2:695875/1‑151 (MQ=255)                      cATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACg                                                                                                                           >  1:478147/1‑151 (MQ=255)                       aTCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGc                                                                                                                          >  2:344558/1‑151 (MQ=255)                         cACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGcc                                                                                                                         >  2:429378/1‑150 (MQ=255)                          aCTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCAc                                                                                                                       >  1:943195/1‑151 (MQ=255)                            tGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTg                                                                                                                     >  1:474115/1‑151 (MQ=255)                              aagaaCCCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTACCCCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACt                                                                                                                      <  2:872189/148‑1 (MQ=255)                                 aacccTTCCAGACGCTGCTGGTCAGTCGTGCGCGGGAGAAATTAGGTTATCCCGTAAACAACCCCAGCGAAGTAGTTTCAAACGTTGGCTACACCTCGCTTGCTTCCGCAGATCCACCCTTCACCGCCGTGAACTGGACGCCACTGCATg                                                                                                                 <  1:736901/147‑1 (MQ=255)                                 aacccTCCCAGCCGCTGCTGGTCAGTCGTGCGCTGGAAAAATTAGGTTATCCCGTCAACAACCCCAGCAAAGTAAATTACAACGTTGGCTACCCCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATg                                                                                                                 <  1:739836/147‑1 (MQ=255)                                    ccTTCCAGACGCTGCTGGTCGGTCGTGCGCTGGAAAAATTAGGTTATCCCGCCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTCCTTCCGGCGATCCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGacaa                                                                                                             <  2:38687/151‑1 (MQ=255)                                        ccAGACGCTGCTGGTCAGTCGTGCGCGGGAGAAATTAGGTTATCCCGTCAACAACCCCAGCGAAGTAGATTACACCGTTGGCTACCCCTCGCTTGCTCCGGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGacaaca                                                                                                           <  1:557388/149‑1 (MQ=255)                                             cGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATCCCGTCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTACCCCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTAc                                                                                                      <  2:829713/149‑1 (MQ=255)                                                  ctgGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTg                                                                                               >  2:949437/1‑151 (MQ=255)                                                       cAGTCGTGCGCTGGAGAAATTAGGTTATCCCGTAAAAAACCCCAGCGAAGTAGATTACAACGTTGGCTACCCCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGt                                                                                          <  2:702669/151‑1 (MQ=255)                                                         gTCGTGCGCTGGAGAAATTAGGTTATCCCTTCAACAACCCCAGCGAAGTAGATTACAACGTGGGCTACCCCCCGCTTGCTCCGGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTg                                                                                         <  1:498358/150‑1 (MQ=255)                                                         gTCGTGCGCTGGAAAAATTAGGTTATCCCGTCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGt                                                                                          <  1:382957/149‑1 (MQ=255)                                                           cGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTGCAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCg                                                                                      >  2:213926/1‑151 (MQ=255)                                                            gTGCGCTGGAAAAATTAGGTTATCCCGTCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTACCCCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCg                                                                                      <  2:295960/150‑1 (MQ=255)                                                              gcgcTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGAt                                                                                    >  2:301292/1‑150 (MQ=255)                                                                  tGGAAAAATTTGTTTAGCCCGTCAACACACCAACGGAGGAAGATTAAAACGTTGGTTACCCCCCTCTGCCTCCCGCCGAACCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATaa                                                                                  <  1:384547/148‑1 (MQ=255)                                                                   ggAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGaa                                                                               >  1:453154/1‑150 (MQ=255)                                                                     aaaaaTTAGGTTATCCCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACCCCCCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAAt                                                                             <  1:280967/148‑1 (MQ=255)                                                                      gaAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAAttt                                                                           >  1:378722/1‑151 (MQ=255)                                                                        aaTTAGGTTATCCCGTCAACAAACCCGGCGAAGTACATGACAACGTTGGATACCCCCCGTTTGCTCCCGGCGATGCACCCTCCGCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTa                                                                         <  2:393697/151‑1 (MQ=255)                                                                         aTTAGGTTATACCATCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTCCCCATCGTTTCCTTCCGGCGAGGAAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTAt                                                                        <  2:214121/151‑1 (MQ=255)                                                                          ttAGGTTATACCGTCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTACCCCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTAt                                                                        <  2:392652/150‑1 (MQ=255)                                                                             ggTTATCCCGTCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTACCCCTCGCTTGCTCCCGGGGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTg                                                                    <  2:576641/151‑1 (MQ=255)                                                                             ggTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTg                                                                    >  2:59438/1‑151 (MQ=255)                                                                                 ataCCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAg                                                                 <  2:361884/150‑1 (MQ=255)                                                                                    ccGTCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTCCCGGCGATGCACCCTCCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAggg                                                               <  1:596913/149‑1 (MQ=255)                                                                                    ccGTAAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACCCCTCTCTTGCTTCCGGGGATGCACCCTCCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAggg                                                               <  2:735572/149‑1 (MQ=255)                                                                                         aacaaACCCAGCGAAGTAGATTCCAACGTTGGCTACCCCTCGCTTCCTCCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTAt                                                           <  1:74709/148‑1 (MQ=255)                                                                                            aaaCCCACCGATGTAGCTTACACCGTTGCCTCCCCCTCCCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTAt                                                           <  2:421376/145‑1 (MQ=255)                                                                                               cccAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACg                                                  >  2:459303/1‑151 (MQ=255)                                                                                                ccAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACgg                                                 >  2:895179/1‑151 (MQ=255)                                                                                                      aaGTAGATTAAAACGTTGGCTACCCCTCGCTTGCTCCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCgg                                            <  2:78425/150‑1 (MQ=255)                                                                                                         tAGATTACAACGTTGGCTACCCCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAGATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGcac                                         <  2:478161/150‑1 (MQ=255)                                                                                                         tAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGcac                                         >  2:435889/1‑150 (MQ=255)                                                                                                           gATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAgg                                      >  2:963930/1‑151 (MQ=255)                                                                                                             ttACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGt                                    >  1:228157/1‑151 (MQ=255)                                                                                                             ttACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGt                                    >  1:673570/1‑151 (MQ=255)                                                                                                              tACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGtt                                   >  1:291985/1‑151 (MQ=255)                                                                                                                 aaCGTTGGCTACCCCTCGCTTGCTCCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTAc                                 <  2:474129/150‑1 (MQ=255)                                                                                                                     tgggCTACCCCTCGCTTGCTCCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCt                               <  1:192072/146‑1 (MQ=255)                                                                                                                       ggCTACCCCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATc                           <  2:453167/150‑1 (MQ=255)                                                                                                                       ggCTACCCCTCGCTTGCTTCCGGCGATGCAACCTTCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGAt                            <  2:314340/149‑1 (MQ=255)                                                                                                                           accccTCGCTTCCTCCCGGGAATGCACCCTCCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTg                              <  1:911394/140‑1 (MQ=255)                                                                                                                            cacCCCGCTTGCTCCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGa                         <  1:192330/147‑1 (MQ=255)                                                                                                                              ccTCGTTTGCTCCGGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGaa                   <  2:375566/151‑1 (MQ=255)                                                                                                                              ccTCGCTTGTTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGa                    <  1:724678/150‑1 (MQ=255)                                                                                                                              ccTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGa                    >  2:631163/1‑150 (MQ=255)                                                                                                                              ccTCGCTTCCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATa                       <  2:943224/147‑1 (MQ=255)                                                                                                                                 cGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGaaaa                 <  2:359360/150‑1 (MQ=255)                                                                                                                                    ttgcttCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAAc                <  2:228161/148‑1 (MQ=255)                                                                                                                                       cttCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAAccgccg           >  1:221221/1‑150 (MQ=255)                                                                                                                                            ggCGATCCACCCTTCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAg      <  1:267364/150‑1 (MQ=255)                                                                                                                                                aTGCACCCTCCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACa  <  1:429366/150‑1 (MQ=255)                                                                                                                                                aTGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACa  >  2:680370/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACA  >  NC_000913/2807221‑2807513

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