Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,886,825	C→G	6.0%	A83P (GCA→CCA)	PP_5161 ←	lysophospholipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,886,825	0	C	G	6.0%	300.3 / 6.6	94	A83P (GCA→CCA)	PP_5161	lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (44/44);  new base G (4/2);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGC  >  NC_002947/5886678‑5886971                                                                                                                                                    |                                                                                                                                                    tGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgcat                                                                                                                                                   >  2:1177866/1‑148 (MQ=255)              gCCGTGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTc                                                                                                                                      >  2:2154979/1‑150 (MQ=255)               ccGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCg                                                                                                                                     >  1:2312130/1‑150 (MQ=255)                 gggCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGaa                                                                                                                                   <  2:1356651/150‑1 (MQ=255)                   gCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAAcc                                                                                                                                 <  1:14612/150‑1 (MQ=255)                     aGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAAcccc                                                                                                                               >  2:853326/1‑150 (MQ=255)                        tCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCt                                                                                                                            >  2:924291/1‑150 (MQ=255)                              cTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGa                                                                                                                      <  1:924291/150‑1 (MQ=255)                              cTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGa                                                                                                                      <  2:2312130/150‑1 (MQ=255)                                gATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGaaa                                                                                                                    <  2:1204808/150‑1 (MQ=255)                                gATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGaaa                                                                                                                    <  1:1320656/150‑1 (MQ=255)                                  tCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGGTTGGCTACTCGGGCAGCCATGGCTTCCCGACCAGGTCGGACCCCGCTGCCTGAAAGc                                                                                                                  >  2:266955/1‑150 (MQ=255)                                    aCGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTa                                                                                                                >  1:1651969/1‑150 (MQ=255)                                      ggggAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTAcc                                                                                                              >  1:2403065/4‑150 (MQ=255)                                      ggCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCt                                                                                                                        >  2:1803735/1‑140 (MQ=255)                                      ggCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCt                                                                                                                        <  1:1803735/140‑1 (MQ=255)                                          aaaCCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACccagc                                                                                                          <  1:2154979/150‑1 (MQ=255)                                              ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGGTCGGGCAGGCATGGCTGCCCGagca                                                                                                                                          >  2:1885682/1‑111 (MQ=39)                                              ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGGTCGGGCAGGCATGGCTGCCCGAGCa                                                                                                                                          <  1:1885682/114‑1 (MQ=39)                                              ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCcagccagc                                                                                                      >  1:1858002/1‑150 (MQ=255)                                              ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCcagccagc                                                                                                      >  2:176339/1‑150 (MQ=255)                                                 gCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                   <  1:722679/150‑1 (MQ=255)                                                     gAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                   <  2:1864955/146‑1 (MQ=255)                                                     gAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                   >  1:1864955/1‑146 (MQ=255)                                                           ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGGTCGGGCAGGCATGGCTGCCCGagca                                                                                                                                          >  1:451292/1‑98 (MQ=38)                                                           ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGGTCGGGCAGGCATGGCTGCCCGAGCa                                                                                                                                          <  2:451292/101‑1 (MQ=38)                                                           ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAg                                                                                        <  2:1545435/150‑1 (MQ=255)                                                             aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGt                                                                                      >  1:515480/1‑150 (MQ=255)                                                             aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGt                                                                                      >  1:260509/1‑150 (MQ=255)                                                              cTCGATCACATGTCGGTAGAGGCCGATGTGGTCGTAATAGCCGTGCAGCAGGAACAGGGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCTGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CTGGtt                                                                                     >  1:337621/1‑150 (MQ=255)                                                                 gATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAg                                                                                  <  2:1548143/150‑1 (MQ=255)                                                                     acaTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCa                                                                              <  1:1996210/150‑1 (MQ=255)                                                                       aTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGaaa                                                                                                                    >  1:1981795/1‑111 (MQ=255)                                                                       aTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGaaa                                                                                                                    <  2:1981795/111‑1 (MQ=255)                                                                       aTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAg                                                                                  >  1:2115021/1‑144 (MQ=255)                                                                       aTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAg                                                                                  <  2:2115021/144‑1 (MQ=255)                                                                       aTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTAGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATa                                                                            >  1:667208/1‑150 (MQ=255)                                                                          tCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGaa                                                                         <  2:745914/150‑1 (MQ=255)                                                                            ggTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTc                                                                                                                                      <  1:1383042/88‑1 (MQ=255)                                                                            ggTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTc                                                                                                                                      >  2:1383042/1‑88 (MQ=255)                                                                             gTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAgcg                                                                      >  2:1104573/1‑150 (MQ=255)                                                                             gTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAgcg                                                                      >  1:2374473/1‑150 (MQ=255)                                                                                gagGCCCATGTAGTCATAATAGCCGTACAGCAGGAACAGCGTCGCGCTACGCGGCTCAGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTCCATGAANGNTNCAAAACCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAAAAGCGCTg                                                                   <  1:1507663/150‑1 (MQ=255)                                                                                   gCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAgg                                                                                 >  1:2310828/1‑133 (MQ=255)                                                                                   gCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAgg                                                                                 <  2:2310828/133‑1 (MQ=255)                                                                                   gCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTa                                                                >  1:1156781/1‑150 (MQ=255)                                                                                     ccATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAgcgcaggacta                                                              >  2:1435423/1‑143 (MQ=255)                                                                                            gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAg                                                                                  <  2:1215225/123‑1 (MQ=255)                                                                                            gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAg                                                                                  >  1:1215225/1‑123 (MQ=255)                                                                                            gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTgcgc                                                       >  2:473607/1‑150 (MQ=255)                                                                                             tCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAg                                                                                        >  1:883544/1‑116 (MQ=255)                                                                                             tCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAg                                                                                        <  2:883544/116‑1 (MQ=255)                                                                                             tCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCt                                                      >  1:79197/1‑150 (MQ=255)                                                                                               gTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtg                                                    >  2:1573025/1‑150 (MQ=255)                                                                                                 aaTAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtgtg                                                  >  2:1709042/1‑150 (MQ=255)                                                                                                    aGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCg                                               >  2:2112125/1‑150 (MQ=255)                                                                                                       cGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACa                                            <  2:2198696/150‑1 (MQ=255)                                                                                                       cGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACa                                            <  1:473607/150‑1 (MQ=255)                                                                                                        gTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCtt                                                     <  2:1380722/140‑1 (MQ=255)                                                                                                        gTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCtt                                                     >  1:1380722/1‑140 (MQ=255)                                                                                                             gcagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggct                                      >  1:407522/1‑150 (MQ=255)                                                                                                             gcagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggct                                      >  2:1483670/1‑150 (MQ=255)                                                                                                              cagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCtt                                                     <  2:1024987/134‑1 (MQ=255)                                                                                                              cagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCtt                                                     >  1:1024987/1‑134 (MQ=255)                                                                                                               agGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctgg                                    <  1:2112125/150‑1 (MQ=255)                                                                                                                 gAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCt                                                                    >  1:869426/1‑116 (MQ=255)                                                                                                                 gAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCt                                                                    <  2:869426/116‑1 (MQ=255)                                                                                                                   aCAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAg                                >  2:1396396/1‑150 (MQ=255)                                                                                                                     aGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtg                              <  2:1156781/150‑1 (MQ=255)                                                                                                                     aGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtg                              <  2:515480/150‑1 (MQ=255)                                                                                                                     aGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtc                              <  2:1651969/150‑2 (MQ=255)                                                                                                                         tCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGa                          >  2:663293/1‑150 (MQ=255)                                                                                                                            gttcgacctggtCGGGCAGGCATGGCTGCCCGAGCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACg                       >  1:2106338/8‑150 (MQ=39)                                                                                                                               ttGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgcg                    >  1:2312073/1‑150 (MQ=255)                                                                                                                                     gctCGGGCAGCCATGCCTGCCCGACCAGGTCGCACCCCGCTGAGTGAAAGCTATCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGc              <  1:1966572/150‑1 (MQ=255)                                                                                                                                       tCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTgg                                                                  <  1:393983/96‑1 (MQ=255)                                                                                                                                       tCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTgg                                                                  >  2:393983/1‑96 (MQ=255)                                                                                                                                       tCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg            <  1:1104573/150‑1 (MQ=255)                                                                                                                                       tCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg            <  1:1483670/150‑1 (MQ=255)                                                                                                                                       tCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg            <  1:1573025/150‑1 (MQ=255)                                                                                                                                        cGGGCCGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGCACACTACCCAGCGTG‑TGGCTGGCGGTCATGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCggg           <  1:285384/150‑1 (MQ=255)                                                                                                                                        cGGGCAGCCATGCCTGCCCGACAAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCggg           <  1:1177866/150‑1 (MQ=255)                                                                                                                                         gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGt          <  2:407522/150‑1 (MQ=255)                                                                                                                                         gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGt          <  1:852674/150‑1 (MQ=255)                                                                                                                                           gCAGCCATGCCTGCCCGAGCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgc        >  2:1362003/1‑150 (MQ=255)                                                                                                                                           gCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgg        >  1:1609749/1‑149 (MQ=255)                                                                                                                                             tctcATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgcgc      <  1:1477301/147‑1 (MQ=255)                                                                                                                                               ccATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAg    <  2:2403065/150‑1 (MQ=255)                                                                                                                                               ccATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAg    <  1:266955/150‑1 (MQ=255)                                                                                                                                               ccATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAg    <  2:1858002/150‑1 (MQ=255)                                                                                                                                               ccATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAg    >  1:1370204/1‑150 (MQ=255)                                                                                                                                                cATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAgg   >  1:1521258/1‑150 (MQ=255)                                                                                                                                                 aTGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCCTAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGGGGGTGCGCAGGc  >  1:2305595/1‑150 (MQ=255)                                                                                                                                                 aTGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGc  <  2:2305595/150‑1 (MQ=255)                                                                                                                                                 aTGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGACTCCAAACCTAAAAGCGCCGGTAGTTCAGCGCGTGTTCTGACATTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCaaga  >  2:1477301/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGG‑CAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGC  >  NC_002947/5886678‑5886971

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