Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,697,374	T→G	12.3%	intergenic (+28/‑49)	hslV → / → hslU	ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,697,374	0	T	G	12.3%	134.5 / 9.6	66	intergenic (+28/‑49)	hslV/hslU	ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base T (27/24);  new base G (2/6);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGC  >  NC_002947/5697225‑5697521                                                                                                                                                       |                                                                                                                                                    tgctgAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGt                                                                                                                                                     >  2:1175806/1‑150 (MQ=255)      aaCAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAgg                                                                                                                                                <  1:685999/150‑1 (MQ=255)         aaGACACACCTGTCGGCCCGCGAAATTGCAGAGATCGCCCTGAACATCGCCGGCGACATATGCGTGTTCACCAACCACACCCTGCCCTTCGAGGAGGAGGACCTGGCCGACTGATTCGGCTGTTCTGGCGCCCCTTC‑ACCGGGAAGGGGc                                                                                                                                             <  2:26126/150‑1 (MQ=255)          ggaccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑ACCGGGAAGGGGcg                                                                                                                                            <  2:6328/149‑2 (MQ=255)           gaccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCtt                                                                                                                                           >  2:802212/1‑150 (MQ=255)           gaccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑ACCGGGAAGGGGcgc                                                                                                                                           <  1:175101/150‑3 (MQ=255)           gaccgaccTCTCGGCCCGCGAAATCGCCGAGACAGCCCTGAAAATCGCCGGCGACATCTGCGTGTCCCCCAACCACAACCTGACCATCGAGGATCCGAACCTGGCCGCCTGAACCACCCGTTCACGCCCTCCTTC‑CCGGTGAAGGGGcgg                                                                                                                                           >  2:561315/1‑148 (MQ=255)            accgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCttt                                                                                                                                          <  1:1175806/150‑1 (MQ=255)             ccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCtttt                                                                                                                                         >  1:277417/1‑150 (MQ=255)                 ccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCAc                                                                                                                                     <  1:802212/150‑1 (MQ=255)                   tctcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGc                                                                                                                                   >  1:1065898/1‑150 (MQ=255)                     tcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGc                                                                                                                                   >  2:472892/1‑148 (MQ=255)                     tcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGc                                                                                                                                   <  1:472892/148‑1 (MQ=255)                         cccGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTc                                                                                                                             >  2:791083/1‑150 (MQ=255)                          ccGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCa                                                                                                                            >  2:909436/1‑150 (MQ=255)                                   cGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCaa                                                                                                                   <  2:176866/150‑1 (MQ=255)                                      cGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGCCTGATTCAGCTGTTCTGGCTCCGCAAC‑CGCGTTAAGTTGCTTTTCCACGCTGACTCACCTGCCCAAGGa                                                                                                                <  1:1114485/150‑1 (MQ=255)                                             gCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGG‑GAAGGGGc                                                                                                                                             <  2:928100/114‑1 (MQ=255)                                             gCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATACAGCTGTTCTGGCGCCCCTTCACCGG‑GAAGGGGc                                                                                                                                             >  1:928100/1‑114 (MQ=255)                                               ccTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGa                                                                                                       <  1:676705/150‑1 (MQ=255)                                               ccTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGa                                                                                                       >  2:630876/1‑150 (MQ=255)                                                cTGAACCTCGCCGGCGCTAGCTTCGTGTGCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGAt                                                                                                      <  2:374490/150‑1 (MQ=255)                                                      aTCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTc                                                                                                >  1:28270/1‑150 (MQ=255)                                                         gCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCAt                                                                                             <  2:277417/150‑1 (MQ=255)                                                         gCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCAt                                                                                             <  1:909436/150‑1 (MQ=255)                                                         gCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATATCCAt                                                                                             <  2:340265/150‑1 (MQ=255)                                                           cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGa                                                                                           <  1:974182/150‑1 (MQ=255)                                                                   tCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc                                                                                   <  2:791516/150‑1 (MQ=255)                                                                      gCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCgag                                                                                <  1:791083/150‑1 (MQ=255)                                                                              accaacaaCACCCTGACCATCTAGGAGAAGGACCTGGCCAACTGGTTCAGGTGTTCTGGCGCCCCTTC‑CCGGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTcaa                                                                        >  2:1200201/1‑148 (MQ=255)                                                                               ccaaccaCAACTTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGAGCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCAc                                                                       >  2:438957/1‑150 (MQ=255)                                                                               ccaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAg                                                                                                                  >  2:247142/1‑107 (MQ=255)                                                                               ccaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAg                                                                                                                  <  1:247142/107‑1 (MQ=255)                                                                                caaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑ACCGGGAAGGGGc                                                                                                                                             <  2:847328/79‑1 (MQ=39)                                                                                caaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑ACCGGGAAGGGGc                                                                                                                                             >  1:847328/1‑79 (MQ=39)                                                                                  accaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTCTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGCTAATGTCCCTAACCGCCCGCAATAGTGTCCAGGaa                                                                    >  2:1163697/1‑150 (MQ=255)                                                                                    cacaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGGTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACt                                                                  >  2:1137134/1‑150 (MQ=255)                                                                                     acaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGAcccccc                                                                                     <  2:538210/130‑1 (MQ=255)                                                                                     acaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGAcccccc                                                                                     >  1:538210/1‑130 (MQ=255)                                                                                        aCCTGACCATCGAGGAGCAGGACCTGACCGACTGATTCAGCTGTTCTGGCGAACCTTC‑CCGGTGAAGGGGCTTTTCGACGCTGACTCACCTGCCCAACGCCCGTATTAATCATCTCCATGACCCCCCGCGACATCCTCCACGGACTCAAc                                                              >  2:295457/1‑150 (MQ=255)                                                                                               cATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCccca                                                       <  1:1137134/150‑1 (MQ=255)                                                                                               cATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTACTGGCACCCATTC‑CAGGAAAAGCGACACTAAAAAGCCGACTAACTATCACAATAACAGTATGGATGAGGGCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCcaca                                                       >  2:439369/1‑150 (MQ=255)                                                                                               cATCGAGGAGCAGGACCCGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCcaca                                                       >  2:619348/1‑150 (MQ=255)                                                                                                 ccGAGGAGCACGACCTCGCCGACTGATTCACCCCTTCTGGCGCCCCTTC‑CCGGGGAAGGGGCTTTTCCACGCTTACTCACCTGCCCAAGGACCGTTTTGTTTTTATCCAACACCCCCCACCGAAACCTCCACGAACTCAACCGCCAcatc                                                     <  2:388732/149‑1 (MQ=255)                                                                                                    aggagAAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACatcatc                                                  <  2:388924/150‑1 (MQ=255)                                                                                                         cAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCa                                             <  2:845176/150‑1 (MQ=255)                                                                                                            gACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGGATTGCTCATGTCCACGACCCCCCGCGAGATCGTCCACTAACTCAACCGCCACATCATCGGCCAggc                                          >  2:436249/1‑149 (MQ=255)                                                                                                                        tGATTCAGATGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcg                              >  2:1201055/1‑150 (MQ=255)                                                                                                                            tCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGAACATGTCCATGACCCCCCGCGAGCTCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGcc                          >  1:640761/1‑150 (MQ=255)                                                                                                                                 tGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCa                                                                                              >  1:787260/1‑77 (MQ=255)                                                                                                                                 tGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCa                                                                                              <  2:787260/77‑1 (MQ=255)                                                                                                                                 tGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGt                                                                           >  1:233994/1‑96 (MQ=255)                                                                                                                                 tGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGt                                                                           <  2:233994/96‑1 (MQ=255)                                                                                                                                       ggCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGc               <  1:903780/150‑1 (MQ=255)                                                                                                                                        gcgcCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGCGGCCATTGcc              >  2:662378/1‑150 (MQ=255)                                                                                                                                        gcgcCCCTTC‑CCGGTGAAGGGGCTGTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTTCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGcc              >  1:226712/1‑150 (MQ=255)                                                                                                                                          gcCCCTTCACCGG‑GAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCaa                                 >  1:1198974/1‑129 (MQ=255)                                                                                                                                          gcCCCTTCACCGG‑GAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCaa                                 <  2:1198974/129‑1 (MQ=255)                                                                                                                                          gcCCCTTCACCGG‑GAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTg                       >  1:890493/1‑139 (MQ=255)                                                                                                                                          gcCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCt            >  1:45589/1‑150 (MQ=255)                                                                                                                                          gcCCCTTC‑ACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTg                       <  2:890493/139‑1 (MQ=255)                                                                                                                                           ccccTTCACCGG‑GAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTg           >  1:183836/1‑150 (MQ=255)                                                                                                                                               ttC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCCCCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCa       >  1:1189133/1‑150 (MQ=255)                                                                                                                                                 c‑ccGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAAc     >  1:959717/1‑150 (MQ=255)                                                                                                                                                     ggTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGc  <  1:438957/150‑1 (MQ=255)                                                                                                                                                      gTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCa                                             <  1:922998/106‑1 (MQ=255)                                                                                                                                                      gTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCa                                             >  2:922998/1‑106 (MQ=255)                                                                                                                                                       |                                                                                                                                                    TGCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACCGC  >  NC_002947/5697225‑5697521

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