Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,831,780	T→C	7.4%	L220P (CTG→CCG)	yfiP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,831,780	0	T	C	7.4%	185.3 / 4.3	67	L220P (CTG→CCG)	yfiP	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base C (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 = 9.42e-01

GCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATT  >  NC_002947/1831639‑1831908                                                                                                                                              |                                                                                                                                  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                                                                                                                    <  5:82607/149‑1 (MQ=255)         cgcgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCa                                                                                                                    <  6:242767/149‑1 (MQ=255)               gatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTg                                                                                                                           >  1:88548/1‑136 (MQ=255)               gatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTg                                                                                                                           <  2:88548/136‑1 (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                                                                            >  4:290522/1‑148 (MQ=255)                                                  ccTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTcc                                                                            >  5:2933/1‑148 (MQ=255)                                                     ggACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgg                                                                        >  6:56518/1‑149 (MQ=255)                                                      gACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAggg                                                                       >  8:31819/1‑149 (MQ=255)                                                      gACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAggg                                                                       >  8:26541/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                                                     <  6:126875/148‑1 (MQ=255)                                                                        cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTg                                                     >  5:126875/1‑148 (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                                                                       >  3:230273/1‑111 (MQ=255)                                                                                            tGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAggg                                                                       <  4:230273/111‑1 (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  <  5:56518/149‑1 (MQ=255)                                                                                                                          ccTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGAtt  <  1:326519/149‑1 (MQ=255)                                                                                                                             gAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgg                                                                        >  4:329213/1‑77 (MQ=255)                                                                                                                             gAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgg                                                                        <  3:329213/77‑1 (MQ=255)                                                                                                                                              |                                                                                                                                  GCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGG‑GCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATT  >  NC_002947/1831639‑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: ‑