Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,425,650	T→C	8.3%	intergenic (+52/‑453)	PP_1249 → / → PP_1250	DUF4223 domain‑containing protein/group II intron‑encoding maturase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,425,650	0	T	C	8.3%	166.2 / 4.8	60	intergenic (+52/‑453)	PP_1249/PP_1250	DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base C (3/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGG  >  NC_002947/1425508‑1425798                                                                                                                                               |                                                                                                                                                       cACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGc                                                                                                                                                   <  5:33157/148‑1 (MQ=255)  aCCTACGACTTCCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGGTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGca‑c                                                                                                                                                <  5:109585/149‑1 (MQ=255)    cTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCAGGCA‑CAt                                                                                                                                              >  3:218259/1‑149 (MQ=255)          cTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCgg                                                                                                                                        <  2:246733/149‑1 (MQ=255)          cTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCgg                                                                                                                                        <  2:251232/149‑1 (MQ=255)           tACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCggg                                                                                                                                       >  8:63161/1‑149 (MQ=255)           tACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCggg                                                                                                                                       >  8:97161/1‑149 (MQ=255)                   ccACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGttttt                                                                                                                                <  6:76709/148‑1 (MQ=255)                    cACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCTATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTt‑gt                                                                                                                             <  1:305663/149‑1 (MQ=255)                    cACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTt‑gt                                                                                                                             >  3:53442/1‑149 (MQ=255)                     aCCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTt‑gtg                                                                                                                            >  7:301072/1‑149 (MQ=255)                       ccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTt‑gtgtg                                                                                                                          >  1:40675/1‑149 (MQ=255)                       ccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTt‑gtgtg                                                                                                                          >  2:251076/1‑149 (MQ=255)                       ccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGT‑ATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTgc                                                                                                                         >  2:259527/1‑149 (MQ=255)                        cTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTG‑CATGACCgggggccgtgga                                                                                                                              >  8:118324/1‑138 (MQ=255)                         ttCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGAACCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTgcg                                                                                                                        >  5:308302/1‑149 (MQ=255)                             gTTTCCATCTCCAACATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACCATGGCCGGGGGTTTTT‑GTGTGCGCCa                                                                                                                     <  1:176000/149‑1 (MQ=255)                              tttCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAgg                                                                                                                   >  3:324891/1‑149 (MQ=255)                               ttCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAgg                                                                                                                   <  8:301072/148‑1 (MQ=255)                                  cATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATg                                                                                                               >  5:15689/1‑149 (MQ=255)                                  aaTCTCCAAGATCATCGGTGTCTGCGGCCCTATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGAGCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATg                                                                                                               <  7:287824/148‑1 (MQ=255)                                   aTCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATgg                                                                                                              >  2:155944/1‑149 (MQ=255)                                      tcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGgcgc                                                                                                           >  5:60352/1‑149 (MQ=255)                                      tcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGgcgc                                                                                                           >  2:194492/1‑149 (MQ=255)                                        cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGt                                                                                                         >  3:9861/1‑149 (MQ=255)                                        cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGt                                                                                                         >  8:9028/1‑149 (MQ=255)                                                  ggTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTAt                                                                                                <  4:9861/148‑1 (MQ=255)                                                  ggTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTCTGCGCCAGGCATGGCGCGTTGCGCGTAt                                                                                                <  1:194492/148‑1 (MQ=255)                                                   gTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATgc                                                                                              <  7:9028/149‑1 (MQ=255)                                                        tGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTG‑CATGACCgggggcctt                                                                                                                                 >  4:64296/1‑106 (MQ=38)                                                        tGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTG‑CATGACCgggggcctt                                                                                                                                 <  3:64296/110‑5 (MQ=38)                                                                  tcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGtt                                                                               <  2:151599/149‑1 (MQ=255)                                                                   cgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtg                                                                              >  8:72426/1‑149 (MQ=255)                                                                    gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgt                                                                             >  5:194519/1‑149 (MQ=255)                                                                       aaaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt                                                                          >  3:323574/1‑149 (MQ=255)                                                                        aaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                         <  5:245566/149‑1 (MQ=255)                                                                        aaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt                                                                          <  4:313735/148‑1 (MQ=255)                                                                        aaCTGCCTCAGCAGCAGT‑ATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTTAATGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                         <  1:259527/149‑1 (MQ=255)                                                                               gcagcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                  <  1:126069/148‑1 (MQ=255)                                                                                 agcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATgcgc                                                                                            <  4:313888/121‑1 (MQ=255)                                                                                  gcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTggg                                                               <  3:91804/149‑1 (MQ=255)                                                                                   cagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                              <  7:97161/149‑1 (MQ=255)                                                                                           aTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGACCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                      <  6:308302/149‑1 (MQ=255)                                                                                                  cGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                               <  2:310891/149‑1 (MQ=255)                                                                                                      aCCCTGATCGCGACGCTTGCGGAGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTcc                                           <  7:134854/149‑1 (MQ=255)                                                                                                        ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GAGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGTAGGTGAAAGTCctc                                         >  5:166413/1‑149 (MQ=255)                                                                                                          tGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAg                                                                   >  2:290470/1‑121 (MQ=255)                                                                                                          tGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAg                                                                   <  1:290470/121‑1 (MQ=255)                                                                                                                    gCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCgg                              >  8:288805/1‑148 (MQ=255)                                                                                                                     cTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCa                            <  2:314645/149‑1 (MQ=255)                                                                                                                        gcggcgTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGTGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGcc                                                                      <  4:232511/104‑1 (MQ=255)                                                                                                                        gcggcgTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGcc                                                                      >  3:232511/1‑104 (MQ=255)                                                                                                                        gcggcgTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                         >  2:83646/1‑149 (MQ=255)                                                                                                                         cggcgTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                         <  1:155944/148‑1 (MQ=255)                                                                                                                           gcgTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgggg                       >  5:44300/1‑148 (MQ=255)                                                                                                                                 tCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                 <  7:72426/148‑1 (MQ=255)                                                                                                                                   cAAGGCCACCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCATGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCATCTACACACCCGGCAAGGGGAAGTGTTa               <  4:53442/148‑1 (MQ=255)                                                                                                                                         ccccGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg        <  4:324891/149‑1 (MQ=255)                                                                                                                                         ccccGGCCATGTG‑CATGACCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg        <  6:138662/149‑1 (MQ=255)                                                                                                                                         ccccGGCCATGTG‑CATGACCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg        >  5:138662/1‑149 (MQ=255)                                                                                                                                               tCATGGA‑CTTGTCCGGGTGTTTTT‑GTGTGCGCCAGGGATGGCGCGTTGCGCGTATGCGCAACCCGCTTGTTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTGCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg  <  2:40675/149‑1 (MQ=11)                                                                                                                                               tCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg  <  7:63161/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                       CACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGG  >  NC_002947/1425508‑1425798

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