Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	956,585	A→C	7.6%	intergenic (+18/‑96)	alaC → / → PP_0818	aminotransferase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	956,585	0	A	C	7.6%	180.1 / 3.2	61	intergenic (+18/‑96)	alaC/PP_0818	aminotransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/25);  new base C (3/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTCGA  >  NC_002947/956438‑956730                                                                                                                                                     |                                                                                                                                                  ggggATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTCTGGCTCGCAAGTAAGCCCAATGGGGCTGCGaa                                                                                                                                                   <  1:237124/149‑1 (MQ=255)               gTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGcc                                                                                                                                    >  2:216972/1‑149 (MQ=255)                 gACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCctc                                                                                                                                  >  3:271712/1‑149 (MQ=255)                      tGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGGCTCTTggg                                                                                                                             >  7:258919/1‑149 (MQ=255)                      tGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTgg                                                                                                                              <  4:271712/148‑1 (MQ=255)                      tGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTgg                                                                                                                              <  1:216972/148‑1 (MQ=255)                        gggACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGt                                                                                                                            >  7:229634/1‑148 (MQ=255)                             gacCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCtt                                                                                                                      >  7:4123/1‑149 (MQ=255)                                  tgtgCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAATTAAGCCCAATGGGGCTGCTAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCgc                                                                                                                 <  4:174428/149‑1 (MQ=255)                                   gtgCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCgcc                                                                                                                >  1:87708/1‑149 (MQ=255)                                    tgCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCgcct                                                                                                               <  1:287405/149‑1 (MQ=255)                                    tgCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCgcct                                                                                                               <  5:313110/149‑1 (MQ=255)                                      cgcTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCgcctgc                                                                                                             <  8:258919/149‑1 (MQ=255)                                         tttGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCattgg                                                                                                                                      <  6:90121/121‑5 (MQ=255)                                         tttGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCattgg                                                                                                                                      >  5:90121/1‑117 (MQ=255)                                             cGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTgg                                                                                                      >  8:225827/1‑149 (MQ=255)                                              gTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGt                                                                                                     >  6:272949/1‑149 (MQ=255)                                               ttGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGt                                                                                                     >  5:188118/1‑148 (MQ=255)                                                  aTCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTAc                                                                                                  <  1:116589/148‑1 (MQ=255)                                                    cGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTAc                                                                                                  >  1:74130/1‑146 (MQ=255)                                                    cGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTAc                                                                                                  <  2:74130/146‑1 (MQ=255)                                                     gagaACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCg                                                                                               <  3:301912/148‑1 (MQ=255)                                                     gagaACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGATT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCg                                                                                               >  1:187427/1‑148 (MQ=255)                                                      agaACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTa                                                                                             >  4:168820/1‑149 (MQ=255)                                                         aCCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACg                                                                                           <  2:87708/148‑1 (MQ=255)                                                         aCCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGc                                                                                          <  6:188118/149‑1 (MQ=255)                                                               aCCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCg                                                                                    <  3:168820/149‑1 (MQ=255)                                                                ccGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGt                                                                                   >  8:225984/1‑149 (MQ=255)                                                                 cGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGt                                                                                   >  3:200309/1‑148 (MQ=255)                                                                   tATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAg                                                                                <  5:272949/149‑1 (MQ=255)                                                                        gCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACt                                                                            >  7:362389/1‑148 (MQ=255)                                                                        gCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACt                                                                            >  5:287438/1‑148 (MQ=255)                                                                            ggCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCAttt                                                                       >  5:372920/1‑149 (MQ=255)                                                                             gCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTCGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCATACCGTCAGCCCTCAttt                                                                       >  4:271942/1‑148 (MQ=255)                                                                                aTCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCa                                                                   <  7:225984/149‑1 (MQ=255)                                                                                 tCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCAttacaaat                                                                  >  1:323396/1‑143 (MQ=255)                                                                                 tCCGCGGGAGCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAATTACGCCCAATGGGGCTGCGACGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCAttt                                                                       <  2:323396/144‑1 (MQ=255)                                                                                        gATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCaaa                                                              >  2:19960/1‑146 (MQ=255)                                                                                             ggCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTtctc                                                                    <  4:339456/135‑1 (MQ=255)                                                                                             ggCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTtctc                                                                    >  3:339456/1‑135 (MQ=255)                                                                                             ggCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATccc                                                      >  5:123971/1‑149 (MQ=255)                                                                                                aaaTGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCa                                                   <  1:19960/149‑1 (MQ=255)                                                                                                aaaTGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCa                                                   >  3:160134/1‑149 (MQ=255)                                                                                                   tGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGa                                                <  8:4123/149‑1 (MQ=255)                                                                                                     ttCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAgt                                              <  8:229634/149‑1 (MQ=255)                                                                                                        cGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTc                                           >  2:291986/1‑149 (MQ=255)                                                                                                                 ggTT‑TGGCTCGCATGTACCCCCAATGGGGCTGCTAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGc                                   <  4:200309/148‑1 (MQ=255)                                                                                                                           cGCAAGTAATCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGTGTGAGCTTTCCGCCTGCCTTCTGGTTAGCCGTACGGCTACCTTCAGCACTCATTTCTCACCAAAGAGATCCCCCTTGAGTGTCTTCACCGCCTACTTCTGc                         <  3:271942/149‑1 (MQ=255)                                                                                                                           cGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTc                                        <  1:281698/134‑1 (MQ=255)                                                                                                                            gCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGc                         <  7:225827/148‑1 (MQ=255)                                                                                                                                gTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCAc                    <  6:123971/149‑1 (MQ=255)                                                                                                                                 tAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACt                   >  5:140596/1‑149 (MQ=255)                                                                                                                                    gCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTtctc                                                                    <  3:346104/97‑1 (MQ=255)                                                                                                                                    gCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTtctc                                                                    >  4:346104/1‑97 (MQ=255)                                                                                                                                      ccAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATcccc                                                     >  7:107951/1‑110 (MQ=255)                                                                                                                                      ccAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATcccc                                                     <  8:107951/110‑1 (MQ=255)                                                                                                                                       cAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCg             >  6:106362/1‑149 (MQ=255)                                                                                                                                        aaTGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCg             >  4:348492/1‑148 (MQ=255)                                                                                                                                        aaTGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCg             >  1:195201/1‑148 (MQ=255)                                                                                                                                           ggggCTGCTTCGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCAcc         >  3:46777/1‑149 (MQ=255)                                                                                                                                           ggggCTGCTTCGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCAcc         >  5:84437/1‑149 (MQ=255)                                                                                                                                                  ctaaGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTcga  <  2:187427/147‑1 (MQ=255)                                                                                                                                                  cGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTcga  >  2:296947/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGGGATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTT‑TGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTCGA  >  NC_002947/956438‑956730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑