Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,790,625	A→T	6.1%	intergenic (‑29/+28)	PP_5073 ← / ← hemE	hypothetical protein/uroporphyrinogen decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,790,625	0	A	T	6.1%	191.9 / 2.4	66	intergenic (‑29/+28)	PP_5073/hemE	hypothetical protein/uroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base T (2/2);  total (32/34)
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.49e-01

CAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTG  >  NC_002947/5790479‑5790769                                                                                                                                                    |                                                                                                                                                 cAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGc                                                                                                                                                <  1:44439/149‑1 (MQ=255)   aaTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCg                                                                                                                                              >  2:221146/1‑149 (MQ=255)     tGCGCCAGTGCTTAACAATGGGGCGAGGATGAGGGCGATATATCGGATTTTCATGATTATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCgg                                                                                                                                             <  3:38688/148‑1 (MQ=255)       cgcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGtt                                                                                                                                          <  6:366594/149‑1 (MQ=255)           aGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCTCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACt                                                                                                                                      <  2:322178/149‑1 (MQ=255)           aGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACt                                                                                                                                      <  5:187758/149‑1 (MQ=255)             tGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACtgt                                                                                                                                    >  7:141614/1‑149 (MQ=255)             tGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAGCGAAGCCGGGttttttg                                                                                                                                     >  6:269878/1‑144 (MQ=255)              gCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCAGTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACtgt                                                                                                                                    >  1:103268/1‑149 (MQ=255)              gCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCAATGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACtgt                                                                                                                                    >  1:98866/1‑149 (MQ=255)               cTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGttttttgtac                                                                                                                                  <  5:269878/149‑8 (MQ=255)                tGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACtgtg                                                                                                                                   >  6:125467/1‑147 (MQ=255)                tGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgca                                                                                                                                 >  3:150615/1‑149 (MQ=255)                   aCAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagtg                                                                                                                              <  5:125467/149‑1 (MQ=255)                     aaGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagtgca                                                                                                                            <  1:221146/149‑1 (MQ=255)                           gCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGt                                                                                                                      <  2:98866/149‑1 (MQ=255)                           gCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGt                                                                                                                      <  2:103268/149‑1 (MQ=255)                                gATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtg                                                                                                                 >  6:296961/1‑149 (MQ=255)                                     gggCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtgg                                                                                                             >  3:214261/1‑148 (MQ=255)                                      ggCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTa                                                                                                           >  6:232011/1‑149 (MQ=255)                                        cGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCTCCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTAc                                                                                                          >  4:178725/1‑148 (MQ=255)                                             tAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtg                                                                                                              >  8:303283/1‑139 (MQ=255)                                             tAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtg                                                                                                              <  7:303283/139‑1 (MQ=255)                                              aaCGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCg                                                                                                   >  3:22498/1‑149 (MQ=255)                                                cGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTAACTGGCGGTACTGCGCCGa                                                                                                  >  3:262190/1‑148 (MQ=255)                                                      ttCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTc                                                                                            <  5:232011/148‑1 (MQ=255)                                                      ttCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTc                                                                                            <  5:296961/148‑1 (MQ=255)                                                        cATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGACATCGAAGCCGGGTTTACTgtgca                                                                                                                                 <  7:285306/109‑1 (MQ=255)                                                        cATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGACATCGAAGCCGGGTTTACTgtgca                                                                                                                                 >  8:285306/1‑109 (MQ=255)                                                         aTGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGCCGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCCGTGCATCGAGTCACTGGTGGTACTGCGCCGCCAgtgcgtgc                                                                                        >  8:228853/1‑149 (MQ=255)                                                             tAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACg                                                                                     <  4:214261/148‑1 (MQ=255)                                                              aaTCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGc                                                                                   >  1:77393/1‑149 (MQ=255)                                                                  gggTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCTGGTTTACTGTGCAGTGCATCAAGTCACTGGTTGTACTGCGCCGACAGTTCGTGCACGGCGTtg                                                                               <  3:343256/149‑1 (MQ=255)                                                                   ggTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGCTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGAGCCGACAGTTAGTGCACGGCGTtga                                                                              >  7:83526/1‑149 (MQ=255)                                                                        cTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCTTAGCCGGGTTTACTGTGCAGTGCATCATGTCACTGGTGGTACTGCGCCCACAGTTCGTGCACGGCGTtgat                                                                             <  7:228553/146‑1 (MQ=255)                                                                        cTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGaa                                                                          <  8:141614/149‑1 (MQ=255)                                                                         tggatcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAAt                                                                         <  7:30710/145‑1 (MQ=255)                                                                         tGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAAt                                                                         >  8:127524/1‑149 (MQ=255)                                                                           tttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGt                                                                                 >  3:66150/1‑139 (MQ=255)                                                                           tttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGt                                                                                 <  4:66150/139‑1 (MQ=255)                                                                           tttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAAt                                                                         <  6:64917/147‑1 (MQ=255)                                                                           tttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATg                                                                        >  5:201714/1‑148 (MQ=255)                                                                              cggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCAcc                                                                    >  6:256014/1‑149 (MQ=255)                                                                                gtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCg                                                                  <  2:77393/149‑1 (MQ=255)                                                                                         ggACCAGATTAGCGACGTTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTcc                                                          <  7:127524/148‑1 (MQ=255)                                                                                          gACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCgg                                                        >  1:1131/1‑149 (MQ=255)                                                                                                      gacgaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCgg                                            >  3:306120/1‑149 (MQ=255)                                                                                                        cgaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCg                                          <  4:150615/149‑1 (MQ=255)                                                                                                             gCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGAt                                      <  2:1131/148‑1 (MQ=255)                                                                                                                  aTTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTg                                <  4:22498/149‑1 (MQ=255)                                                                                                                      aaaCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGcc                             >  2:94725/1‑148 (MQ=255)                                                                                                                         cTTAACGTACAAAAAAGCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAgg                         <  3:178725/149‑1 (MQ=255)                                                                                                                          ttAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACCGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGt                                 <  8:342954/140‑1 (MQ=255)                                                                                                                          ttAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACCGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGt                                 >  7:342954/1‑140 (MQ=255)                                                                                                                            aaCGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAgg                         <  7:242376/146‑1 (MQ=255)                                                                                                                            aaCGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAgg                         >  8:242376/1‑146 (MQ=255)                                                                                                                              cGTAGAAAAATCGCGGCATCGAAGTCGGGTTTACTGTGCATTGCATCAAGTCATTGGTGGTAATGCTCCGACAGTTCGTGGAAGGCGTTGATGAATGCTCACGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCAAAGGTTTaa                    <  8:83526/149‑1 (MQ=255)                                                                                                                                 aCAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAaca                  >  3:284228/1‑148 (MQ=255)                                                                                                                                  cAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAacac                 <  1:94725/148‑1 (MQ=255)                                                                                                                                     aaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAAcaca                >  6:334366/1‑146 (MQ=255)                                                                                                                                     aaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAAcaca                <  5:334366/146‑1 (MQ=255)                                                                                                                                       aaCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGc            >  5:53685/1‑148 (MQ=255)                                                                                                                                        aCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGc            >  8:58596/1‑147 (MQ=255)                                                                                                                                          ccGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGt        >  3:120388/1‑149 (MQ=255)                                                                                                                                           cGGCATCGACGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGCGGGACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTg       <  4:262190/149‑1 (MQ=255)                                                                                                                                                 cGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTg  <  7:298546/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTC‑CTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTG  >  NC_002947/5790479‑5790769

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