Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,581	0	A	C	44.6%	110.4 / 21.3	65	V314G (GTA→GGA)	thiP	thiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (32/4);  new base C (0/29);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.83e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGC  >  NC_000913/73435‑73723                                                                                                                                                    |                                                                                                                                               gCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCCACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTAccc                                                                                                                                             >  2:301771/1‑150 (MQ=255)     cgcCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATc                                                                                                                                         <  1:222587/150‑1 (MQ=255)       cccGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATcgac                                                                                                                                      <  1:1052539/151‑1 (MQ=255)       cccGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATcgac                                                                                                                                      >  1:809144/1‑151 (MQ=255)        ccGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATcgacg                                                                                                                                     <  1:369677/151‑1 (MQ=255)          gCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGAt                                                                                                                                   >  2:499033/1‑151 (MQ=255)           cAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATc                                                                                                                                  >  2:992289/1‑151 (MQ=255)           cAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATc                                                                                                                                  >  2:955316/1‑151 (MQ=255)                  cgACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGcc                                                                                                                            >  1:231729/1‑150 (MQ=255)                  cgACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCa                                                                                                                           <  2:218666/151‑1 (MQ=255)                         tCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACgg                                                                                                                    >  2:788664/1‑151 (MQ=255)                            aTAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTgg                                                                                                                 >  2:608797/1‑151 (MQ=255)                                       atGGTCAGCACTCCGCACAATACCCCGCCCCCCAGCGAAAAACCCAGGGAGGCCCACAGCGCCGCCCACAGCCCCGGTTGGG‑CCAGCACTCCCGGCAACTGGCATTTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCaa                                                                                                           <  2:509111/146‑1 (MQ=255)                                         ggTCACCCCTACGCAAAACCCCCCGGCCGCCGGCGCAATACGAAACGAGGCCCACAGCGCCTCCCACAGCCCCGGTTGGG‑CCACCCCTCCCGGCAACGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAcagca                                                                                                      <  2:348720/149‑1 (MQ=255)                                            cAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACagcagcagc                                                                                                 >  2:349874/1‑151 (MQ=255)                                             acccccaCGCACAATCCCCCCCCCGCCAGCACAATCCGCAACGAGGCCCACAGCGCCTGCCACAGCACGGGTTGTG‑CCAGCCCTCCCGGAAATTGGCGTTTTCCCCCTTCGAGGATCACCGCCAGTAACGGTGGCAGCAAcagcagcag                                                                                                  <  2:777794/143‑1 (MQ=255)                                                   aCGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGc                                                                                           >  2:894753/1‑150 (MQ=255)                                                    cGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGcac                                                                                         >  1:250041/1‑151 (MQ=255)                                                      cacaATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGcaca                                                                                        >  1:456145/1‑150 (MQ=255)                                                       acaacaCCCCCGCCGCCAGCGAAACCCGCACGGAGGCCCACGGCGCCTGCCACAGCCCCGGTTGGG‑CCAGCACTCCGGGAAACGGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACaa                                                                                       <  2:296657/150‑1 (MQ=255)                                                       acaATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAt                                                                                      >  2:8139/1‑151 (MQ=255)                                                         aatcccccggcccccAGCGCAATCCGAAACGGGTCCCCCAGCCCCGCCCACAGCCCCGGTTGCC‑CCAGCACTTCCGGCAACGGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAt                                                                                      <  2:1138418/140‑1 (MQ=255)                                                         aatccccccccccccAGCCAAATACGCAACGGGTCCCACAGCGCTTCCCCCAGC‑CCGGTTGGGCCCAGCCCTCCCGGCAACGGGGGATTTCCCCCACCGACGTTCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAtt                                                                                     <  2:842777/137‑1 (MQ=255)                                                         aatcccccTGCCGCCAGCGCAATCCGCAACGAGGTCCACGGCGCCTCCCACAGCCCCGTTTGTG‑CCACCCCTCCCGGCAACTGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAt                                                                                      <  2:83901/143‑1 (MQ=255)                                                         aacaccccTCCCCCCAGCGAAATCCGAAACGAGGCCCACAGCGCCTGCCACAGCCCCGGTTGTG‑CCAGCCCTTCCGGCAACGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGcaca                                                                                        <  2:27467/141‑1 (MQ=255)                                                         aaTACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAtt                                                                                     >  2:628613/1‑150 (MQ=255)                                                          aTACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTa                                                                                    >  2:842859/1‑150 (MQ=255)                                                            acacCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAac                                                                                  >  1:1030176/1‑150 (MQ=255)                                                             ccccTGCCGCCAGCGCAATACGCAACGAGGCCCACAGCGCCTGCCCCAGCCCCGGTTGGG‑CCAGCACTCCCGGCAACGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAac                                                                                  <  2:525889/147‑1 (MQ=255)                                                                  gccgccAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCgtgt                                                                           >  2:48052/1‑151 (MQ=255)                                                                   ccgccAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTc                                                                          >  1:820325/1‑151 (MQ=255)                                                                    cgccAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCg                                                                         >  2:257263/1‑151 (MQ=255)                                                                          cgattaaCGCAACGGGGCCCACAGCGCCGGCCACAGCCCCGGTTGTG‑CCAGCCCTCCCGGAAACGGGCGATTTCCCCCAGCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCaa                                                                      <  2:943124/142‑1 (MQ=255)                                                                            cAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCg                                                                 >  2:145897/1‑151 (MQ=255)                                                                            cAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACAGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCg                                                                 >  1:518685/1‑151 (MQ=255)                                                                             agtcctccaCGCGGTCCACGGCGCCTCCCACACCGCCGGGTGTC‑CCAGCCCTCCCGGCACGGGGAGATTGCCCCCATCAACTTTCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGc                                                                  <  2:297509/141‑1 (MQ=255)                                                                             aaTCCGCAAGGAGGCCCACGCCCCCGCCCACAGCACCGGTTGTG‑CCAGCACTCCCGGAAAGTGGCGATTTCCCCCACCGACGATCGCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCg                                                                 <  1:1142158/150‑1 (MQ=255)                                                                                 cGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTAt                                                            >  1:1162641/1‑151 (MQ=255)                                                                                       gggccccACGGCCCCTGCCACAGCCCCGGTTGGG‑CCAGCACTTCCGGCAACGGGCGTTTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGAc                                                      <  2:178169/146‑1 (MQ=255)                                                                                       gAGTTCCCCAGCCCCTGCCACGCCCCGGGTGGGG‑CCAGCACTCCCGGCAACGGGCGATTTCCCCCAGCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGa                                                       <  1:378045/150‑1 (MQ=255)                                                                                       gAGGCCCACAGCGCCTGCCACAGCCCGGGTTGCG‑CCAGCACTCCCGGCAACGGGCGATTCCCCCCACCGCCGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCt                                                              <  1:584298/143‑1 (MQ=255)                                                                                         ggcccACAGCGCCTGCCACAGCCCGGGTTGGG‑CCAGCACTCCCGGCAACGGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCt                                                              <  1:140597/138‑1 (MQ=255)                                                                                            ccacggccccgcccACAGCCCCGGTTGGG‑CCAGCCCTCCCGGCAACGGGCGTTTTCCCCCACCGACAATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCg                                                 <  2:543745/139‑1 (MQ=255)                                                                                            ccACAGCCCTTGCCACCGCCCCGGTGGTC‑CCCGCCCCTCCGGCCGGGGGCATTTTCCCCCCTGGACGAACACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATc                                                  <  2:514960/150‑1 (MQ=255)                                                                                            caacccgccgtGCAACAGACCCGTTTCTG‑CCGGCCCTTCCGCCAGCTGGTATTTTCCCCGAACTAAGTTAACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGAc                                                      <  1:249058/136‑1 (MQ=255)                                                                                             cacaGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGt                                                >  2:899169/1‑151 (MQ=255)                                                                                                    ccTGCCACAGCCCCGGTTGGG‑CCAGCCCTTCCGGCAACTGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTc                                         <  2:147253/151‑1 (MQ=255)                                                                                                    ccTCCCCCAGCCCCGGTTGGC‑CCAGCCCTTCCGGCAGCTGGCTTTTTCCCCCACCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCg                                             <  1:888928/147‑1 (MQ=255)                                                                                                    ccTCCCACAGCCCCGGTTGCC‑CCACCACTCCCGGCAACGGGCGATTCCCCCCACCCACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGAt                                                   <  1:145886/141‑1 (MQ=255)                                                                                                     cTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTcg                                        >  1:54245/1‑151 (MQ=255)                                                                                                        ccccaGCCCCGGTTGCC‑CCACCCCTCCCGGAAACGGGCGATTTCCCCCATCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTcgc                                       <  2:772170/146‑1 (MQ=255)                                                                                                          acaGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTc                                         >  2:190182/1‑145 (MQ=255)                                                                                                                  ggTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAGCACCGTGTCGCAGATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTg                             >  1:435019/1‑149 (MQ=255)                                                                                                                   gTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgc                            >  2:374116/1‑149 (MQ=255)                                                                                                                      gtg‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgcagc                         >  2:425717/1‑149 (MQ=255)                                                                                                                       gtcccAGCCCTCCCGGCAACGGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGcagcag                       <  2:766262/148‑1 (MQ=255)                                                                                                                          ccAGCCCTCCCGGCACCTGGCGTTTTCCCCCATCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCg                     <  1:323183/150‑1 (MQ=255)                                                                                                                              ccctcccGGCAACGGGCGATTTCCCCCACCGACGCTCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtggtg                <  1:57358/146‑1 (MQ=255)                                                                                                                                  ccggacAACTGGCGATTTCCCCCACCGACAAACACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGCCGCGCCAGCCTTGCAGCAGCGTGGTGc               <  2:1096099/143‑1 (MQ=255)                                                                                                                                   ccGGCACCGGGCGATTCCCCCCACCGACAACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtggt                 <  1:769431/145‑1 (MQ=255)                                                                                                                                    cGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCggg           >  2:950171/1‑150 (MQ=255)                                                                                                                                       cAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgcgc        >  1:1014125/1‑150 (MQ=255)                                                                                                                                        aaCTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgcgc        >  1:38956/1‑149 (MQ=255)                                                                                                                                           tGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATgg   >  2:1020038/1‑151 (MQ=255)                                                                                                                                            ggCGTTTTCCCCCACCGACGACCCCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGc  <  1:788623/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               GCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTG‑CCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGC  >  NC_000913/73435‑73723

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: ‑