Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,505,443	A→G	6.3%	T15A (ACT→GCT)	PP_2199 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,505,443	0	A	G	6.3%	181.6 / 2.4	63	T15A (ACT→GCT)	PP_2199	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/26);  new base G (2/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-01

TTACCTGTATGTACTGACCGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCGGCAAGACCTGGCTGATG  >  NC_002947/2505296‑2505583                                                                                                                                                    |                                                                                                                                             ttACCTGTATGTACTGACCGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCa                                                                                                                                              >  6:24240/1‑148 (MQ=255)   aCCTGTATGTACTGACCGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTc                                                                                                                                           >  8:273685/1‑149 (MQ=255)      tgtatgtaCTGACCGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCggg                                                                                                                                        >  5:83209/1‑149 (MQ=255)           gtaCTGACCGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTAc                                                                                                                                   >  3:151770/1‑149 (MQ=255)            taCTGACCGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCCGAGTGCTCAGGGCttgc                                                                                                                                  <  3:122072/149‑4 (MQ=255)                gACCGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGc                                                                                                                              >  5:169630/1‑149 (MQ=255)                gACCGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGc                                                                                                                              >  3:11505/1‑149 (MQ=255)                   cGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGa                                                                                                                           >  3:189517/1‑149 (MQ=255)                   cGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCCGAGTGCTCAGGGCttgc                                                                                                                                  >  5:177231/1‑139 (MQ=255)                   cGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCCGAGTGCTCAGGGCttgc                                                                                                                                  <  6:177231/142‑4 (MQ=255)                    gATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGAc                                                                                                                          >  2:375808/1‑149 (MQ=255)                     aTGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGAcc                                                                                                                         <  6:136560/149‑1 (MQ=255)                        aCAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCCAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTg                                                                                                                      >  7:63403/1‑149 (MQ=255)                           aGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGcc                                                                                                                    >  2:81329/1‑148 (MQ=255)                            ggATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGcc                                                                                                                    >  5:352474/1‑147 (MQ=255)                             gATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAg                                                                                                                 >  7:327959/1‑149 (MQ=255)                               tGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGGCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGc                                                                                                               >  5:246352/1‑149 (MQ=255)                                 gcgcgcTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGccc                                                                                                             >  5:23639/1‑149 (MQ=255)                                   gcgcTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCcgc                                                                                                           >  7:70349/1‑149 (MQ=255)                                         cTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCgtcgtc                                                                                                     >  3:276521/1‑149 (MQ=255)                                             cagTCGGCCAGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAgg                                                                                                 >  3:175031/3‑149 (MQ=255)                                               gTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGcgtc                                                                                                        <  1:315516/140‑1 (MQ=255)                                               gTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGcgtc                                                                                                        >  2:315516/1‑140 (MQ=255)                                                 cGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCtt                                                                                             <  4:276521/149‑1 (MQ=255)                                                  ggCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTg                                                                                            <  6:23639/149‑1 (MQ=255)                                                    ccTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGCGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGc                                                                                           >  6:81515/1‑148 (MQ=255)                                                           gACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGc                                                                                   <  5:24240/149‑1 (MQ=255)                                                            aCTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGG‑CTTGCAAggcctcga                                                                                 <  6:83209/149‑8 (MQ=255)                                                              tGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTc                                                                                >  6:117881/1‑149 (MQ=255)                                                              tGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTc                                                                                >  6:152625/1‑149 (MQ=255)                                                                     cgcCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTccc                                                                         <  5:117881/149‑1 (MQ=255)                                                                             ccGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAg                                                                 <  6:24656/149‑1 (MQ=255)                                                                              cGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAgg                                                                <  1:375808/149‑1 (MQ=255)                                                                                 gggTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGgcc                                                             <  7:273685/149‑1 (MQ=255)                                                                                    tGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGgccgcc                                                          >  2:223197/1‑149 (MQ=255)                                                                                       cAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCccca                                                       >  6:238225/1‑149 (MQ=255)                                                                                        aGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCcccac                                                      >  3:366646/1‑149 (MQ=255)                                                                                         gTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCcccacc                                                     <  6:246352/149‑1 (MQ=255)                                                                                               tAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGc                                               <  4:175031/149‑1 (MQ=255)                                                                                                  aaTGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCct                                             >  8:48971/1‑148 (MQ=255)                                                                                                   aTGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCACGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCctct                                           <  3:44034/149‑1 (MQ=255)                                                                                                          aTGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTg                                    <  1:81329/149‑1 (MQ=255)                                                                                                               ccTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGcc                               >  6:84815/1‑149 (MQ=255)                                                                                                               ccTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGcc                               <  5:84815/149‑1 (MQ=255)                                                                                                               ccTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGcc                               >  2:159159/1‑149 (MQ=255)                                                                                                                cTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTTCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCATTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGccc                              >  5:21662/1‑149 (MQ=255)                                                                                                                    aCCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTc                          <  5:152625/149‑1 (MQ=255)                                                                                                                    aCCCCCAAACCCGCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTc                          <  8:63403/149‑1 (MQ=255)                                                                                                                      ccccAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCgg                        <  4:189517/149‑1 (MQ=255)                                                                                                                      ccccAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCgg                        <  8:327959/149‑1 (MQ=255)                                                                                                                       cccAAACCCTCTACCAGCACGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCggg                       <  5:13304/149‑1 (MQ=255)                                                                                                                        ccAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTc                          <  8:14938/145‑1 (MQ=255)                                                                                                                        ccAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTc                          >  7:14938/1‑145 (MQ=255)                                                                                                                           aaCCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGcgcg                   >  7:45710/1‑149 (MQ=255)                                                                                                                            aCCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGTCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCgg                  >  5:152915/1‑149 (MQ=255)                                                                                                                             cccTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCGGc                 <  1:229987/149‑1 (MQ=255)                                                                                                                                  tACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCGGCAAGAc            >  1:149352/1‑149 (MQ=255)                                                                                                                                    ccAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCGGCAAGAc            <  8:45710/147‑1 (MQ=255)                                                                                                                                      aGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCGGCAAGACCTgg        <  4:11505/149‑1 (MQ=255)                                                                                                                                       gCAAGCACTGAGCACTCGGGGCTAAGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGAGTGCTTCCAGTCCCTTGAGCACGGCCGCCCCACCCAGGGGCTCTACCTCTGGGGCCCGGTCGGGCGCGGCAAGAACTGGc       >  4:365897/1‑149 (MQ=255)                                                                                                                                       cgtaGCCCCGAGTGCTCAGGGCTACGTTGCAGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGACCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCGGCAAGACCTGGc       >  1:31179/4‑149 (MQ=255)                                                                                                                                            cccTGAGCACTCGGGGCTACGTTGCTGTCCATGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCTCCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCGGCAAGACCTGGCTGATg  <  6:169630/149‑1 (MQ=255)                                                                                                                                            cccTGAGCACTCGGGGCTACGTTGCTGACCCTTCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTAAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCGGCAAGACCTGGCTGATg  <  8:70349/149‑1 (MQ=255)                                                                                                                                            cccTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGACCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCGGCAAGACCTGGCTGATg  >  1:42979/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                             TTACCTGTATGTACTGACCGATGACAAGGATGGCGCGCTACTCAAGGTCGGCCTGGCCGACTGACGCGCGCCTGCACCGAGGGTGACAGTCGGGTAGAATGCGGCATGGCCCTCGACCCCCAAACCCTCTACCAGCAAGCCCTGAGCACTCGGGGCTACGTTGCTGACCCTGCCCAGGCCCGCGTCGTCGAGGCCTTGCAAGCCTGCTTCCAGTCCCTTGAGCAGGGCCGCCCCACCCAGGGCCTCTACCTCTGGGGCCCGGTCGGGCGCGGCAAGACCTGGCTGATG  >  NC_002947/2505296‑2505583

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