Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,425,656	A→G	7.9%	intergenic (+58/‑447)	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,656	0	A	G	7.9%	176.4 / 3.8	63	intergenic (+58/‑447)	PP_1249/PP_1250	DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base A (25/33);  new base G (3/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGT  >  NC_002947/1425509‑1425800                                                                                                                                                    |                                                                                                                                                   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                                                                                                                            >  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                                                                                                             >  2:194492/1‑149 (MQ=255)                                     tcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGgcgc                                                                                                             >  5:60352/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                                                                                                                                   <  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          >  5:138662/1‑149 (MQ=255)                                                                                                                                        ccccGGCCATGTG‑CATGACCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg          <  6:138662/149‑1 (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)                                                                                                                                                    |                                                                                                                                                   ACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCA‑CATGGCCGGGGGTTTTT‑GTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGT  >  NC_002947/1425509‑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: ‑