Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,831,777	G→A	7.4%	R219Q (CGG→CAG)	yfiP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,831,777	0	G	A	7.4%	195.8 / 3.9	67	R219Q (CGG→CAG)	yfiP	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base A (2/3);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

CCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATT  >  NC_002947/1831631‑1831908                                                                                                                                                   |                                                                                                                                     ccTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCggcggc                                                                                                                                    >  8:135424/1‑149 (MQ=255)         gcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCtt                                                                                                                            <  7:307368/149‑1 (MQ=255)         gcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCtt                                                                                                                            >  3:93376/1‑149 (MQ=255)          cTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAGCCGCCGCTTGCCGCtca                                                                                                                           <  7:184098/149‑3 (MQ=255)          cTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTg                                                                                                                           >  8:276696/1‑149 (MQ=255)          cTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTg                                                                                                                           <  7:276696/149‑1 (MQ=255)          cTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTg                                                                                                                           <  2:275758/149‑1 (MQ=255)           cggtaGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGa                                                                                                                          <  2:149325/145‑1 (MQ=255)             cAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTGAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGAc                                                                                                                         <  6:132017/148‑1 (MQ=255)             cAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTgcc                                                                                                                         <  4:93376/148‑3 (MQ=255)             cAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGAAAGCGGCGGCTGCCGCTTGAc                                                                                                                         >  8:304536/1‑148 (MQ=255)                 cgcgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCa                                                                                                                    <  6:242767/149‑1 (MQ=255)                 cgcgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCa                                                                                                                    <  5:82607/149‑1 (MQ=255)                       gatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTg                                                                                                                           <  2:88548/136‑1 (MQ=255)                       gatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTg                                                                                                                           >  1:88548/1‑136 (MQ=255)                         tCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGaca                                                                                                            >  3:15876/1‑149 (MQ=255)                                       cAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                              >  5:313888/1‑149 (MQ=255)                                         gAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCa                                                                                            <  5:350076/149‑1 (MQ=255)                                           aGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCAt                                                                                           <  6:313888/148‑1 (MQ=255)                                             ttGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGGTTAGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAATCGGGGGCTGCCGCTTGACGAGCAGGACGACACCCACAAGCGTTTGCATAc                                                                                         <  4:68733/148‑1 (MQ=255)                                               gcggcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcacccac                                                                                                       <  8:96888/132‑1 (MQ=255)                                               gcggcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcacccac                                                                                                       >  7:96888/1‑132 (MQ=255)                                                cggcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTc                                                                                     <  7:135424/149‑1 (MQ=255)                                                         gCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCtt                                                                                                                            <  4:6776/101‑1 (MQ=255)                                                         gCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCtt                                                                                                                            >  3:6776/1‑101 (MQ=255)                                                          ccTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTcc                                                                            >  5:2933/1‑148 (MQ=255)                                                          ccTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTcc                                                                            >  4:290522/1‑148 (MQ=255)                                                             ggACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgg                                                                        >  6:56518/1‑149 (MQ=255)                                                              gACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAggg                                                                       >  8:26541/1‑149 (MQ=255)                                                              gACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAggg                                                                       >  8:31819/1‑149 (MQ=255)                                                                     cTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAAc                                                               <  1:102665/149‑1 (MQ=255)                                                                      tGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAAcc                                                              >  1:242597/1‑149 (MQ=255)                                                                          gATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTg                                                          >  1:187088/1‑149 (MQ=255)                                                                              gCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCttt                                                      >  1:33408/1‑149 (MQ=255)                                                                                cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTg                                                     >  5:126875/1‑148 (MQ=255)                                                                                cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTg                                                     <  6:126875/148‑1 (MQ=255)                                                                                   gCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTc                                                 <  7:26541/149‑1 (MQ=255)                                                                                    cTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTc                                                 <  7:31819/148‑1 (MQ=255)                                                                                    cTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCa                                                <  4:15876/149‑1 (MQ=255)                                                                                     ttCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCa                                                >  8:305914/1‑148 (MQ=255)                                                                                       cGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCa                                                <  7:305914/146‑1 (MQ=255)                                                                                       cGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATa                                              <  7:304536/148‑1 (MQ=255)                                                                                          aGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGa                                          <  2:242597/149‑1 (MQ=255)                                                                                                    tGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAggg                                                                       <  4:230273/111‑1 (MQ=255)                                                                                                    tGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAggg                                                                       >  3:230273/1‑111 (MQ=255)                                                                                                     gCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTgg                                 >  2:326519/1‑147 (MQ=255)                                                                                                         aTCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGc                             >  2:280627/1‑147 (MQ=255)                                                                                                          tCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAGCCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGc                             >  2:243317/1‑146 (MQ=255)                                                                                                              gATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGCGCACAGCCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCa                           <  1:243317/145‑1 (MQ=255)                                                                                                              gATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTg                      <  1:280627/149‑1 (MQ=255)                                                                                                               aTGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAgg                         >  8:223235/1‑145 (MQ=255)                                                                                                               aTGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAgg                         <  7:223235/145‑1 (MQ=255)                                                                                                                 gtgtTCAGCCTGCATTACCTGATCGGCAAGCGGCGACTGCCGCTTGACGAGCAGGACGACACACACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGc                   >  2:39120/1‑149 (MQ=255)                                                                                                                  tgtTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGc                   <  3:290522/148‑1 (MQ=255)                                                                                                                  tgtTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGc                   <  2:187088/148‑1 (MQ=255)                                                                                                                  tgtTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGc                   <  2:33408/148‑1 (MQ=255)                                                                                                                    ttCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑Gca                                                                   >  4:31694/1‑98 (MQ=255)                                                                                                                       aGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACg             >  1:8963/1‑149 (MQ=255)                                                                                                                       aGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACg             >  2:283147/1‑149 (MQ=255)                                                                                                                           tGCTTTACCTCCGCGCCATGCGCGGGATGCC‑CTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAgg        <  8:159982/149‑1 (MQ=255)                                                                                                                            gCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAgg        >  7:260039/1‑149 (MQ=255)                                                                                                                             cATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAg         <  1:229021/147‑1 (MQ=255)                                                                                                                              aTTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGc      >  1:322650/1‑149 (MQ=255)                                                                                                                                  ccTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGAtt  >  1:162389/1‑149 (MQ=255)                                                                                                                                  ccTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGAtt  <  1:326519/149‑1 (MQ=255)                                                                                                                                  ccTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGAtt  <  5:56518/149‑1 (MQ=255)                                                                                                                                     gAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgg                                                                        >  4:329213/1‑77 (MQ=255)                                                                                                                                     gAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgg                                                                        <  3:329213/77‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                     CCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATT  >  NC_002947/1831631‑1831908

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