Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,033,040	C→T	5.9%	G226G (GGC→GGT)	PP_3555 →	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,033,040	0	C	T	5.9%	203.1 / 2.3	67	G226G (GGC→GGT)	PP_3555	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base C (36/27);  new base T (2/2);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GTGCGTGCCGGCAGTGGCCAGACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCG  >  NC_002947/4032894‑4033188                                                                                                                                                    |                                                                                                                                                     gtgcgtgcCGGCAGTGGCCAGACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGACGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTGTAGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGc                                                                                                                                                      <  1:247414/148‑1 (MQ=255)     gtgcCGGCAGTGGCCAGACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAg                                                                                                                                                >  6:114154/1‑149 (MQ=255)        ccGGCAGTGGCCAGACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGcc                                                                                                                                             <  1:18966/149‑1 (MQ=255)        ccGGCAGTGGCCAGACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGc                                                                                                                                              >  4:69977/1‑148 (MQ=255)                 gCCAGACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACcac                                                                                                                                    >  6:216377/1‑149 (MQ=255)                   cAGACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCACGACCAGGTGAAGGTCACCCTGGTGGAAGGTTAAGTCCTGGTGTCCAGTG‑ACGGCAGTGACGGTGGCTACCGCCTTGGCACTGGCAAGCAGGCAAGCTACcac                                                                                                                                    >  8:217152/1‑147 (MQ=255)                    aGACCCGGGTCACCGGCACCCAGTTCAACGTTTGGAAGTACCAGGAACAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGTATGCAGTCCAGCTACCACAAg                                                                                                                                 >  5:183186/1‑149 (MQ=255)                    aGACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGTATGCAGGCCAGCTACCACAAg                                                                                                                                 >  5:188581/1‑149 (MQ=255)                     gACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAgg                                                                                                                                <  4:348769/149‑1 (MQ=255)                     gACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAgg                                                                                                                                <  2:139901/149‑1 (MQ=255)                      aCCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAggg                                                                                                                               >  3:319301/1‑149 (MQ=255)                         cGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGACTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGa                                                                                                                            >  5:317846/1‑149 (MQ=255)                               aCCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTc                                                                                                                        >  2:136493/1‑147 (MQ=255)                                 cGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACtt                                                                                                                         <  8:121913/144‑1 (MQ=255)                                 cGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACtt                                                                                                                         >  7:121913/1‑144 (MQ=255)                                 cGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGc                                                                                                                    >  8:283053/1‑149 (MQ=255)                                   gCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCg                                                                                                                  <  5:114154/149‑1 (MQ=255)                                      cccAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGATCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAg                                                                                                               >  5:45765/1‑149 (MQ=255)                                       ccAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGc                                                                                                              >  5:63495/1‑149 (MQ=255)                                              aaCGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGc                                                                                                        >  5:345353/1‑148 (MQ=255)                                                  tgtgGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGc                                                                                                   <  6:177339/149‑1 (MQ=255)                                                    tgGAAGTACCAGGACCAGGTGAAGGTCACCCTTGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCAGGAGCAAAGcga                                                                                                 >  5:159870/1‑149 (MQ=255)                                                         gTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAACGCgacgac                                                                                             >  2:348303/1‑148 (MQ=255)                                                            ccaggaccaggTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATg                                                                                         <  4:293230/149‑1 (MQ=255)                                                            ccaggaccaggTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATg                                                                                         >  3:109856/1‑149 (MQ=255)                                                             caggaccaggTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGCCTATgg                                                                                        >  5:299236/1‑149 (MQ=255)                                                             caggaccaggTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATgg                                                                                        <  5:216377/149‑1 (MQ=255)                                                              aggaccaggTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGcgacg                                                                                               >  5:86011/1‑141 (MQ=255)                                                              aggaccaggTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGcgacg                                                                                               <  6:86011/141‑1 (MQ=255)                                                              aggaccaggTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATgg                                                                                        >  5:34177/1‑148 (MQ=255)                                                              aggaccaggTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGACTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGc                                                                                       <  6:317846/149‑1 (MQ=255)                                                                gaccaggTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCaa                                                                                     <  6:159870/149‑1 (MQ=255)                                                                     ggTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGcc                                                                                <  8:101766/149‑1 (MQ=255)                                                                     ggTGAAGGTCACCCTGGTAGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCaaa                                                                                                     <  6:27935/128‑1 (MQ=255)                                                                     ggTGAAGGTCACCCTGGTAGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCaaa                                                                                                     >  5:27935/1‑128 (MQ=255)                                                                           ggTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCAt                                                                          >  2:121562/1‑149 (MQ=255)                                                                               aCCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGgcg                                                                      >  3:19228/1‑149 (MQ=255)                                                                                 ccTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGgcg                                                                      >  6:74821/1‑147 (MQ=255)                                                                                  cTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAg                                                                   <  8:309852/149‑1 (MQ=255)                                                                                        gAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGc                                                               >  4:197824/1‑147 (MQ=255)                                                                                        gAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGc                                                               >  1:243443/1‑147 (MQ=255)                                                                                            gTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTg                                                         <  6:45765/149‑1 (MQ=255)                                                                                             ttCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTgg                                                        >  6:48678/1‑149 (MQ=255)                                                                                             ttCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTgg                                                        <  4:319301/149‑1 (MQ=255)                                                                                                aGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTg                                                      <  6:88953/148‑1 (MQ=255)                                                                                                  tGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTc                                                   <  1:136493/149‑1 (MQ=255)                                                                                                   gCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTc                                                   >  1:270181/1‑148 (MQ=255)                                                                                                       gtgtCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCCACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGAc                                                <  1:121562/147‑1 (MQ=255)                                                                                                        tgtCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAAc                                             >  6:100953/1‑149 (MQ=255)                                                                                                            cAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGc                                                           >  5:89026/1‑131 (MQ=255)                                                                                                            cAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGc                                                           <  6:89026/131‑1 (MQ=255)                                                                                                                  aGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTg                                    <  6:345353/149‑1 (MQ=255)                                                                                                                  aGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTg                                    <  7:283053/149‑1 (MQ=255)                                                                                                                    ggCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTgg                                   >  4:327332/1‑148 (MQ=255)                                                                                                                            cGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGgcgc                          >  2:176491/1‑149 (MQ=255)                                                                                                                              gTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTg                        >  2:220324/1‑149 (MQ=255)                                                                                                                                    aCCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGCTGCTCGACAACCTGCGCCTGGAACAGGCGCTGCCGGt                   >  7:245157/1‑148 (MQ=255)                                                                                                                                     ccGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGt                   <  3:69977/147‑1 (MQ=255)                                                                                                                                     ccGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTg                  <  4:109856/148‑1 (MQ=255)                                                                                                                                     ccGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTg                  <  3:197824/148‑1 (MQ=255)                                                                                                                                      cGCCTTGGGCCTGGTATGCAGGCCAGCTACCACAAGGGCGAGTTCGAGGCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGAt                <  6:188581/149‑1 (MQ=255)                                                                                                                                      cGCCTTGGCCCTGGTATGCAGGCCAGCTACCACAAGGGCGACTTCTAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGTCATGTCGCAGCGGCAATCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTTCCGTTGAt                <  6:183186/149‑1 (MQ=255)                                                                                                                                      cGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGAt                <  1:176491/149‑1 (MQ=255)                                                                                                                                       gCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAg                                                                   <  8:304233/97‑1 (MQ=255)                                                                                                                                              cccTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGc         >  7:286014/1‑148 (MQ=255)                                                                                                                                                 tGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTAcc     >  5:341161/1‑149 (MQ=255)                                                                                                                                                 tGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCACTAc      >  8:323145/1‑148 (MQ=255)                                                                                                                                                   gCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTc   >  2:322441/1‑149 (MQ=255)                                                                                                                                                    cATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCg  <  5:100953/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTGCGTGCCGGCAGTGGCCAGACCCGGGTCACCGGCACCCAGTTCAACGTGTGGAAGTACCAGGACCAGGTGAAGGTCACCCTGGTGGAAGGTTCAGTGCTGGTGTCCAGTG‑AGGGCAGTGCCGGTGGCTACCGCCTTGGCCCTGGCATGCAGGCCAGCTACCACAAGGGCGACTTCGAGCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCG  >  NC_002947/4032894‑4033188

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: ‑