Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,886,602	C→T	8.3%	D423D (GAC→GAT)	rmuC →	DNA recombination limiting inversions protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,886,602	0	C	T	8.3%	174.3 / 4.2	60	D423D (GAC→GAT)	rmuC	DNA recombination limiting inversions protein
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base T (3/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

TGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACCT  >  NC_002947/1886454‑1886748                                                                                                                                                      |                                                                                                                                                   tGATTGTCAGCCCTACCACGTTGCTGGCCACCTTTCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGAc                                                                                                                                                    >  6:148878/1‑149 (MQ=255)       tCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACt                                                                                                                                               <  6:229355/148‑1 (MQ=255)       tCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTg                                                                                                                                              >  2:60666/1‑149 (MQ=255)         aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTggg                                                                                                                                            <  8:224507/149‑1 (MQ=255)         aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTggg                                                                                                                                            >  6:193022/1‑149 (MQ=255)         aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTgg                                                                                                                                             <  5:174845/148‑1 (MQ=255)          gCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGc                                                                                                                                           >  4:760/1‑149 (MQ=255)            ccTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAg                                                                                                                                         >  2:72975/1‑149 (MQ=255)                  aCGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCtt                                                                                                                                   <  4:83722/149‑1 (MQ=255)                    gTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGc                                                                                                                                     <  4:150679/145‑1 (MQ=255)                    gTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGc                                                                                                                                     >  3:150679/1‑145 (MQ=255)                    gTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTgc                                                                                                                                 <  2:280323/149‑1 (MQ=255)                         tGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGcagcag                                                                                                                            >  8:324014/1‑149 (MQ=255)                         tGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGcagcag                                                                                                                            >  2:278960/1‑149 (MQ=255)                            ccACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTg                                                                                                                         <  5:173798/149‑1 (MQ=255)                            ccACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTg                                                                                                                         <  1:2244/149‑1 (MQ=255)                                 ttGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCtt                                                                                                                                   >  1:309209/1‑134 (MQ=255)                                 ttGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCtt                                                                                                                                   <  2:309209/134‑1 (MQ=255)                                      ggTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCtt                                                                                                               <  4:101380/149‑1 (MQ=255)                                      ggTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCt                                                                                                                >  4:340260/1‑148 (MQ=255)                                        ttATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTAt                                                                                                             <  7:324014/149‑1 (MQ=255)                                              aCAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCa                                                                                                       <  5:148878/149‑1 (MQ=255)                                                   cTGTGGAAGCAGGAGCGTCAAGGGCATAACGCCCGCGAAATTGCCGCTCGGGCCGGCTGGCTCTACGACAATTTCGTGCTGTTCATCCAGGACCTGGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAgcgc                                                                                                   <  1:278960/149‑1 (MQ=255)                                                                    tCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAg                                                                                 <  6:209813/149‑1 (MQ=255)                                                                      aaGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAgg                                                                                <  1:72975/148‑1 (MQ=255)                                                                          gCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTg                                                                           <  7:220705/149‑1 (MQ=255)                                                                             gAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGc                                                                         >  5:330912/1‑148 (MQ=255)                                                                                    cgcgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCTGCa                                                                                                       <  4:65448/111‑1 (MQ=255)                                                                                    cgcgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCa                                                                                                       >  3:65448/1‑111 (MQ=255)                                                                                     gcgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                                 >  5:330295/1‑148 (MQ=255)                                                                                      cgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTa                                                               <  6:348016/149‑1 (MQ=255)                                                                                      cgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTa                                                               <  4:170946/149‑1 (MQ=255)                                                                                          aTTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTAGTGCTGTTCATCCAGGACCT‑GGACGGACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCg                                                           >  7:311930/1‑149 (MQ=255)                                                                                              ccGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGc                                                       <  5:65226/149‑1 (MQ=255)                                                                                              ccGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGc                                                       >  7:148864/1‑149 (MQ=255)                                                                                                gAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGa                                                     <  3:760/149‑1 (MQ=255)                                                                                                       ccGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCt                                               >  6:193132/1‑148 (MQ=255)                                                                                                           ctggctCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGagcag                                                 >  6:299989/1‑142 (MQ=255)                                                                                                           ctggctCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGagcag                                                 <  5:299989/142‑1 (MQ=255)                                                                                                                tctACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGctgctg                                     >  3:133753/1‑149 (MQ=255)                                                                                                                     gACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTgg                                    <  4:63583/145‑1 (MQ=255)                                                                                                                     gACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgcgc                                <  5:193022/149‑1 (MQ=255)                                                                                                                     gACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgcgc                                >  3:99292/1‑149 (MQ=255)                                                                                                                     gACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTCGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTgg                                    >  3:63583/1‑145 (MQ=255)                                                                                                                           ttCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc                          <  5:193132/149‑1 (MQ=255)                                                                                                                           ttCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc                          <  2:151198/149‑1 (MQ=255)                                                                                                                           ttCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGAAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc                          >  3:282610/1‑149 (MQ=255)                                                                                                                           ttCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGCCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc                          >  1:151198/1‑149 (MQ=255)                                                                                                                                gCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc                      >  1:139831/1‑148 (MQ=255)                                                                                                                                ccagTTCGTCCAGGTCCT‑GGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                                   >  2:59330/4‑103 (MQ=39)                                                                                                                                ccagTTCGTCCAGGTCCT‑GGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                                   >  2:53867/4‑103 (MQ=39)                                                                                                                                ccagTTCGTCCAGGTCCT‑GGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                                   <  1:59330/100‑1 (MQ=39)                                                                                                                                ccagTTCGTCCAGGTCCT‑GGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                                   <  1:53867/100‑1 (MQ=39)                                                                                                                                  tGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc                      >  1:8378/1‑146 (MQ=255)                                                                                                                                   gTTCGTCCAGGTCCT‑GGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCaa                    >  5:122059/1‑147 (MQ=255)                                                                                                                                   gTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGGAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCaa                    >  7:44228/1‑147 (MQ=255)                                                                                                                                      cATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGc                <  6:330295/148‑1 (MQ=255)                                                                                                                                       aTCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc                      >  4:225613/1‑141 (MQ=255)                                                                                                                                       aTCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc                      <  3:225613/141‑1 (MQ=255)                                                                                                                                       aTCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCtt              <  4:99292/149‑1 (MQ=255)                                                                                                                                                    gACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACct  >  7:193961/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCT‑GGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACCT  >  NC_002947/1886454‑1886748

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