Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,366,486	A→T	6.8%	R83R (CGT→CGA)	PP_3843 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,366,486	0	A	T	6.8%	170.0 / 3.1	59	R83R (CGT→CGA)	PP_3843	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/29);  new base T (2/2);  total (28/31)
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.59e-01

ATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCT  >  NC_002947/4366341‑4366634                                                                                                                                                   |                                                                                                                                                     aTAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGt                                                                                                                                                   <  3:219106/149‑1 (MQ=255)  tAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTc                                                                                                                                                  <  1:311523/149‑1 (MQ=255)           gtgtGGCGACGCAGTCTTTTTACCCACTGCTTCTTCTCTGATTACTCTAAGTCGCCGCCCGACCACCGGCGCTGCGCGCCTCCGCTATTTCTAGGCAGCCAGAAACCCTGCCGCAGGCGCACCTAG‑AGCGGCGACACGTCGGGCCCgg                                                                                                                                          <  7:222429/148‑1 (MQ=255)              tgGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCg                                                                                                                                           >  8:294495/1‑144 (MQ=255)              tgGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCg                                                                                                                                           <  7:294495/144‑1 (MQ=255)                gCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCCCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCg                                                                                                                                                 <  5:53864/136‑1 (MQ=255)                gCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCCCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCg                                                                                                                                                 >  6:53864/1‑136 (MQ=255)                    cGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGTCCAACGcccc                                                                                                                                <  2:242784/149‑1 (MQ=255)                         gCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAAc                                                                                                                                     <  6:179665/139‑1 (MQ=255)                         gCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAAc                                                                                                                                     >  5:179665/1‑139 (MQ=255)                         gCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGc                                                                                                                            >  4:299037/1‑148 (MQ=255)                               cACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGcc                                                                                                                     >  5:300506/1‑149 (MQ=255)                                aCGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGccc                                                                                                                    >  7:331239/1‑149 (MQ=255)                                  gCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTg                                                                                                                  <  6:91432/149‑1 (MQ=255)                                  gCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTg                                                                                                                  <  2:99782/149‑1 (MQ=255)                                  gCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTg                                                                                                                  >  5:91432/1‑149 (MQ=255)                                     cTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCGAGAAAGCCTGCCGCAGGCGCACCTTC‑AGTGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTc                                                                                                               <  3:145691/149‑1 (MQ=255)                                        cGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATc                                                                                                            <  5:280380/149‑1 (MQ=255)                                              cGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCg                                                                                                      >  2:269925/1‑149 (MQ=255)                                              cGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCg                                                                                                      >  6:216927/1‑149 (MQ=255)                                                 tgatgaCAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCTGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGc                                                                                                     >  3:5624/1‑147 (MQ=255)                                                   atgaCAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAgg                                                                                                 <  3:69848/149‑1 (MQ=255)                                                        aGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTTTAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGc                                                                                                                                         <  2:312012/104‑1 (MQ=255)                                                        aGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGc                                                                                                                                         >  1:312012/1‑104 (MQ=255)                                                        aGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGc                                                                                               >  6:64387/1‑146 (MQ=255)                                                        aGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGc                                                                                               >  1:57546/1‑146 (MQ=255)                                                            aGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAAt                                                                                         <  3:299037/148‑1 (MQ=255)                                                            aGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATc                                                                                        >  1:294377/1‑149 (MQ=255)                                                                gccgcTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCg                                                                                       >  4:109861/1‑146 (MQ=255)                                                                gccgcTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGaa                                                                                     <  3:109861/148‑1 (MQ=255)                                                                gccgcTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGaa                                                                                     <  2:57546/148‑1 (MQ=255)                                                                gccgcTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGaa                                                                                     <  8:103208/148‑1 (MQ=255)                                                                gccgcTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGaa                                                                                     <  8:331239/148‑1 (MQ=255)                                                                     accaCCTCCGGCGCTGCGCGCCTACCCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGCCGCACCTTC‑AGCCGCGACACGTCGGGCCCGGCCAACGCCCCACCCCCGCGCCCTGGACATCGACAAGCCAGGGCGACAAACGAAGGccac                                                                               >  2:311901/4‑147 (MQ=255)                                                                       gACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGcctcct                                                                             >  2:83681/1‑149 (MQ=255)                                                                       gACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGcctcc                                                                              <  4:5624/148‑1 (MQ=255)                                                                          cTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGa                                                                          <  6:300506/149‑1 (MQ=255)                                                                              ggCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACTCGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGa                                                                       >  7:65391/1‑148 (MQ=255)                                                                              ggCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACTCGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGa                                                                       >  7:70205/1‑148 (MQ=255)                                                                                   tGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTCGAGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGc                                                                   >  1:265868/1‑148 (MQ=255)                                                                                           tCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGccac                                                           >  4:89445/1‑147 (MQ=255)                                                                                             cGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAgg                                                       >  6:42848/1‑149 (MQ=255)                                                                                              gCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGcc                                                             >  2:38964/1‑142 (MQ=255)                                                                                              gCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGcc                                                             <  1:38964/142‑1 (MQ=255)                                                                                              gCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAggg                                                      <  5:216927/149‑1 (MQ=255)                                                                                              gCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAggg                                                      <  5:42848/149‑1 (MQ=255)                                                                                              gCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAggg                                                      <  1:269925/149‑1 (MQ=255)                                                                                               cTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGa                                                                          >  1:189391/1‑128 (MQ=255)                                                                                               cTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGa                                                                          <  2:189391/128‑1 (MQ=255)                                                                                                  tttCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAGGGGc                                                    <  1:83681/147‑1 (MQ=255)                                                                                                               aGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGc                                                              <  6:246649/124‑1 (MQ=255)                                                                                                               aGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGc                                                              >  5:246649/1‑124 (MQ=255)                                                                                                                      cTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGc                              >  8:145802/1‑149 (MQ=255)                                                                                                                         ccGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGtt                           >  8:9415/1‑149 (MQ=255)                                                                                                                             aGGCGCTCCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCGAGCGCGCTGGTCATCGTCTCGCCAGGGCGAGAATCGAAGGCCTGCTCGAGGTTAGCCATGCCACCAGGGGGGTACACCCCGATGAACATCTGCGTTTTGc                       <  1:311901/149‑1 (MQ=255)                                                                                                                             aGGCGCACCTTC‑AGCGGCGACTCGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGc                       <  8:65391/149‑1 (MQ=255)                                                                                                                             aGGCGCACCTTC‑AGCGGCGACTCGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGc                       <  8:70205/149‑1 (MQ=255)                                                                                                                                     cTTC‑AGCGGCGACACGTCGGGCCCGGCCAACTCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCcgc               <  3:89445/149‑1 (MQ=255)                                                                                                                                     cTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAACCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCcgc               >  5:216603/1‑149 (MQ=255)                                                                                                                                                   acGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCt  <  2:294377/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     ATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTC‑AGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGATAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCT  >  NC_002947/4366341‑4366634

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