Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,674,178	A→C	12.2%	E89A (GAG→GCG) ‡	tqsA →	autoinducer 2 exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,674,178	0	A	C	12.2%	199.0 / 18.5	73	E89A (GAG→GCG) ‡	tqsA	autoinducer 2 exporter
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base C (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 = 1.00e+00

GCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTC  >  NC_000913/1674030‑1674323                                                                                                                                                     |                                                                                                                                                   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                                                                             <  1:1287868/150‑1 (MQ=255)                                                                      cTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGttgttg                                                                             <  2:908427/150‑1 (MQ=255)                                                                        ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                           >  1:1366074/1‑150 (MQ=255)                                                                        ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                           >  1:787868/1‑150 (MQ=255)                                                                        ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                           >  1:945232/1‑150 (MQ=255)                                                                           ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCaa                                                                                               >  2:1241953/1‑127 (MQ=255)                                                                           ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCaa                                                                                               <  1:1241953/127‑1 (MQ=255)                                                                           ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCaa                                                                                               <  1:1236875/127‑1 (MQ=255)                                                                           ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCaa                                                                                               >  2:1236875/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)                                                                                                                                                     |                                                                                                                                                   GCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTC  >  NC_000913/1674030‑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: ‑