Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,647,085	A→T	5.9%	L91Q (CTG→CAG)	gcd ←	quinoprotein glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,647,085	0	A	T	5.9%	207.3 / 2.7	68	L91Q (CTG→CAG)	gcd	quinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base T (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

TGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAG  >  NC_002947/1646939‑1647233                                                                                                                                                    |                                                                                                                                                          tggcCAGCGCGCAAGCGCCCGAAGCGACGACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGt                                                                                                                                                         >  1:206800/1‑149 (MQ=255) tggcCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGt                                                                                                                                                         >  1:201144/1‑149 (MQ=255)       gcgcgcAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGccac                                                                                                                                               >  8:138719/1‑149 (MQ=255)           gcAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGAGGGAGCACCAGAAAGCCAGCCGTGGCACCAGG‑‑‑‑TGGCAACA‑Gt                                                                                                                                          >  1:307830/1‑149 (MQ=255)            cAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCAACTGG‑‑‑‑TGccacggctg                                                                                                                                          <  4:256614/149‑6 (MQ=255)                           aCCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCGCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGc                                                                                                                                    >  1:266820/1‑140 (MQ=255)                           aCCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGc                                                                                                                                    <  2:266820/139‑1 (MQ=255)                           aCCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGaa                                                                                                                          <  7:138719/149‑1 (MQ=255)                                tGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCg                                                                                                                            >  6:357752/1‑142 (MQ=255)                                tGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCg                                                                                                                            <  5:357752/142‑1 (MQ=255)                                tGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAggg                                                                                                                     >  8:177828/1‑149 (MQ=255)                                  ccACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAggg                                                                                                                     <  8:330620/147‑1 (MQ=255)                                  ccACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAggg                                                                                                                     <  1:141224/147‑1 (MQ=255)                                               tgcagtgTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGc                                                                                                                                                  <  6:121270/106‑1 (MQ=255)                                               tgcagtgTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGc                                                                                                                                                  >  5:121270/1‑106 (MQ=255)                                                      ttGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCcagcag                                                                                               <  5:44103/149‑1 (MQ=255)                                                            ttGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAg                                                                                         >  7:136498/1‑149 (MQ=255)                                                             tGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAgg                                                                                        >  5:79199/1‑149 (MQ=255)                                                                 ggctggGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCa                                                                                     >  1:94585/1‑148 (MQ=255)                                                                      ggC‑CGATCAGCGGGCGACGGGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAgg                                                                               <  1:354840/149‑1 (MQ=255)                                                                          cGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGccc                                                                            <  8:332210/149‑1 (MQ=255)                                                                           gATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTAAGGGGGCACCATTTACAGGCCAAATTAGAAAAGGGCCCCACCGGTGGTGCCCACCAGCCCCAGGCCaaaaag                                                                                >  1:114933/1‑144 (MQ=255)                                                                               aGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAggg                                                                        <  6:79199/148‑1 (MQ=255)                                                                                  gggCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGAt                                                                    >  6:280386/1‑149 (MQ=255)                                                                                   ggCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATc                                                                   >  3:105243/1‑149 (MQ=255)                                                                                   ggCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATc                                                                   >  3:301555/1‑149 (MQ=255)                                                                                            cccACGGCAGCAGCAACACCACGCCGATGGCGAACCATATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTCCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                          >  6:84987/1‑149 (MQ=255)                                                                                            cccACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                          <  2:3762/149‑1 (MQ=255)                                                                                            cccACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                          <  2:94585/149‑1 (MQ=255)                                                                                             ccACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTggg                                                         <  5:280386/149‑1 (MQ=255)                                                                                              cACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACTCGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGc                                                        <  2:307830/149‑1 (MQ=255)                                                                                              cACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCTGGGCAGCAcc                                                                                           <  2:347305/114‑1 (MQ=255)                                                                                              cACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCTGGGCAGCAcc                                                                                           >  1:347305/1‑114 (MQ=255)                                                                                                  gcagcagcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagc                                                    >  3:15055/1‑149 (MQ=255)                                                                                                     gcagcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagc                                                 <  7:177828/149‑1 (MQ=255)                                                                                                       agcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagcag                                               <  8:136498/149‑1 (MQ=255)                                                                                                          aaCACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGAc                                            >  7:29395/1‑149 (MQ=255)                                                                                                          aaCACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGAc                                            >  4:78235/1‑149 (MQ=255)                                                                                                            caccacGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGc                                          >  5:225582/1‑149 (MQ=255)                                                                                                               cacGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGCTCTGGCGTTGGGCCAGCAGCAGGACGccc                                        >  4:262371/1‑148 (MQ=255)                                                                                                                acGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCTTAACGCCCg                                       >  1:128202/1‑148 (MQ=255)                                                                                                                    cGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAgg                                  >  5:266801/1‑149 (MQ=255)                                                                                                                      aTGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAgg                                  >  4:1850/1‑147 (MQ=255)                                                                                                                      aTGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAgg                                  <  3:1850/147‑1 (MQ=255)                                                                                                                           gAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGc                           >  7:278182/1‑149 (MQ=255)                                                                                                                           gAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGc                           <  4:105243/149‑1 (MQ=255)                                                                                                                              ccAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCaa                        >  6:121227/1‑149 (MQ=255)                                                                                                                               cAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCaa                        >  1:234399/1‑148 (MQ=255)                                                                                                                                 gATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagc                                                 >  7:102041/1‑121 (MQ=255)                                                                                                                                 gATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagc                                                 <  8:102041/121‑1 (MQ=255)                                                                                                                                     gccagccGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGc                 <  2:234399/149‑1 (MQ=255)                                                                                                                                     gccacccGTTGCACCAGA‑‑‑‑TGACACCA‑GCCCACGCCTACTTCGAACCGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCTATCTGGCGTTGGTCCAGCAGCAGGACGCCCGAAAGCGCAAAGCCAATACCGGc                 <  8:282580/149‑1 (MQ=255)                                                                                                                                      ccagccGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCGAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCg                >  7:229939/1‑149 (MQ=255)                                                                                                                                      ccagccGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCGGGACGCCCGAAAGGGCAAACCCAATACCGGc                 >  4:225296/1‑148 (MQ=255)                                                                                                                                        agccGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCCGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGAt              <  2:119699/149‑1 (MQ=255)                                                                                                                                          tggTGGCAACTGG‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagca                                                   >  4:107269/3‑110 (MQ=255)                                                                                                                                          tggTGGCAACTGG‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagca                                                   <  3:107269/108‑1 (MQ=255)                                                                                                                                          tggTGGCAACTGG‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGAt              >  3:167128/3‑147 (MQ=255)                                                                                                                                            gTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                          >  2:41663/1‑101 (MQ=255)                                                                                                                                            gTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                          <  1:41663/101‑1 (MQ=255)                                                                                                                                            gTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAgg          <  5:84987/149‑1 (MQ=255)                                                                                                                                              ggCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgta        <  2:206800/149‑1 (MQ=255)                                                                                                                                               gCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgtag       >  6:28514/1‑149 (MQ=255)                                                                                                                                               gCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgta        <  5:121227/148‑1 (MQ=255)                                                                                                                                               gCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgta        <  2:201144/148‑1 (MQ=255)                                                                                                                                                   cAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACa   >  4:335177/1‑149 (MQ=255)                                                                                                                                                   cAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACa   <  8:29395/149‑1 (MQ=255)                                                                                                                                                   cAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACa   <  5:28514/149‑1 (MQ=255)                                                                                                                                                    aGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAg  <  3:335177/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                          TGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGC‑CGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCA‑GTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAG  >  NC_002947/1646939‑1647233

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