Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,352,174	A→G	7.7%	intergenic (+64/+78)	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,174	0	A	G	7.7%	180.0 / 5.0	65	intergenic (+64/+78)	PP_2947/PP_2948	transcriptional regulator MvaT/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base G (2/3);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

AGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACG  >  NC_002947/3352028‑3352320                                                                                                                                                    |                                                                                                                                                   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                                                  >  1:352025/1‑148 (MQ=255)                                                                                                  gCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                  >  3:63783/1‑148 (MQ=255)                                                                                                  gCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                  >  5:69301/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                                                                   >  8:351680/3‑99 (MQ=37)                                                                                                                                   tgtttCGAGGAGGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTaa                                                                   <  7:351680/97‑1 (MQ=37)                                                                                                                                    ttttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGCTGGGCGGAATGAGCTGAAGTGGTTAAAAAGGAAGTTTAGCCAATAAGCGCTCAGTGTCCGTTGAAATACTTCGTAAATTTTGGCATCAGCCGATTgg                >  8:240880/1‑149 (MQ=255)                                                                                                                                     tttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTg                 <  8:32293/147‑1 (MQ=255)                                                                                                                                       tCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGc              <  5:354943/148‑1 (MQ=255)                                                                                                                                       tCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGc              <  5:76490/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)                                                                                                                                                    |                                                                                                                                                   AGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTC‑TTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACG  >  NC_002947/3352028‑3352320

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