Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,674,177	G→T	12.3%	E89* (GAG→TAG) ‡	tqsA →	autoinducer 2 exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,674,177	0	G	T	12.3%	224.5 / 18.0	73	E89* (GAG→TAG) ‡	tqsA	autoinducer 2 exporter
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base T (4/5);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

TCGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTC  >  NC_000913/1674028‑1674323                                                                                                                                                      |                                                                                                                                                    tcGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACg                                                                                                                                                     <  1:344195/150‑1 (MQ=255)   gCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAg                                                                                                                                                   <  1:21821/150‑1 (MQ=255)    cATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGt                                                                                                                                                  <  1:818312/150‑1 (MQ=255)      ttatttatTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTg                                                                                                                                                <  2:608366/150‑1 (MQ=255)           tattGCTGTTATTCTTATCCCGCTGGTGCAACACATGGTCCGCTGTCGTGTGCCGCGGGTACTTGCGGTGTCGATTTTGATGACCATCATCGTTATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑cgcg                                                                                                                                          <  2:122418/150‑1 (MQ=255)           tattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑cgcg                                                                                                                                          <  2:1108160/150‑1 (MQ=255)           tattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGagcgcg                                                                                                                                          <  1:1155049/150‑6 (MQ=255)            attGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGAGCGCgg                                                                                                                                         <  2:954647/150‑1 (MQ=255)            attGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTACGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCgg                                                                                                                                         <  2:619703/150‑1 (MQ=255)             ttGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGAGCGCGGa                                                                                                                                        <  1:1259125/150‑1 (MQ=255)               gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACg                                                                                                                                      >  1:1396821/1‑150 (MQ=255)               gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGAGCGCGGaac                                                                                                                                      <  1:88702/150‑3 (MQ=255)                cTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGt                                                                                                                                     >  1:844198/1‑150 (MQ=255)                    tATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTAcc                                                                                                                                 <  1:37878/150‑1 (MQ=255)                            aCCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATc                                                                                                                         >  1:49502/1‑150 (MQ=255)                                 cTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAAc                                                                                                                    >  1:1125584/1‑150 (MQ=255)                                   ggTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAActc                                                                                                                  >  2:1043415/1‑150 (MQ=255)                                                 tCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgccgc                                                                                                    <  1:1250521/150‑1 (MQ=255)                                                   cGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTGGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTg                                                                                                  <  2:1003711/150‑1 (MQ=255)                                                    gCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATCGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGc                                                                                                 <  2:118778/150‑1 (MQ=255)                                                     cTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCa                                                                                                <  1:127949/150‑1 (MQ=255)                                                        gCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAgc                                                                                             >  1:1232386/1‑150 (MQ=255)                                                        gCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAgc                                                                                             <  2:612638/150‑1 (MQ=255)                                                           tgtgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                          >  1:197233/1‑150 (MQ=255)                                                             tgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTg                                                                                        <  2:48601/150‑1 (MQ=255)                                                             tgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTg                                                                                        <  2:548243/150‑1 (MQ=255)                                                              gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGa                                                                                       >  1:796396/1‑150 (MQ=255)                                                                   gtgtACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTCCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgt                                                                                  <  1:260124/150‑1 (MQ=255)                                                                    tgtACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgtt                                                                                 >  1:523285/1‑150 (MQ=255)                                                                        cTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGttgttg                                                                             <  2:908427/150‑1 (MQ=255)                                                                        cTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGttgttg                                                                             <  1:1287868/150‑1 (MQ=255)                                                                          ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                           >  1:945232/1‑150 (MQ=255)                                                                          ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                           >  1:787868/1‑150 (MQ=255)                                                                          ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                           >  1:1366074/1‑150 (MQ=255)                                                                             ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCaa                                                                                               >  2:1236875/1‑127 (MQ=255)                                                                             ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCaa                                                                                               <  1:1241953/127‑1 (MQ=255)                                                                             ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCaa                                                                                               <  1:1236875/127‑1 (MQ=255)                                                                             ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCaa                                                                                               >  2:1241953/1‑127 (MQ=255)                                                                             ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAAcg                                                                        >  2:1044177/1‑150 (MQ=255)                                                                                gtCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGt                                                                     >  1:907329/1‑150 (MQ=255)                                                                                   gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAgg                                                                  >  1:562051/1‑150 (MQ=255)                                                                                    aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                 >  2:969809/1‑150 (MQ=255)                                                                                         gatgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCGAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGa                                                            >  2:1104849/1‑150 (MQ=255)                                                                                           tgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACg                                                          <  1:605496/150‑1 (MQ=255)                                                                                             aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtc                                                        <  2:796396/150‑1 (MQ=255)                                                                                              ccATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtct                                                       >  1:371293/1‑150 (MQ=255)                                                                                                  catcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                   >  1:665358/1‑150 (MQ=255)                                                                                                  catcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                   <  2:785513/150‑1 (MQ=255)                                                                                                     cGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGa                                                <  2:1396821/150‑1 (MQ=255)                                                                                                       tGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGAcc                                              >  2:1098360/1‑150 (MQ=255)                                                                                                          tGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGc                                           <  1:969809/150‑1 (MQ=255)                                                                                                             cGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTgg                                        >  2:682215/1‑150 (MQ=255)                                                                                                             cGATGGTGTTGCTATTAGCTTATCTCGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTAACGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGAAGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCtcc                                        >  2:159075/1‑148 (MQ=255)                                                                                                                        cTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATAtt                             >  1:143829/1‑150 (MQ=255)                                                                                                                         tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTg                            <  2:406054/150‑1 (MQ=255)                                                                                                                         tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTg                            <  2:401075/150‑1 (MQ=255)                                                                                                                            tAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATc                         >  1:887336/1‑150 (MQ=255)                                                                                                                             aGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATcc                        >  2:182633/1‑150 (MQ=255)                                                                                                                                ttATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGaa                     >  2:340833/1‑150 (MQ=255)                                                                                                                                  aTCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcg                   <  2:1125584/150‑1 (MQ=255)                                                                                                                                      gggTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcggc               <  1:682215/150‑1 (MQ=255)                                                                                                                                        gTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGa             <  2:1232386/150‑1 (MQ=255)                                                                                                                                        cgtCCGCG‑TCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGa             >  1:1034814/3‑150 (MQ=255)                                                                                                                                         ttCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt            <  1:1043415/150‑1 (MQ=255)                                                                                                                                         gtCCGCG‑TCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt            >  1:1372547/2‑150 (MQ=255)                                                                                                                                            cgcgCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGAc         >  1:1268853/1‑150 (MQ=255)                                                                                                                                            ccgcgtCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcggc               <  1:1432785/140‑1 (MQ=255)                                                                                                                                            ccgcgtCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcggc               >  2:1432785/6‑145 (MQ=255)                                                                                                                                             gcgcTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACg        >  1:1021073/1‑150 (MQ=255)                                                                                                                                                 tCAACTCGTTGAGCGCGGAAGTTACAGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGAAACGCGTAGGGATTGACGTCTCAGTTTACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTg     >  2:96786/1‑150 (MQ=255)                                                                                                                                                  cAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt   >  1:1184440/1‑150 (MQ=255)                                                                                                                                                   aaCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTc  >  2:1265172/1‑150 (MQ=255)                                                                                                                                                   aaCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTc  >  1:164422/1‑150 (MQ=255)                                                                                                                                                   aaCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTc  <  2:787868/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    TCGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTC  >  NC_000913/1674028‑1674323

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