Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,425,655	C→T	8.0%	intergenic (+57/‑448)	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,655	0	C	T	8.0%	174.8 / 4.2	63	intergenic (+57/‑448)	PP_1249/PP_1250	DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base T (3/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGT  >  NC_002947/1425508‑1425800                                                                                                                                                    |                                                                                                                                                    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:251232/149‑1 (MQ=255)          cTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCgg                                                                                                                                          <  2:246733/149‑1 (MQ=255)           tACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCggg                                                                                                                                         >  8:97161/1‑149 (MQ=255)           tACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCggg                                                                                                                                         >  8:63161/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                                                                                                                            >  2:251076/1‑149 (MQ=255)                       ccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTt‑gtgtg                                                                                                                            >  1:40675/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                                                                                                           >  8:9028/1‑149 (MQ=255)                                        cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGt                                                                                                           >  3:9861/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                                                                                                                                   <  3:64296/110‑5 (MQ=38)                                                        tGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTG‑CATGACCgggggcctt                                                                                                                                   >  4:64296/1‑106 (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                                                                     <  1:290470/121‑1 (MQ=255)                                                                                                          tGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAg                                                                     >  2:290470/1‑121 (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)                                                                                                                                               ataTGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGCAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg    <  1:89507/147‑1 (MQ=255)                                                                                                                                                 aTGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGATAGTCCTCTAAACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgt  <  1:251076/149‑1 (MQ=255)                                                                                                                                                 aTGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgt  <  6:15689/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGT  >  NC_002947/1425508‑1425800

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