Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,960,099	T→A	6.0%	intergenic (+50/‑78)	bolA → / → PP_1758	DNA‑binding transcriptional dual regulator/rhodanese‑like domain‑containing sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,960,099	0	T	A	6.0%	197.3 / 2.5	67	intergenic (+50/‑78)	bolA/PP_1758	DNA‑binding transcriptional dual regulator/rhodanese‑like domain‑containing sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base A (2/2);  total (33/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 = 9.11e-01

ATCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTG  >  NC_002947/1959954‑1960242                                                                                                                                                   |                                                                                                                                                aTCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGc                                                                                                                                              <  7:7382/149‑1 (MQ=255)     aTGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGcc                                                                                                                                          >  5:173936/1‑149 (MQ=255)      tGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGccc                                                                                                                                         <  5:77222/149‑1 (MQ=255)        cccTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCt                                                                                                                                       >  5:139765/1‑149 (MQ=255)        cccTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCt                                                                                                                                       >  8:192163/1‑149 (MQ=255)         ccTGGCCCTCCCCACCTACACCGCCGAAGAGTGGGCGAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCtt                                                                                                                                      >  4:43627/1‑149 (MQ=255)          cTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTg                                                                                                                                     >  7:34052/1‑149 (MQ=255)             gCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGGA‑GAATCCGCATCACCCAAAGGGCGCCGCTTAGCGGCCCCTTGCtg                                                                                                                                  >  8:66761/1‑148 (MQ=255)              ccATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAg                                                                                                                                               >  6:165354/1‑135 (MQ=255)              ccATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAg                                                                                                                                               <  5:165354/135‑1 (MQ=255)                  ccACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTAt                                                                                                                             >  5:49588/1‑149 (MQ=255)                   cacaCCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATg                                                                                                                            >  4:272018/1‑149 (MQ=255)                   cacaCCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAAGGGGCCGCTAAGCGGccccctttgg                                                                                                                                  <  4:202529/143‑7 (MQ=255)                   cacaCCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAAGGGGCCGCTAAGCGGccccctttgg                                                                                                                                  >  3:202529/1‑137 (MQ=255)                    acacCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGc                                                                                                                           <  6:139765/149‑1 (MQ=255)                        cTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTATTGGCCCCTTGCTTTTTATGCAAAc                                                                                                                       <  7:51313/149‑1 (MQ=255)                        cTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAAc                                                                                                                       <  8:39718/149‑1 (MQ=255)                            acCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTc                                                                                                                    >  5:336839/1‑148 (MQ=255)                              cgccgAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCaaa                                                                                                                 >  1:72619/1‑149 (MQ=255)                                ccgAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAAcc                                                                                                               <  6:49588/149‑1 (MQ=255)                                          ggCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGccc                                                                                                     >  4:53208/1‑149 (MQ=255)                                          ggCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGccc                                                                                                     <  3:272018/149‑1 (MQ=255)                                          ggCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGcc                                                                                                      <  6:336839/148‑1 (MQ=255)                                           gCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTACGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCttttt                                                                                                                               <  1:75795/123‑1 (MQ=255)                                           gCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCttttt                                                                                                                               >  2:75795/1‑122 (MQ=255)                                            ccAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGc                                                                                                                                              >  8:2710/1‑106 (MQ=255)                                            ccAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGc                                                                                                                                              <  7:2710/106‑1 (MQ=255)                                             cAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCGCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCt                                                                                                   >  5:309510/1‑148 (MQ=255)                                                 gTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCAGGTCCGCCCCTTAc                                                                                                >  3:220277/1‑147 (MQ=255)                                                  tCGGCGTTGCACCGGCCTCGCCCGGGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAACGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAg                                                                                             <  7:66761/149‑1 (MQ=255)                                                           cACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAaccac                                                                                    >  5:228422/1‑149 (MQ=255)                                                              cGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTg                                                                                 >  1:16510/1‑149 (MQ=255)                                                                   tCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGAt                                                                            >  1:308582/1‑149 (MQ=255)                                                                   tCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTATATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGAt                                                                            >  2:232101/1‑149 (MQ=255)                                                                    cGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGAtt                                                                           >  3:332803/1‑149 (MQ=255)                                                                     gCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGAtt                                                                           <  6:173936/148‑1 (MQ=255)                                                                     gCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTATATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGAt                                                                            <  1:232101/147‑1 (MQ=255)                                                                      cccGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAACGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTcc                                                                         <  3:43627/149‑1 (MQ=255)                                                                          tgtgtgCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTa                                                                     <  8:34052/149‑1 (MQ=255)                                                                           gtgtgCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTAc                                                                    <  2:16510/149‑1 (MQ=255)                                                                            tgtgCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTAc                                                                    <  7:192163/148‑1 (MQ=255)                                                                             gtgCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACAt                                                                  <  6:228422/149‑1 (MQ=255)                                                                                    cggcggGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCaa                                                           <  1:338699/149‑1 (MQ=255)                                                                                          gCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCg                                                                        <  1:338399/130‑1 (MQ=255)                                                                                          gCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCAtc                                                     <  4:332803/149‑1 (MQ=255)                                                                                                   ccGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCg                                                                        <  2:273164/121‑1 (MQ=255)                                                                                                   ccGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCg                                                                        >  1:273164/1‑121 (MQ=255)                                                                                                             gTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTAc                                   >  2:230504/1‑148 (MQ=255)                                                                                                                 tGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTAAGAGGGCAACCACCTGGGGATTCCTTTACATTTATCAAGCAATCGTCTTCTCGGAGATGTAAAAGt                               >  8:317459/1‑148 (MQ=255)                                                                                                                   gTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTcgc                                             >  6:47544/1‑132 (MQ=255)                                                                                                                   gTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTcgc                                             <  5:47544/132‑1 (MQ=255)                                                                                                                   gTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCg                            >  5:58289/1‑149 (MQ=255)                                                                                                                         aTCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCAcc                       <  6:58289/149‑1 (MQ=255)                                                                                                                               aTCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGaa                 <  2:72619/149‑1 (MQ=255)                                                                                                                                  aCCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTAc                                                                    >  1:290975/1‑95 (MQ=255)                                                                                                                                  aCCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTAc                                                                    <  2:290975/95‑1 (MQ=255)                                                                                                                                    ccAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGt                                                                       <  4:235362/90‑1 (MQ=255)                                                                                                                                    ccAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGt                                                                       >  3:235362/1‑90 (MQ=255)                                                                                                                                     cAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTAc           <  6:309510/149‑1 (MQ=255)                                                                                                                                     cAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTAc           >  3:207869/1‑149 (MQ=255)                                                                                                                                      aaaGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACg          <  8:264799/149‑1 (MQ=255)                                                                                                                                        aGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTc        <  2:308582/149‑1 (MQ=255)                                                                                                                                          ggggCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTgg                   <  1:147213/136‑1 (MQ=255)                                                                                                                                          ggggCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTgg                   >  2:147213/1‑136 (MQ=255)                                                                                                                                          ggggCCGCTAAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTAc           >  6:123931/1‑144 (MQ=255)                                                                                                                                          ggggCCGCTAAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTAc           <  5:123931/144‑1 (MQ=255)                                                                                                                                              ccGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTg  >  1:73339/1‑149 (MQ=255)                                                                                                                                                  ttAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGaa                 >  2:227415/1‑130 (MQ=255)                                                                                                                                                  ttAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGaa                 <  1:227415/130‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                ATCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTA‑GAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTG  >  NC_002947/1959954‑1960242

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