Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,175,797	(C)6→7	coding (210/582 nt)	yqiA ←	esterase YqiA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,175,791	1	.	C	98.3%	202.7 / ‑3.3	58	coding (216/582 nt)	yqiA	esterase YqiA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (32/25);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

ACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TAT‑‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTC  >  NC_000913/3175652‑3175925                                                                                                                                              |                                                                                                                                         aCTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAa‑tat                                                                                                                                           >  2:857366/1‑140 (MQ=255)  ctctAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TAT‑‑c                                                                                                                                        <  1:1094311/140‑1 (MQ=255)  ctctAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TAT‑‑c                                                                                                                                        >  1:877208/1‑140 (MQ=255)      aGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TAT‑‑ccc‑cc                                                                                                                                   >  1:282201/1‑140 (MQ=255)             atTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGa                                                                                                                             >  1:899682/1‑140 (MQ=255)                 ctgcCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGaa                                                                                                                         >  2:965306/1‑140 (MQ=255)                 ctgcCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGaa                                                                                                                         >  2:6917/1‑140 (MQ=255)                  tgcCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAAc                                                                                                                        <  1:1055833/140‑1 (MQ=255)                    cccGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCg                                                                                                                      >  2:1003309/1‑140 (MQ=255)                    cccGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCg                                                                                                                      <  1:935667/140‑1 (MQ=255)                    cccGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TAT‑‑ccc‑ccc                                                                                                                                  <  2:356758/127‑1 (MQ=255)                    cccGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TAT‑‑ccc‑ccc                                                                                                                                  >  1:356758/1‑127 (MQ=255)                     ccGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGa                                                                                                                     <  2:270321/140‑1 (MQ=255)                     ccGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TAT‑‑ccc‑ccccagt                                                                                                                              >  2:45207/1‑126 (MQ=255)                     ccGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TAT‑‑ccc‑ccccagt                                                                                                                              <  1:45207/130‑5 (MQ=255)                       ggTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACa                                                                                                                   <  2:98829/140‑1 (MQ=255)                        gtgtAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACaa                                                                                                                  <  1:524937/140‑1 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  1:1014738/1‑140 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  2:52914/1‑140 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  1:1020952/1‑140 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  1:1108079/1‑140 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  2:536666/1‑140 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  2:1095530/1‑140 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  1:809314/1‑140 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  1:1168715/1‑140 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  1:548265/1‑140 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  2:540484/1‑140 (MQ=255)                                   ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGaa                                                                                                       >  2:1103872/1‑140 (MQ=255)                                          tGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATCAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGc                                                                                                <  2:845148/140‑1 (MQ=255)                                           gACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCTTAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGcc                                                                                               <  1:155966/140‑1 (MQ=255)                                           gACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGcc                                                                                               <  2:282201/140‑1 (MQ=255)                                           gACCGAGATAGCCCGTCAGCAGTTCAAACGGGCGCCCCGCCGGGTTTACCCCCCCTGCGGGCAGCAAAAAACATTGTGACAACCAGGGGGCGTAA‑TATCCCCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGc                                                                                                <  2:898874/140‑1 (MQ=255)                                            aCCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCa                                                                                              <  1:965306/140‑1 (MQ=255)                                            aCCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCa                                                                                              <  1:540484/140‑1 (MQ=255)                                                           cagcagTTCAAACGGGCGTACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATg                                                                               >  1:398439/1‑140 (MQ=255)                                                                gTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTc                                                                          >  2:1038593/1‑140 (MQ=255)                                                                   cAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCcag                                                                       <  2:1020952/140‑1 (MQ=255)                                                                           cgcACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTg                                                               >  2:562812/1‑140 (MQ=255)                                                                                cgccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCcgcg                                                          >  1:43990/1‑140 (MQ=255)                                                                                  ccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTc                                                        >  2:895331/1‑140 (MQ=255)                                                                                               cacTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcggc                                           >  1:257716/1‑140 (MQ=255)                                                                                                acTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGt                                                         >  1:37231/1‑125 (MQ=255)                                                                                                acTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGt                                                         <  2:37231/125‑1 (MQ=255)                                                                                                 cTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcggcgg                                         <  1:417236/140‑1 (MQ=255)                                                                                                    cGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCaa                                      >  2:1111359/1‑140 (MQ=255)                                                                                                     gggCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAAc                                     >  1:1207645/1‑140 (MQ=255)                                                                                                      ggCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACt                                    >  2:750482/1‑140 (MQ=255)                                                                                                           caTAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGa                               <  1:1103872/140‑1 (MQ=255)                                                                                                            aTAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGa                               <  2:665808/139‑1 (MQ=255)                                                                                                            aTAAAACATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGa                               >  1:665808/1‑139 (MQ=255)                                                                                                                aagaTTGAGAAACCCGGGTGGCGTAATTCC‑‑CCC‑CCCAATGCCGAACCGACAAACCCCGACGAACCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGat                          <  2:676286/137‑1 (MQ=255)                                                                                                                 aCATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGatc                         <  2:1168715/140‑1 (MQ=255)                                                                                                                 aCATTGTGACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGatc                         <  1:410184/140‑1 (MQ=255)                                                                                                                     tgtgACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtct                     >  2:330032/1‑140 (MQ=255)                                                                                                                     tgtgACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtct                     >  2:18515/1‑140 (MQ=255)                                                                                                                     tgtgACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtct                     <  1:330032/140‑1 (MQ=255)                                                                                                                     tgtgACAACCAGGTGGCGTAA‑TAT‑‑CCCACCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtct                     >  1:149624/1‑140 (MQ=255)                                                                                                                        gACAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCaa                  >  1:1031818/1‑140 (MQ=255)                                                                                                                          cAACCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACg                <  1:807835/140‑1 (MQ=255)                                                                                                                           aaCCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGt               <  1:562812/140‑1 (MQ=255)                                                                                                                           aaCCAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGt               <  1:52914/140‑1 (MQ=255)                                                                                                                             ccAGGTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCa             >  1:1051077/1‑140 (MQ=255)                                                                                                                                 gTGGCGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGt         <  2:257716/140‑1 (MQ=255)                                                                                                                                     cGTAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGGGATg     >  2:198324/1‑140 (MQ=255)                                                                                                                                       tAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCaa                  >  2:393225/1‑125 (MQ=255)                                                                                                                                       tAA‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCaa                  <  1:393225/125‑1 (MQ=255)                                                                                                                                        aa‑TATC‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTc  <  2:809314/140‑1 (MQ=255)                                                                                                                                              |                                                                                                                                         ACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAA‑TAT‑‑CCC‑CCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTC  >  NC_000913/3175652‑3175925

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑