Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	45,965	T→C	6.6%	Q29R (CAA→CGA)	PP_0042 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	45,965	0	T	C	6.6%	172.8 / 2.3	60	Q29R (CAA→CGA)	PP_0042	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/25);  new base C (2/2);  total (33/27)
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.79e-01

CTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAGCAATGCATAGGCATCACCGG  >  NC_002947/45820‑46111                                                                                                                                                  |                                                                                                                                                   cTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGaa                                                                                                                                                 <  5:271866/149‑1 (MQ=255)              gggTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCtctc                                                                                                                                    <  8:137206/149‑1 (MQ=255)                      cGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGaaa                                                                                                                             <  6:242275/148‑1 (MQ=255)                        aCTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATCGAAGCTGCTTGCTCTCGGAGAAAgtg                                                                                                                          >  2:314542/1‑149 (MQ=255)                         cTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAgtgt                                                                                                                         <  5:225354/149‑1 (MQ=255)                            aCGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGaaa                                                                                                                             >  8:286236/1‑142 (MQ=255)                            aCGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGaaa                                                                                                                             <  7:286236/142‑1 (MQ=255)                                        cGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGctc                                                                                                                                      >  4:258484/1‑121 (MQ=255)                                        cGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGctc                                                                                                                                      <  3:258484/121‑1 (MQ=255)                                          cAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCAcc                                                                                                        >  5:148848/1‑149 (MQ=255)                                               tCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGaa                                                                                                   <  1:359436/149‑1 (MQ=255)                                                 tATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGaaa                                                                                                  >  3:12338/1‑148 (MQ=255)                                                   tCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGaaaa                                                                                                 <  6:45020/147‑1 (MQ=255)                                                   tCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAaac                                                                                                >  7:108117/1‑148 (MQ=255)                                                        aTGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAg                                                                                          >  2:53154/1‑149 (MQ=255)                                                             atatTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTgcag                                                                                     <  1:358422/149‑1 (MQ=255)                                                                tttGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCgaga                                                                                                    <  1:82998/131‑1 (MQ=255)                                                                tttGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCgaga                                                                                                    >  2:82998/1‑131 (MQ=255)                                                                tttGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGcagcag                                                                                  >  8:278305/1‑149 (MQ=255)                                                                tttGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGcagcag                                                                                  >  3:191585/1‑149 (MQ=255)                                                                     aaTCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATCGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGAt                                                                              <  1:314542/148‑1 (MQ=255)                                                                      aTCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATaa                                                                            >  6:266631/1‑149 (MQ=255)                                                                      aTCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTACTTGCTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATaa                                                                            >  6:261711/1‑149 (MQ=255)                                                                          tGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAg                                                                        >  8:3616/1‑149 (MQ=255)                                                                                tctGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGtt                                                                                                                        <  5:113748/95‑1 (MQ=255)                                                                                tctGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGtt                                                                                                                        >  6:113748/1‑95 (MQ=255)                                                                                  tGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCa                                                                >  7:45709/1‑149 (MQ=255)                                                                                  tGGCATTGGCTTGCCATTTATCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCttt                                                                    <  3:53014/145‑1 (MQ=255)                                                                                    gCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGATTTTTCac                                                               >  7:237550/1‑147 (MQ=255)                                                                                      aTTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATaaa                                                            >  8:110854/1‑149 (MQ=255)                                                                                      aTTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATaa                                                             >  1:47963/1‑148 (MQ=255)                                                                                          gCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAAcc                                                         <  8:108117/148‑1 (MQ=255)                                                                                                 aTTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGaaaa                                                                                                 <  7:272495/101‑1 (MQ=255)                                                                                                 aTTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGaaaa                                                                                                 >  8:272495/1‑101 (MQ=255)                                                                                                 aTTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAAt                                                 <  4:191585/149‑1 (MQ=255)                                                                                                       ctcCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCaa                                            <  4:12338/148‑1 (MQ=255)                                                                                                       ctcCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCaa                                            <  7:278305/148‑1 (MQ=255)                                                                                                          cgcgTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCaaa                                           <  6:148848/146‑1 (MQ=255)                                                                                                           gcgTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACt                                        >  2:169598/1‑148 (MQ=255)                                                                                                             gTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTtct                                     >  6:70715/1‑149 (MQ=255)                                                                                                                cTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATCGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCgg                                  >  5:343674/1‑149 (MQ=255)                                                                                                                 ttGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCggg                                 >  3:106138/1‑149 (MQ=255)                                                                                                                      cACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAAt                                                 >  7:272071/1‑128 (MQ=255)                                                                                                                      cACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAAt                                                 <  8:272071/128‑1 (MQ=255)                                                                                                                      cACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTaa                            >  2:287608/1‑149 (MQ=255)                                                                                                                      cACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTaa                            >  6:161798/1‑149 (MQ=255)                                                                                                                           cGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCt                         >  4:183756/1‑147 (MQ=255)                                                                                                                            ggTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAg                      <  7:3616/149‑1 (MQ=255)                                                                                                                                  gtgtTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCctc                                                      >  3:113498/1‑111 (MQ=255)                                                                                                                                  gtgtTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCctc                                                      <  4:113498/111‑1 (MQ=255)                                                                                                                                  gtgtTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAGCaa                   <  5:161798/146‑1 (MQ=255)                                                                                                                                  gtgtTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAGCAATg                 >  5:134316/1‑148 (MQ=255)                                                                                                                                  gtgtTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAGCAATg                 >  6:176277/1‑148 (MQ=255)                                                                                                                                  gtgtTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAGCAATGc                >  6:181470/1‑149 (MQ=255)                                                                                                                                     ttCCTTGGTGTATCGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAGCAATGCATa             <  6:343674/149‑1 (MQ=255)                                                                                                                                       ccTTGGTGTATTGAAGCTGCTTTCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAGCAATGCATAg            <  5:70715/148‑1 (MQ=255)                                                                                                                                         ttGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTa                       >  3:172498/1‑135 (MQ=255)                                                                                                                                         ttGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTa                       <  4:172498/135‑1 (MQ=255)                                                                                                                                             tgtATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAAGTTCTCGGGCGTAACCTTAGCAATGCATAGGCATCAc     >  5:344904/1‑149 (MQ=255)                                                                                                                                               tATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAGCAATGCATAGGCATCAcc    <  3:183756/148‑1 (MQ=255)                                                                                                                                                 ttGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAGCAATGCATAGGCATCACCgg  >  3:181031/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CTGGCAGGTACGGGGGTCCTGCGACTGACGTACGAATTTCGCAACGTCTATCTTCATGCCATATTTGTAATCCTGAATCTCTGGCATTGGCTTGCCATTTTTCTCCGCGTACTTGGCCACTGCGGTCTGGTGTTCCTTGGTGTATTGAAGCTGCTTGCTCTCGGAGAAAGTGTTGGCCTGCGAAGCCACCGAGAAAACAACAAGTGCAGCAGCGATAATTAGCTTTTTCATAAAACCCTCTCAATTTCAAAACTTCTCGGGCGTAACCTTAGCAATGCATAGGCATCACCGG  >  NC_002947/45820‑46111

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