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%	115.0 / 11.4	65	V314G (GTA→GGA)	thiP	thiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (30/6);  new base C (0/29);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 6.53e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGG  >  NC_000913/73437‑73722                                                                                                                                                 |                                                                                                                                              cgcgCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCAt                                                                                                                                         >  1:609267/1‑151 (MQ=255)     cccGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATcgac                                                                                                                                     <  2:437389/151‑1 (MQ=255)          aGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATc                                                                                                                                 >  1:1609/1‑150 (MQ=255)             tCGCGACTGCTCCATATCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCg                                                                                                                             >  2:397778/1‑151 (MQ=255)                  aCTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGt                                                                                                                        >  1:948318/1‑151 (MQ=255)                   cTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGt                                                                                                                        >  2:627725/1‑150 (MQ=255)                    tGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTa                                                                                                                       >  1:506540/1‑150 (MQ=255)                        ccATAGCAGCATCAGGGCCAGCACTCCCCAAAATCCCCCGGCCGCCAGCGAAATACGCACCGGGGCCCACGGCGCCTCCCACAGCCCCGGTTGTGCCAGCCCTTCCGGCAACTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACg                                                                                                                    <  1:434865/149‑1 (MQ=255)                        ccATAACAGCATCATGGTCACCCCCCCGCAAAATCCCCCTGCCGCCGGCGCAATCCGCAACGGGGCCCACAGCGCCTGCCACAGCCCGGGTTGCGCCAGCACTTCCGGCAACTGGCTATTTCCCCCATCGACGATCACCGCCAGTAACg                                                                                                                    <  1:142684/149‑1 (MQ=255)                         cATACCAGCATCATGGCCAGCCCTCCCCAAAATCCCCCTCCCCCCAGCGCAATCCCCAACGAGGCCCACAGCGCCTCCCACAGCCCGGGTTGTCCCACCCCTCCCGCCAACGGGCGTTTTCCCCCACCGACGATCACCGCCAGTAAc                                                                                                                     <  1:170806/147‑1 (MQ=255)                           aaccagcaTCATGTTCAGCCCTCCCCACAATCCCCCTCCCCCCAGCACAATACGCAACGAGGCCCACAGCGCCTGCCACAGCCCCGGTTGGGCCAGCACTTCCGGGAACTGGCGATTTACCCCACCGACGATCACCGCCAGTAACGGTg                                                                                                                 <  1:596051/146‑1 (MQ=255)                               agcaTCAGGTCCACCCCTCCGCACAACACCCCGGCCGCCGGCGAAATCCGCAACGGGGCCCCCGGCGCCTGCCACAGCCCCGTTTGTCCCAGCCCTTCCGGCAACTGGCGATTTCCCCAACCGACGATCACCGCCAGTAACGGTGgcag                                                                                                             <  1:668834/149‑1 (MQ=255)                                 catcatGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAc                                                                                                         >  2:253612/1‑151 (MQ=255)                                 catcatGGTCAGCACTACGCAAAATCCCCCTGCCGCCAGCGCAATACGCAACGAGGCCCACAGCGCCTGCCACAGCCCCGGTTGGGCCAGCACTTCCGGCAACTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGgc                                                                                                               <  2:252081/145‑1 (MQ=255)                                  atcatGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAc                                                                                                         >  2:604138/1‑150 (MQ=255)                                    catGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAcagc                                                                                                      >  2:377142/1‑151 (MQ=255)                                       ggTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAcagcagc                                                                                                   >  2:623749/1‑151 (MQ=255)                                           accACTCCGCAAAATCCCCCCCCCGCCAGCCCAATCCGAAAGGAGCTCCACAGCCCCTGCCACGCCACGGGTGGCGCCAGCCCTTCGGGCAACTGGCGATTTCCCCCATCTACCATCACCGCCAGTAACGGTGGCAGCAACagcagcagc                                                                                                <  2:754417/148‑1 (MQ=255)                                             cACTCCGCACAATCCCCCTGCCGCCAGCGCAATCCGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTCCCAGCACTTCCGGCAACTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGcc                                                                                             <  2:496541/151‑1 (MQ=255)                                                  cccacaATCCCCCCGCCCCCAGCGCAATACGCAAGGGGGCCCACGGCGCCTCCCACACCCCCGGTTGCGCCAGCCCTCCCGGCAACGGGCGATTCCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGc                                                                                          <  1:851871/147‑1 (MQ=255)                                                  cccaaaatcccccTGCCCCCGCCGCAATCCGCAACGAGTCCCACAGCCCCTGCCACAGCCCCGGTTGTGCCAGCCCTTCCGGCAACGGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAg                                                                                           <  1:262417/137‑1 (MQ=255)                                                     acaATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAt                                                                                     >  2:441491/1‑151 (MQ=255)                                                      caATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAtt                                                                                    >  2:1020827/1‑151 (MQ=255)                                                        aTACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTa                                                                                   >  2:4859/1‑150 (MQ=255)                                                               tGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCgtg                                                                           >  1:796346/1‑151 (MQ=255)                                                                  cgccAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCg                                                                        >  2:282604/1‑151 (MQ=255)                                                                    ccAGCGCAATCCGCAACGAGGCCCACAGCGCCTGCCACAGCCCCGGTTGCGCCAGCCCTTCCGGCACCTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGc                                                                       <  2:268686/150‑1 (MQ=255)                                                                          cAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCg                                                                >  1:275485/1‑151 (MQ=255)                                                                           aacccGCAACGAGGCCCCCGGCGCCTGCCCCAGCCCCGGTTGTGCCAGCCCTTCCGGCAACTGGCGTTTTCCCCCATCAACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCgg                                                               <  2:165060/147‑1 (MQ=255)                                                                           aaaacaaaaCGAGGCCCACAGCGCCGGCCACAGCGCCGTTGGTGCCAGCTTTTCCGGCAACGGGAAATTTACCCCATCGCCGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGc                                                                       <  1:580514/136‑1 (MQ=255)                                                                            aTACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGc                                                              >  1:949122/1‑151 (MQ=255)                                                                               cGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTAt                                                           >  2:98832/1‑151 (MQ=255)                                                                               cGCAACGAGGCCCCCAGCCCCTGCCACAGCCCCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTa                                                            <  2:462921/150‑1 (MQ=255)                                                                                  acggAGGCCCCCGGCCCCTCCCACAGCCCCGGTTGTGCCAGCACTCCCGGCAGCGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTAt                                                           <  1:770394/145‑1 (MQ=255)                                                                                  aaCGAGGCCCACGGCGCCTGCCACACCCCCGGTTGGCCCAGCCCTCCCGGCACCGGGCGATTTCCCCCATCGCCGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTAt                                                           <  1:91760/148‑1 (MQ=255)                                                                                   aCGGGGCCCACCGCCCCTGCCACACCACCGGTTGTCCCAGCCTTCCCGTCAACTGGGGATTTCCCCCACCGCCGATCACCGCCAGTAACGGTTGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCa                                                        <  2:218147/150‑1 (MQ=255)                                                                                      aGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACg                                                    >  2:179638/1‑151 (MQ=255)                                                                                       ggcccACAGCGCCTGCCACAGCCCCGGTTGTGCCAGCCCTCCCGGAAACTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAg                                                       <  2:597664/144‑1 (MQ=255)                                                                                       ggcccACAGCGCCTGCCACAGCACCGGTTGTGCCAGCCCTTCCGGCAACTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAg                                                       <  2:535258/144‑1 (MQ=255)                                                                                          ccagccccccTGCACCAGCACTGTTTGTCCCAGCACTGCGGGCAGCTGGCGATTTACCCCACCGAGAATCTCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCg                                                <  2:801396/143‑1 (MQ=255)                                                                                          ccACAGCGCCTCCCACAGCCCCGGTTGTGCCAGCCCTCCCGGCAACTGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCg                                                <  1:253594/151‑1 (MQ=255)                                                                                           cacgGCCCCTGCCACAGCCCCGGTTGTGCCCGCCCTCCCGGCAACTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGt                                               <  1:949417/151‑1 (MQ=255)                                                                                                  ccTGCCCCAGCCCCGGTTGTGCCAGCCCTCCCGGCAACTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCg                                            <  1:129333/147‑1 (MQ=255)                                                                                                  ccTGCCACAGCACCGGTTGGGCCAGCACTCCCGGCAACGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTc                                        <  2:121410/151‑1 (MQ=255)                                                                                                     gCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTcgcg                                     >  2:1009545/1‑151 (MQ=255)                                                                                                      ccACAGCACCGGTTGTGCCAGCACTCCCGGCAACGGGCGATTTACCCCAGCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTcg                                       <  1:441457/148‑1 (MQ=255)                                                                                                       cacaGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGcc                                   >  2:61467/1‑151 (MQ=255)                                                                                                          acccccGGTTGTGCCAGCCCTCCCGGCAACTGGCGATTTCCCCCACCGCCGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAg                                 <  2:240008/146‑1 (MQ=255)                                                                                                           gCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGcc                               >  2:461119/1‑151 (MQ=255)                                                                                                              ccGGTTGTCCCAGCACTTCCGGCAACGGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCtt                             <  2:837752/150‑1 (MQ=255)                                                                                                                   tgggCCAGCCCTCCCGGCAACGGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgcagc                        <  1:377115/147‑1 (MQ=255)                                                                                                                        cAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtg                  >  1:374777/1‑151 (MQ=255)                                                                                                                         aGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtgg                 >  2:413287/1‑151 (MQ=255)                                                                                                                           ccctcccGGCACCGGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtggtg               <  1:407509/146‑1 (MQ=255)                                                                                                                            aCTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGc              >  2:205304/1‑151 (MQ=255)                                                                                                                             cTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGcc             >  1:3692/1‑151 (MQ=255)                                                                                                                                ccGGCAACTGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCgg           <  2:740043/150‑1 (MQ=255)                                                                                                                                ccGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCggg          >  2:73948/1‑151 (MQ=255)                                                                                                                                ccGGCAACGGGCGTTTTCCCCCATCGACAACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCggg          <  2:948412/151‑1 (MQ=255)                                                                                                                                  ggCACCGGGCGTTTCACCCCGGCGACGTTCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgcg        <  2:129301/151‑1 (MQ=255)                                                                                                                                   gCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgcgc       >  2:800683/1‑151 (MQ=255)                                                                                                                                    cAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgcgc       >  2:33860/1‑150 (MQ=255)                                                                                                                                    aatggggCGATTCCCCCCAGCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCa      <  1:282582/146‑1 (MQ=255)                                                                                                                                        tGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATgg  >  2:603116/1‑151 (MQ=255)                                                                                                                                         ggCGATTCCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAAt    <  2:838786/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                              CGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGG  >  NC_000913/73437‑73722

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