Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,244,394	G→A	6.4%	P31S (CCG→TCG)	rnt ←	ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,244,394	0	G	A	6.4%	181.8 / 2.1	63	P31S (CCG→TCG)	rnt	ribonuclease T
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base A (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTCGCCAGCGCGCAGCGGGACTACGG  >  NC_002947/1244248‑1244542                                                                                                                                                     |                                                                                                                                                     gcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCa                                                                                                                                                    >  3:235264/1‑149 (MQ=255)         aCGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCcg                                                                                                                                            >  1:309821/1‑149 (MQ=255)           ggCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCATTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCcgcg                                                                                                                                          >  1:74394/1‑149 (MQ=255)             cTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGCCCGGCAGATAGCCGCgg                                                                                                                                         >  6:167109/1‑148 (MQ=255)                 aCTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCa                                                                                                                                                    <  6:159471/133‑1 (MQ=255)                 aCTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCa                                                                                                                                                    >  5:159471/1‑133 (MQ=255)                   tCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCt                                                                                                                                  >  8:259571/1‑149 (MQ=255)                     ggTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTc                                                                                                                                 <  5:70457/148‑1 (MQ=255)                                gTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTg                                                                                                                     <  4:235264/149‑1 (MQ=255)                                 tATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTgg                                                                                                                    >  6:256143/1‑149 (MQ=255)                                          ggAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCg                                                                                                           >  6:244460/1‑149 (MQ=255)                                            aaCAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCgct                                                                                                         <  7:259571/149‑1 (MQ=255)                                               aGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCgctgct                                                                                                      >  5:15763/1‑149 (MQ=255)                                               aGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCgctgct                                                                                                      >  2:78634/1‑149 (MQ=255)                                                ggAAACCTTTCTCGTCCATGCCAAAGGAGACGGCGGCAAATTACAGCAATGCATCGGTAGCGCTGTTGAAACCACAGGTTTCCACCTCC‑‑ACCACGCCCGCCAGATAGCCCCCGAAGCCCTCGGCCATCGGGGGGCGAGGAACgctgcag                                                                                                     >  2:289617/1‑147 (MQ=255)                                                    aCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAg                                                                                                                <  4:157223/134‑1 (MQ=255)                                                    aCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAg                                                                                                                >  3:157223/1‑134 (MQ=255)                                                    aCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCtt                                                                                                 <  5:30970/149‑1 (MQ=255)                                                    aCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCtt                                                                                                 <  6:15763/149‑1 (MQ=255)                                                      cTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTTCTGCCTTgg                                                                                               <  2:249944/149‑1 (MQ=255)                                                      cTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                                               <  7:5317/149‑1 (MQ=255)                                                      cTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                                               <  5:256143/149‑1 (MQ=255)                                                       tttCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTggg                                                                                              >  1:164985/1‑149 (MQ=255)                                                             gTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTc                                                                                          <  5:244460/147‑1 (MQ=255)                                                                 aTGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtc                                                                                    <  6:305331/149‑1 (MQ=255)                                                                     cAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGTCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCt                                                                                <  1:78634/149‑1 (MQ=255)                                                                       aTGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCGGTTGAAACCACCGGTTTCAACTTCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                                              <  5:167109/149‑1 (MQ=255)                                                                             aCGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAg                                                                        <  2:74394/149‑1 (MQ=255)                                                                                  ggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCg                                                                                                           >  2:250227/1‑109 (MQ=255)                                                                                  ggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCg                                                                                                           <  1:250227/109‑1 (MQ=255)                                                                                  ggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCtt                                                                   >  4:167357/1‑149 (MQ=255)                                                                                  ggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCt                                                                    >  8:226931/1‑148 (MQ=255)                                                                                          ccAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCAcg                                                           >  3:231629/1‑149 (MQ=255)                                                                                           cAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                                               <  6:279504/112‑1 (MQ=255)                                                                                           cAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                                               >  5:279504/1‑112 (MQ=255)                                                                                            aGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGACAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCAcgcg                                                         >  8:142802/1‑149 (MQ=255)                                                                                            aGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGACAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCAcgcg                                                         >  8:147506/1‑149 (MQ=255)                                                                                             gCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGt                                                        >  7:225357/1‑149 (MQ=255)                                                                                              cAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGAAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGACTTCGTAAAGGTCTTCGCTCACGCGtt                                                       >  8:302829/1‑149 (MQ=255)                                                                                                attcATCGGGAGCGGTTGTGAAACCAACGGTTTCAATATGCTGTCCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGtt                                                       <  1:289617/146‑1 (MQ=255)                                                                                                    aTCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTtctc                                                 >  1:206480/1‑149 (MQ=255)                                                                                                    aTCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTtctc                                                 >  5:205347/1‑149 (MQ=255)                                                                                                            ctcTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTtct                                                  <  5:343085/138‑1 (MQ=255)                                                                                                            cgcTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTtct                                                  >  6:343085/1‑140 (MQ=255)                                                                                                                 ttGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGACAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCa                                     <  7:142802/148‑1 (MQ=255)                                                                                                                 ttGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGACAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCa                                     <  7:147506/148‑1 (MQ=255)                                                                                                                    aaaCCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAgcg                                  <  8:225357/148‑1 (MQ=255)                                                                                                                        caccGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCca                                               >  1:349771/1‑131 (MQ=255)                                                                                                                        caccGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCca                                               <  2:349771/131‑1 (MQ=255)                                                                                                                        caccGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAgg                                          >  4:109714/1‑136 (MQ=255)                                                                                                                        caccGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAgg                                          <  3:109714/136‑1 (MQ=255)                                                                                                                            ggTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTc                         <  6:205347/149‑1 (MQ=255)                                                                                                                            ggTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTc                         <  3:167357/149‑1 (MQ=255)                                                                                                                            ggTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTc                         <  7:226931/149‑1 (MQ=255)                                                                                                                                  aaCATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCACCAGCGCAGTTTTTCGCCAgc                   >  8:130460/1‑149 (MQ=255)                                                                                                                                  aaCATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCACCAGCGCAGTTTTTCGCCAgc                   >  8:110359/1‑149 (MQ=255)                                                                                                                                           accacGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTCGCCAGCGCGCAGc            >  5:287452/1‑149 (MQ=255)                                                                                                                                           accacGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTCGCCAGCGCGCAGc            >  2:343836/1‑149 (MQ=255)                                                                                                                                            ccacGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTCGCCAGCGCGCAGCg           >  2:287389/1‑149 (MQ=255)                                                                                                                                             cacGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTtctc                                                 <  5:228862/110‑1 (MQ=255)                                                                                                                                             cacGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTtctc                                                 >  6:228862/1‑110 (MQ=255)                                                                                                                                             cacGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTc                         <  6:152329/134‑1 (MQ=255)                                                                                                                                             cacGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTCGCCAGCGCGCAGCgg          <  6:287452/149‑1 (MQ=255)                                                                                                                                                     gCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTCGCCAgcgcg                >  8:154334/1‑135 (MQ=255)                                                                                                                                                     gCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTCGCCAgcgcg                <  7:154334/135‑1 (MQ=255)                                                                                                                                                     gCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTCGCCAGCGCGCAGCGGGACTACgg  <  2:164985/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCC‑‑ACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTTCGCCAGCGCGCAGCGGGACTACGG  >  NC_002947/1244248‑1244542

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