Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,352,176	C→T	7.7%	intergenic (+66/+76)	PP_2947 → / ← PP_2948	transcriptional regulator MvaT/GntR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,352,176	0	C	T	7.7%	187.9 / 4.3	65	intergenic (+66/+76)	PP_2947/PP_2948	transcriptional regulator MvaT/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base T (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

AGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAG  >  NC_002947/3352028‑3352323                                                                                                                                                      |                                                                                                                                                    agaCCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGc                                                                                                                                                     >  2:63703/1‑149 (MQ=255)    ccAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAATCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTc                                                                                                                                                  <  5:7727/149‑1 (MQ=255)    ccAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTc                                                                                                                                                  <  4:244916/149‑1 (MQ=255)       aGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGcc                                                                                                                                               >  7:193027/1‑149 (MQ=255)           tGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGcctcct                                                                                                                                           <  8:167675/149‑1 (MQ=255)              cAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGa                                                                                                                                        <  5:239973/149‑1 (MQ=255)               aaCCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAGGCGACGCTGTCGCTTCCCCGaa                                                                                                                                       <  5:118520/149‑1 (MQ=255)               aaCCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGa                                                                                                                                        >  2:220428/1‑148 (MQ=255)                 ccACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAGGCGACGCTGTCGCTTCCTCGaaaa                                                                                                                                     <  8:132982/149‑2 (MQ=255)                  cacaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGaaa                                                                                                                                      >  5:146805/1‑147 (MQ=255)                   acaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACa                                                                                                                                    <  1:63703/149‑1 (MQ=255)                      aGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGc                                                                                                                                >  5:350075/1‑149 (MQ=255)                               aGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATctt                                                                                                                       >  5:240013/1‑147 (MQ=255)                                 gAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGTGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGATACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtgg                                                                                                                     <  8:91558/149‑1 (MQ=255)                                  aaTGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtggg                                                                                                                    <  1:144534/149‑1 (MQ=255)                                      gAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGtgggtggg                                                                                                                >  7:31570/1‑149 (MQ=255)                                      gAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGtgggtggg                                                                                                                >  7:251552/1‑149 (MQ=255)                                        aaGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCg                                                                                                              <  6:217145/149‑1 (MQ=255)                                             cAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAAt                                                                                                          >  6:354943/1‑148 (MQ=255)                                                    gggCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGAc                                                                                                  <  4:34961/149‑1 (MQ=255)                                                      gCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGt                                                                                                <  4:19663/149‑1 (MQ=255)                                                      gCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGt                                                                                                >  1:368740/1‑149 (MQ=255)                                                        agaagaAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGc                                                                                                      >  1:249391/1‑141 (MQ=255)                                                        agaagaAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGc                                                                                                      <  2:249391/141‑1 (MQ=255)                                                            gaagTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTaa                                                                                          <  7:290016/149‑1 (MQ=255)                                                             aagTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGAGGAATGAGCTGACGTGGTTaaa                                                                                         >  1:374984/1‑149 (MQ=255)                                                                tGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAg                                                                                      >  7:32293/1‑149 (MQ=255)                                                                  gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGa                                                                                    <  7:70059/149‑1 (MQ=255)                                                                          ggTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCg                                                                                                                                                    >  6:119457/1‑77 (MQ=255)                                                                          ggTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCg                                                                                                                                                    <  5:119457/77‑1 (MQ=255)                                                                              tCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACt                                                                        <  8:251552/149‑1 (MQ=255)                                                                                 tGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAg                                                                     >  6:66969/1‑149 (MQ=255)                                                                                        tgCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACCGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                               <  6:350075/148‑1 (MQ=255)                                                                                              tGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCg                                                        <  2:33663/149‑1 (MQ=255)                                                                                                  gCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                     >  5:69301/1‑148 (MQ=255)                                                                                                  gCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                     >  3:63783/1‑148 (MQ=255)                                                                                                  gCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                     >  1:352025/1‑148 (MQ=255)                                                                                                   cccTCATTGCGCTATCTGAAACGACTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGTATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCa                                                   <  1:354190/149‑1 (MQ=255)                                                                                                    ccTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACa                                                                                       >  7:345171/1‑112 (MQ=255)                                                                                                    ccTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACa                                                                                       <  8:345171/112‑1 (MQ=255)                                                                                                        aTTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACt                                              >  3:70416/1‑149 (MQ=255)                                                                                                           gcgcTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCg                                           >  6:76490/1‑149 (MQ=255)                                                                                                            cgcTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCg                                           >  2:19105/1‑148 (MQ=255)                                                                                                             gcTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGt                                          >  7:356604/1‑148 (MQ=255)                                                                                                                aTCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGAGACAGCGTCGCCTCCTCGAAACGGGGCGTTCATCGTGGGTGTGCGGTATGAGAGGACGTGGTTATACAGTATGTTTCGGCACTAAGCGGTCAGTGTACGTTGGAATACttgatatat                                      >  8:341893/1‑142 (MQ=255)                                                                                                                   tGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAAtttt                                   >  2:100798/1‑149 (MQ=255)                                                                                                                     aaaCGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTgg                                 <  6:146805/149‑1 (MQ=255)                                                                                                                     aaaCGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTg                                  <  4:63783/148‑1 (MQ=255)                                                                                                                          ccTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtgg                                                                                                                     <  4:120698/60‑1 (MQ=255)                                                                                                                          ccTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtgg                                                                                                                     >  3:120698/1‑60 (MQ=255)                                                                                                                          ccTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCa                            <  2:120906/149‑1 (MQ=255)                                                                                                                            ttGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGc                          <  2:374984/149‑1 (MQ=255)                                                                                                                                   tttttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTg                    >  4:366347/1‑149 (MQ=255)                                                                                                                                   tgtttCGAGGAGGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTaa                                                                      <  7:351680/97‑1 (MQ=37)                                                                                                                                   tgtttCGAGGAGGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTaa                                                                      >  8:351680/3‑99 (MQ=37)                                                                                                                                    ttttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGCTGGGCGGAATGAGCTGAAGTGGTTAAAAAGGAAGTTTAGCCAATAAGCGCTCAGTGTCCGTTGAAATACTTCGTAAATTTTGGCATCAGCCGATTgg                   >  8:240880/1‑149 (MQ=255)                                                                                                                                     tttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTg                    <  8:32293/147‑1 (MQ=255)                                                                                                                                       tCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGc                 <  5:76490/148‑1 (MQ=255)                                                                                                                                       tCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGc                 <  5:354943/148‑1 (MQ=255)                                                                                                                                         gaggagGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGa               >  5:127100/1‑148 (MQ=255)                                                                                                                                         gAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGa               >  2:262188/1‑148 (MQ=255)                                                                                                                                         gAGGAAGCGACAGCGTCGCCTACTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACACGAAGTTTCGAAACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTCGCATCAGCCGATTGGGCGa               >  5:137391/1‑148 (MQ=255)                                                                                                                                             aaGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTg          <  7:280939/149‑1 (MQ=255)                                                                                                                                                 gACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAAc      >  3:5156/1‑149 (MQ=255)                                                                                                                                                  aCAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACg     >  4:379000/1‑149 (MQ=255)                                                                                                                                                     gCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACggag  <  1:220428/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGAG  >  NC_002947/3352028‑3352323

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