Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,888,479	C→A	9.9%	A513D (GCT→GAT)	gcl →	glyoxylate carboligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,888,479	0	C	A	9.9%	188.8 / 8.0	64	A513D (GCT→GAT)	gcl	glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base A (3/4);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTACGTCCACGTACTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTGATTCTCGAGCGTGTGACCAACATTTCC  >  NC_002947/4888334‑4888627                                                                                                                                                   |                                                                                                                                                     ccGTACGTCCACGTACTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCtg                                                                                                                                                    <  3:228709/148‑1 (MQ=255) ccGTACGTCCACGTACTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCtg                                                                                                                                                    >  4:228709/1‑148 (MQ=255) ccGTACGTCCACGTACTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCt                                                                                                                                                     <  2:225112/147‑1 (MQ=255)  cGTACGTCCACTTAGTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACAAGGATAACTGTGTACAACTGGAGTTCGAGAACATCAACTC‑TACAGAAGCTGCCACATACGGTTTCGACCACGTCGAtgtg                                                                                                                                                  >  8:317955/1‑149 (MQ=255)    tACGTCCACGTACTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGt                                                                                                                                                >  1:174373/1‑149 (MQ=255)     aCGTCCACGTACTGGTGAACAACGCCTACCTGGGCCGGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACAGCGACGTGGTc                                                                                                                                               <  5:223912/149‑1 (MQ=255)      cGTCCACGTACTGGTGAACAACGCCTAACTGGGCATGATCCGCCAGGCGCAGCTTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACAGCGACGTGGTc                                                                                                                                               <  5:218881/148‑1 (MQ=255)        tCCACGTACTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGAAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACAGCGACGTGGTCGac                                                                                                                                            <  3:43011/149‑2 (MQ=255)         ccACGTACTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGAGTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAg                                                                                                                                           <  1:180551/149‑1 (MQ=255)          cACGTACTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAgg                                                                                                                                          <  2:174373/149‑1 (MQ=255)             gTACTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCt                                                                                                                                       >  2:207530/1‑149 (MQ=255)                cTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGGCGACCACGGCGCTGTGGTCGAAGGCCTgg                                                                                                                                     <  3:1856/148‑1 (MQ=255)                cTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACATACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTggc                                                                                                                                    >  7:107408/1‑148 (MQ=255)                 tGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGt                                                                                                                                   >  5:190785/1‑149 (MQ=255)                   gTGCACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCt                                                                                                                                       <  5:46181/143‑1 (MQ=255)                   gTGCACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCt                                                                                                                                       >  6:46181/1‑143 (MQ=255)                   gTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGtt                                                                                                                                  <  5:37334/148‑1 (MQ=255)                    tGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTAAGAGAACATCAACTC‑GACAGCCGCTGCCACCTACGGTGTCGACAACGTCGCTGTGGTCGAAGGCCTGGGTTg                                                                                                                                 >  1:214565/1‑148 (MQ=255)                         aacGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGc                                                                                                                           >  6:232807/1‑149 (MQ=255)                              cTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATcc                                                                                                                      <  7:36750/149‑1 (MQ=255)                                 ccTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCgt                                                                                                                    <  2:145106/148‑1 (MQ=255)                                 ccTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCgt                                                                                                                    <  1:207530/148‑1 (MQ=255)                                  cTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCgtgt                                                                                                                  >  1:338407/1‑149 (MQ=255)                                         tGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAg                                                                                                           >  3:299986/1‑149 (MQ=255)                                                ccAGGCTCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCcagcag                                                                                                    <  7:292417/149‑1 (MQ=255)                                                ccAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCcagcag                                                                                                    <  8:110114/149‑1 (MQ=255)                                                  aGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCgtgt                                                                                                                  >  1:137452/1‑133 (MQ=255)                                                             ggCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTCGGACCGACCATGCAACAAACGTTTACTACAACTTCCATGTGTTCTAAGGCATGGGTTCCAAGGCCAACCTTTTGTTCGAGCCACCAGAATACGCccatg                                                                                        >  1:258767/1‑149 (MQ=255)                                                                  cGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCccctgccct                                                                                    <  2:258767/147‑1 (MQ=255)                                                                    aCATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGAt                                                                                 >  2:307017/1‑148 (MQ=255)                                                                      aTGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCaa                                                                              <  6:190785/149‑1 (MQ=255)                                                                        ggATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAg                                                                             <  4:299986/148‑1 (MQ=255)                                                                               tgtgTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCaa                                                                              <  6:213359/140‑1 (MQ=255)                                                                               tgtgTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCaa                                                                              >  5:213359/1‑140 (MQ=255)                                                                                 tgtACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACag                                                                       >  2:61433/1‑145 (MQ=255)                                                                                   tACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATg                                                                 >  3:305428/1‑149 (MQ=255)                                                                                   tACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATg                                                                 >  6:163596/1‑149 (MQ=255)                                                                                      aaCTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGcc                                                              >  1:114101/1‑149 (MQ=255)                                                                                       aCTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCccctgcc                                                                                      <  1:141850/124‑1 (MQ=255)                                                                                       aCTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCccctgcc                                                                                      >  2:141850/1‑124 (MQ=255)                                                                                       aCTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCg                                                             <  5:232807/149‑1 (MQ=255)                                                                                           gCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAag                                                          <  1:61433/148‑1 (MQ=255)                                                                                           gCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGCGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGCTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGaa                                                           <  1:90190/147‑1 (MQ=255)                                                                                             gTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGt                                                       >  8:103163/1‑149 (MQ=255)                                                                                                  agaACATCCACTC‑GCCCGACGCTGCCACCTCCGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCcgc                                                  <  4:308300/149‑1 (MQ=255)                                                                                                       aTCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCACTGCCCTGATAAAGGCACAGAAGATGGACGAAGAGTTCCGCGTGcc                                             >  6:330771/1‑149 (MQ=255)                                                                                                            cTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCggtgg                                        >  5:103037/1‑149 (MQ=255)                                                                                                            cTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCggtgg                                        >  4:103262/1‑149 (MQ=255)                                                                                                                gaccgacGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCggtggtg                                      <  2:114101/148‑1 (MQ=255)                                                                                                                   cgacGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGtt                                   >  4:77448/1‑148 (MQ=255)                                                                                                                      cGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAg                               <  6:103037/149‑1 (MQ=255)                                                                                                                      cGCTGACACCTTCGGTGTCGACCACGTCGCTGTGGTAGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCACAAATCGCCCCGGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTGCGCGTGCCGGTGGTGGTTGAAg                               <  2:214565/149‑1 (MQ=255)                                                                                                                        cTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTg                             <  1:307017/149‑1 (MQ=255)                                                                                                                            cACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTTCCGGTGGTGGTTGAAGTGATtc                         >  3:192099/1‑149 (MQ=255)                                                                                                                             aCCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAg                               <  1:317164/142‑1 (MQ=255)                                                                                                                             aCCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAg                               >  2:317164/1‑142 (MQ=255)                                                                                                                             aCCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTGATtct                        >  4:295571/1‑149 (MQ=255)                                                                                                                              ccTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTGATtctc                       <  4:305428/149‑1 (MQ=255)                                                                                                                                tACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCCGTGGTGGTTGAAGTGATTCTCgc                     >  8:114094/1‑148 (MQ=255)                                                                                                                                caggccttCGACCACAGCGACGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTGATTCTCGa                     <  2:232762/142‑1 (MQ=255)                                                                                                                                  cGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTGATTCTCGa                     <  5:163596/147‑1 (MQ=255)                                                                                                                                     cttCGACCACAGCGACGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTGATTCTCGAGCgtg                >  2:132496/3‑149 (MQ=255)                                                                                                                                         gagcACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTGATTCTCGAGCGTGTGa             <  8:107408/145‑1 (MQ=255)                                                                                                                                               gTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTGATTCTCGAGCGTGTGACCAACAt      <  3:77448/149‑1 (MQ=255)                                                                                                                                                   cTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGTTCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTGATTCTCGAGCGTGTGACCAACATTTcc  <  7:317955/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CCGTACGTCCACGTACTGGTGAACAACGCCTACCTGGGCCTGATCCGCCAGGCGCAGCGTGGCTTCGACATGGATTACTGTGTACAACTGGCGTTCGAGAACATCAACTC‑GACCGACGCTGCCACCTACGGTGTCGACCACGTCGCTGTGGTCGAAGGCCTGGGTTGCAAGGCCATCCGTGTGTTCGAGCCAGCAGAAATCGCCCCTGCCCTGATCAAGGCACAGAAGATGGCCGAAGAGTTCCGCGTGCCGGTGGTGGTTGAAGTGATTCTCGAGCGTGTGACCAACATTTCC  >  NC_002947/4888334‑4888627

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