Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,642,077	G→A	6.5%	S32S (AGC→AGT)	PP_1442 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,642,077	0	G	A	6.5%	183.6 / 2.1	61	S32S (AGC→AGT)	PP_1442	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base A (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

CAGGCGCGCTGTAAACCTCGGGGCCAAGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAA  >  NC_002947/1641935‑1642223                                                                                                                                                |                                                                                                                                                   cAGGCGCGCTGTAAACCTCGGGGCCAAGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTc                                                                                                                                               <  5:194964/148‑1 (MQ=255)  aGGCGCGCTGTAAACCTCGGGGCCAAGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCg                                                                                                                                              <  2:35256/148‑1 (MQ=255)     cgcgcTGTAAACCTCGGGGCCAAGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGt                                                                                                                                          >  3:243034/1‑149 (MQ=255)         cTGTAAACCTCGGCGCCAAGGCTATCCAAAACGTTGCGCGCCGCGTTCTCTACAAGCCACCGGTAGCATGAGGTGGGGTTCC‑TGGGCAGGGTAAGCGTGAGGATCTGGACGACCTGCAGGCGCTTGGGGCATTCGATTTCGGCGACCAg                                                                                                                                      >  8:192971/1‑149 (MQ=255)            tAAACCTCGGGGCCAAGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGt                                                                                                                                    >  2:69292/1‑148 (MQ=255)                  tCGGGGCCAAGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACggg                                                                                                                             >  1:174611/1‑149 (MQ=255)                  tCGGGGCCAAGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCACTTTCGGCGTCCAGGTCAACggg                                                                                                                             >  7:17954/1‑149 (MQ=255)                   cGGGGCCAAGGCTTTGCAATAAGTAGGGAGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGAGGTGGGGTTGC‑TGGGCAGGGTAACCGTGAGGATCTGGACGACCTGCACGCGCTTGGGGCATTCGCTTTCGGCGTCCAGTTCAACGGGc                                                                                                                            >  7:102165/1‑149 (MQ=255)                    ggggCCAAGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCt                                                                                                                           <  3:144524/149‑1 (MQ=255)                        ccAAGGCTTTGCAATACCTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGGGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGa                                                                                                                       <  7:254193/149‑1 (MQ=255)                         cAAGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAg                                                                                                                      >  5:269289/1‑149 (MQ=255)                          aaGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAgg                                                                                                                     <  5:120802/149‑1 (MQ=255)                             gCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGATCAACGGGCTGCGAGGAtt                                                                                                                  <  3:220522/149‑1 (MQ=255)                                    aaTACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACg                                                                                                                               <  3:159513/129‑1 (MQ=255)                                    aaTACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACg                                                                                                                               >  4:159513/1‑129 (MQ=255)                                         gTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGCGTGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTgcgcgc                                                                                                         >  1:66990/1‑147 (MQ=255)                                              gcgcCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGa                                                                                                  >  3:56550/1‑148 (MQ=255)                                                     gTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAggg                                                                                          <  7:192971/149‑1 (MQ=255)                                                     gTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCACTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAggg                                                                                          <  8:17954/149‑1 (MQ=255)                                                      ttCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGc                                                                                                     >  4:247881/1‑137 (MQ=255)                                                      ttCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGc                                                                                                     <  3:247881/137‑1 (MQ=255)                                                          gTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAg                                                                                     >  8:129221/1‑149 (MQ=255)                                                              cagccagcGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAg                                                                                     >  6:207934/1‑145 (MQ=255)                                                              cagccagcGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAg                                                                                     <  5:207934/145‑1 (MQ=255)                                                              cagccagcGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCt                                                                                   >  1:146587/1‑147 (MQ=255)                                                                     cGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAAc                                                                          <  7:25637/149‑1 (MQ=255)                                                                       gTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTgcgcg                                                                                                          <  7:20490/115‑1 (MQ=255)                                                                       gTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTgcgcg                                                                                                          >  8:20490/1‑115 (MQ=255)                                                                       gTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGa                                                                        <  7:129221/149‑1 (MQ=255)                                                                       gTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGa                                                                        <  5:255267/149‑1 (MQ=255)                                                                          gCCGGTGGTGGGGTTTC‑TGGGCAGGGTAAGCGTGTGGATCTGGGCGACTTGCAGGCGCGTGGGGCATTGGCTTTCGGCGTCCAGGTCAACGGGCTGGGAGGATTTTTGCGCGCTAGCAGAAAGCAGGGCGAAGCTCAGGGGATCGAGc                                                                      <  8:102165/148‑1 (MQ=255)                                                                           ccGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCTTCCAGGTCAAGGGGCTGCGATGATTGTTGCGCGCAAGCAGAAAGCAGGGGGAAGCTCAGGGGAACGAGCa                                                                     <  4:243034/148‑1 (MQ=255)                                                                           ccGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCACTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCa                                                                                                    >  7:108615/1‑117 (MQ=255)                                                                           ccGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCACTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCa                                                                                                    <  8:108615/117‑1 (MQ=255)                                                                            cGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGa                                                                   >  8:148653/1‑149 (MQ=255)                                                                            cGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGa                                                                   >  8:233954/1‑149 (MQ=255)                                                                            cGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGa                                                                   >  8:32957/1‑149 (MQ=255)                                                                               tggtggGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACga                                                                <  2:66990/149‑1 (MQ=255)                                                                                 gtggGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACgagg                                                              <  2:174611/149‑1 (MQ=255)                                                                                     ggTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGagag                                                          <  7:233954/149‑1 (MQ=255)                                                                                      gTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAgg                                                         >  2:190350/1‑149 (MQ=255)                                                                                       ttGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAgggg                                                                             <  1:21320/128‑1 (MQ=255)                                                                                       ttGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAgggg                                                                             >  2:21320/1‑128 (MQ=255)                                                                                             gggCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCAtgg                                                   >  7:108264/1‑149 (MQ=255)                                                                                              ggCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCAtg                                                    <  2:146587/147‑1 (MQ=255)                                                                                                 aGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGgaga                                                           <  1:114199/137‑1 (MQ=255)                                                                                                 aGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGgaga                                                           >  2:114199/1‑137 (MQ=255)                                                                                                           tGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATCGTGGCTCCGAAGGCt                                     >  5:209004/1‑149 (MQ=255)                                                                                                                aTCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCg                                 >  5:201624/1‑148 (MQ=255)                                                                                                                aTCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCg                                 <  6:201624/148‑1 (MQ=255)                                                                                                                    ggCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAg                                        >  3:232283/1‑137 (MQ=255)                                                                                                                    ggCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAg                                        <  4:232283/137‑1 (MQ=255)                                                                                                                         aCCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAg                       <  7:148653/149‑1 (MQ=255)                                                                                                                            tGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAAc                                                                          <  2:71312/95‑1 (MQ=255)                                                                                                                            tGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAAc                                                                          >  1:71312/1‑95 (MQ=255)                                                                                                                                     ttGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCt                                     >  8:124850/1‑123 (MQ=255)                                                                                                                                       ggggCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATc         >  1:137234/1‑149 (MQ=255)                                                                                                                                        gggCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATc         <  1:190350/148‑1 (MQ=255)                                                                                                                                          gCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAg      >  1:119139/1‑149 (MQ=255)                                                                                                                                          gCATTCGCGTTCCGCGTCCAGGGCAACGGGCTGCGAGGCTTTTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAg      <  8:222081/149‑1 (MQ=255)                                                                                                                                              tCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGaaa  >  3:184844/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                   CAGGCGCGCTGTAAACCTCGGGGCCAAGGCTTTGCAATACGTTGGGCGCCGGGTTCTGTACCAGCCAGCGGTAGCCGGTGGTGGGGTTGC‑TGGGCAGGGTAAGCGTGAGGATCTGGCCGACCTGCAGGCGCTTGGGGCATTCGCTTTCGGCGTCCAGGTCAACGGGCTGCGAGGATTGTTGCGCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAA  >  NC_002947/1641935‑1642223

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