Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,693,803	G→A	12.4%	A32T (GCC→ACC)	PP_2362 →	usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,693,803	0	G	A	12.4%	147.2 / 7.2	59	A32T (GCC→ACC)	PP_2362	usher protein
Reads supporting (aligned to +/- strand):  ref base G (28/22);  new base A (2/5);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

TGGTCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGCGGGCCGG  >  NC_002947/2693658‑2693951                                                                                                                                                  |                                                                                                                                                     tGGTCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCa                                                                                                                                                   >  6:108160/1‑149 (MQ=255)   gTCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATg                                                                                                                                                 <  5:220196/149‑1 (MQ=255)    tCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCGGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCCGGGGCTGACCATgg                                                                                                                                                >  8:88139/1‑149 (MQ=255)    tCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATg                                                                                                                                                 >  5:44908/1‑148 (MQ=255)    tCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATg                                                                                                                                                 <  6:44908/148‑1 (MQ=255)         acTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGc                                                                                                                                           >  6:271334/1‑149 (MQ=255)           tCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGccc                                                                                                                                         >  7:266547/1‑149 (MQ=255)             ggcggACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCcg                                                                                                                                       >  2:234146/1‑149 (MQ=255)             ggcggACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGACCATGGCCAGcccca                                                                                                                                       <  5:259622/149‑2 (MQ=255)             ggcggACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGACGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCcg                                                                                                                                       <  1:234146/149‑1 (MQ=255)              gcgtACGCCATGCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCATGGGCCAGGGAGGACCCGGCAATGGTTGTAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGACCATGGCCAGccccag                                                                                                                                      <  6:104022/149‑3 (MQ=255)              gcggACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGACCATGGCCAGccccag                                                                                                                                      <  7:168717/149‑3 (MQ=255)               cggACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCCGGGGCTGACCATGGCCAGccccagg                                                                                                                                     <  7:88139/149‑4 (MQ=255)                ggACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCcgcac                                                                                                                                    <  7:280909/149‑3 (MQ=255)                ggACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCgcgc                                                                                                                                    <  7:160275/149‑1 (MQ=255)                  aCGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTg                                                                                                                                  <  5:271334/149‑1 (MQ=255)                      cATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGACCATGGCCAGccc                                                                                                                                         >  8:10360/1‑138 (MQ=255)                      cATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGACCATGGCCAGccc                                                                                                                                         <  7:10360/138‑1 (MQ=255)                             cacggCCAATGATCCAACGGCCTAATGACGTGCAGGGGCCCGGGCGGCCCCGGCAATGGTGGGAAACCGTGGGTGAAATGGCTCCCGGCGACGACGTGCCGCGGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCcga                                                                                                                       <  8:153438/146‑1 (MQ=255)                              aaGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCcgcg                                                                                                                                     <  3:199377/134‑1 (MQ=255)                              aaGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCcgcg                                                                                                                                     >  4:199377/1‑134 (MQ=255)                              aaGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTg                                                                                                                            <  4:305718/143‑1 (MQ=255)                              aaGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTg                                                                                                                            >  3:305718/1‑143 (MQ=255)                              aaGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCcgac                                                                                                                      <  2:133597/149‑1 (MQ=255)                                   cAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACtt                                                                                                                 >  3:61926/1‑149 (MQ=255)                                   cAATGAGCCAACGGACGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGAGAAATGGCTCCCGGCGACGACGTGCCGCTCGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGCTACTGGCCGACGACtt                                                                                                                 >  3:274292/1‑149 (MQ=255)                                         gCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTgg                                                                                                                                 <  1:285350/127‑1 (MQ=255)                                         gCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTgg                                                                                                                                 >  2:285350/1‑127 (MQ=255)                                          ccAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCa                                                                                                          >  1:241439/1‑149 (MQ=255)                                           cAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCAc                                                                                                         >  7:134314/1‑149 (MQ=255)                                            aaCGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCAcc                                                                                                        >  2:314934/1‑149 (MQ=255)                                                        aaGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGgaga                                                                                             <  4:61926/148‑1 (MQ=255)                                                         aGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCGGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGAt                                                                                            <  3:274074/148‑1 (MQ=255)                                                                  gCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCACACGGAGATGTCCGCCATc                                                                                  >  6:246086/1‑149 (MQ=255)                                                                     aGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGcc                                                                               >  7:303091/1‑149 (MQ=255)                                                                       ggaggaCCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCGTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGa                                                                             <  1:201633/149‑1 (MQ=255)                                                                                ggCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTgg                                                                                                                                 >  7:170600/1‑88 (MQ=255)                                                                                ggCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCt                                                                    >  1:9828/1‑149 (MQ=255)                                                                                    ttGGTTGTAAACGTTGTTTTAATTGGGTCCGGGCGAGGAGGTGCCGGTGGTGCCTGGGGCTGCCATTGGGCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTAc                                                                <  1:314934/148‑1 (MQ=255)                                                                                      ggTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCt                                                              >  2:229651/1‑149 (MQ=255)                                                                                             aaCCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCtg                                                       >  6:250408/1‑149 (MQ=255)                                                                                              aCCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCtgg                                                      <  8:134314/149‑1 (MQ=255)                                                                                              aCCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCtgg                                                      <  2:9828/149‑1 (MQ=255)                                                                                                      tgAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAAc                                              <  5:250408/149‑1 (MQ=255)                                                                                                          aggccTCCCGCCGCCGCCGTGCGGCGGGTGCCTGGGGCGGCCCCTGGGCCCCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGa                                          <  8:266547/145‑1 (MQ=255)                                                                                                            ggCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATg                                        >  2:266592/1‑149 (MQ=255)                                                                                                             gCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGc                                       <  5:108160/149‑1 (MQ=255)                                                                                                               tCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCtggtg                                                    <  3:20067/134‑1 (MQ=255)                                                                                                               tCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCtggtg                                                    >  4:20067/1‑134 (MQ=255)                                                                                                               tCCCGGCGACGACCTCCACCTGGTCCCACGGGCTGGCCATCCTCAGCCCCGCGCTCATACTGGCCGACGACTTGCCCCCACCGCCCACCCAGATGTCCGCCATCGCCGACGCCACCCTATACCTCGACCTGGTGGTGAACCAGATGcc                                      >  6:327521/1‑148 (MQ=255)                                                                                                                  cGGCGACGCCGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGCGGCCACCCTTTACCTCGACCTGGTGGTGCACCAGATGCCCAgg                                  >  8:113731/1‑149 (MQ=255)                                                                                                                      gacgacGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGtt                               >  1:215230/1‑148 (MQ=255)                                                                                                                          acgTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCt                          <  1:229651/149‑1 (MQ=255)                                                                                                                               ccGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATAc                     <  8:303091/149‑1 (MQ=255)                                                                                                                                   tGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGt                 >  2:20712/1‑149 (MQ=255)                                                                                                                                        ccTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTG‑CCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTgcag             >  8:278711/1‑147 (MQ=255)                                                                                                                                         cTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGa                                          >  7:64898/1‑118 (MQ=255)                                                                                                                                         cTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGa                                          <  8:64898/118‑1 (MQ=255)                                                                                                                                         cTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTgcag             >  3:2199/1‑147 (MQ=255)                                                                                                                                                 ggCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAggg                                 <  6:195016/119‑1 (MQ=255)                                                                                                                                                 ggCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAggg                                 >  5:195016/1‑119 (MQ=255)                                                                                                                                                  gCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGcgggccgg  >  1:13117/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TGGTCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGCGGGCCGG  >  NC_002947/2693658‑2693951

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