Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,201,440	G→C	6.6%	F25L (TTC→TTG)	PP_4582 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,201,440	0	G	C	6.6%	198.7 / 4.3	61	F25L (TTC→TTG)	PP_4582	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base C (2/2);  total (28/33)
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.98e-01

GAAAGCGACATCCAGCGCAGTGATGGCGATTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTG  >  NC_002947/5201296‑5201587                                                                                                                                                 |                                                                                                                                                    gAAAGCGACATCCAGCGCAGTGATGGCGATTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGaa                                                                                                                                                   <  2:26546/147‑1 (MQ=255) gAAAGCGACATCCAGCGCAGTGATGGCGATTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAg                                                                                                                                                  <  6:143955/148‑1 (MQ=255) gAAAGCGACATCCAGCGCAGTGATGGCGATTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAg                                                                                                                                                  <  8:235127/148‑1 (MQ=255)       gACATCCAGCGCAGTGATGGCGATTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCccac                                                                                                                                             >  5:255535/1‑147 (MQ=255)                       aTGGCGATTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAgg                                                                                                                            <  6:310463/148‑1 (MQ=255)                          gCGATTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAggcggc                                                                                                                        >  2:89776/1‑149 (MQ=255)                          gCGATTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAggcggc                                                                                                                        <  6:255535/149‑1 (MQ=255)                          gCGATTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAggcgg                                                                                                                         <  8:167427/148‑1 (MQ=255)                             aTTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTAc                                                                                                                     <  3:233788/149‑1 (MQ=255)                                 cAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCg                                                                                                                 >  8:283381/1‑149 (MQ=255)                                  aGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGCTGAGCTTGAAGCCCACCAGCAGGCACCCAAAGGCGGCTACTGcc                                                                                                                 >  2:274796/1‑147 (MQ=255)                                   ggggCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTa                                                                                                               <  4:161473/149‑1 (MQ=255)                                       cGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAg                                                                                                           <  7:126947/149‑1 (MQ=255)                                          aGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGccc                                                                                                        <  1:83159/149‑1 (MQ=255)                                          aGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGccc                                                                                                        >  7:244389/1‑149 (MQ=255)                                             tGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAg                                                                                                     >  5:218434/1‑149 (MQ=255)                                              gTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGACTACTGCGTACAAGCCCCAgg                                                                                                    >  4:183451/1‑149 (MQ=255)                                               tcattGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCAACGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCAACAGCAGGCACCCTAAGGCGGCTACTGCGTAAAAccc                                                                                                         >  4:283619/4‑143 (MQ=255)                                               tCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGcc                                                                                                         <  3:283619/143‑1 (MQ=255)                                                ccTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGc                                                                                                   >  4:184225/1‑148 (MQ=255)                                                  ttGACCGGGTCGACGATGGTCGGGGTCAGCAGAAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGACAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTTCGTACAAGCCCCAGGCg                                                                                                  >  1:273246/1‑147 (MQ=255)                                                    gcccGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAg                                                                                               <  6:218434/146‑1 (MQ=255)                                                    gACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGaa                                                                                                <  1:56491/147‑1 (MQ=255)                                                    gACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAg                                                                                               <  5:38367/148‑1 (MQ=255)                                                    gACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGACTACTGCGTACAAGCCCCAGGCGaa                                                                                                <  3:183451/147‑1 (MQ=255)                                                     aCCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGgt                                                                                             >  8:243766/1‑149 (MQ=255)                                                         ggTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAg                                                                                         <  7:15331/149‑1 (MQ=255)                                                               cgaTGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTtc                                                                                   <  3:184225/149‑1 (MQ=255)                                                                     tCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCgg                                                                               >  6:130775/1‑147 (MQ=255)                                                                     tCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCgg                                                                               <  5:130775/147‑1 (MQ=255)                                                                     tCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGt                                                                              >  5:45151/1‑148 (MQ=255)                                                                     tCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGt                                                                              >  1:153266/1‑148 (MQ=255)                                                                        gggTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATg                                                                          <  1:274796/149‑1 (MQ=255)                                                                         ggTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATgg                                                                         >  3:290576/1‑149 (MQ=255)                                                                             agcagcagcACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGtt                                                                     >  6:184905/1‑149 (MQ=255)                                                                              gcagcagcaCCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTc                                                                    <  5:52775/149‑1 (MQ=255)                                                                                agcagcaCCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGt                                                                  <  6:45151/149‑1 (MQ=255)                                                                                   agcaCCGCCAGCACCACGCGTAGCGGTTGGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGc                                                                >  1:195006/1‑148 (MQ=255)                                                                                               accacGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACTAGCCCCAGGCGAAGGTGTAGt                                                                                        >  2:307616/1‑112 (MQ=255)                                                                                               accacGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACTAGCCCCAGGCGAAGGTGTAGt                                                                                        <  1:307616/112‑1 (MQ=255)                                                                                                ccacGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCaaa                                                    >  6:237926/1‑147 (MQ=255)                                                                                                ccacGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAaga                                                  <  7:243766/149‑1 (MQ=255)                                                                                                  aatcgTAGCGGTTCGGGGATCTAGCGCCACATTCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCAAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAgaga                                                <  6:240948/146‑1 (MQ=255)                                                                                                   ctcgTAGCGGTTCGCGCAGCCTGCGTCACATCCAGCCGGTGAGCTTGATGCGCAACAGCAGGCTCCCTAAGGCGGCTACTGCGTACAATCCCCAGGCGAATGTGTAGTCGTTCTGGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGAt                                               <  2:273246/147‑1 (MQ=255)                                                                                                     cgTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGc                                             <  2:195006/149‑1 (MQ=255)                                                                                                         gCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTAt                                         <  5:238787/149‑1 (MQ=255)                                                                                                                  tcaggcagcGCCACATCAAGCCGGTGAGCTTGAAGCCCACCATCAGGCACCCTAAGGCGGCTACTGCGTACAAGCACCTGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCaaa                                                    <  4:58314/128‑1 (MQ=255)                                                                                                                  gcAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCaaa                                                    >  3:58314/1‑129 (MQ=255)                                                                                                                  gcAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGGCGTACTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAacac                                <  1:89776/149‑1 (MQ=255)                                                                                                                                 tCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGa                 >  7:166220/1‑149 (MQ=255)                                                                                                                                 tCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGa                 >  7:171127/1‑149 (MQ=255)                                                                                                                                 tCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGa                 <  8:244389/149‑1 (MQ=255)                                                                                                                                  ccAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGa                 >  5:76116/1‑148 (MQ=255)                                                                                                                                       tggTGGGCTTCAAGCTCACCAGCAGGCACCGTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATaaa                                   >  1:27292/2‑125 (MQ=255)                                                                                                                                       tggTGGGCTTCAAGCTCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATaaa                                   <  2:27292/124‑1 (MQ=255)                                                                                                                                         gTGGGCTTCAAGCTCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATg                                        >  4:278477/1‑118 (MQ=255)                                                                                                                                         gTGGGCTTCAAGCTCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATg                                        <  3:278477/118‑1 (MQ=255)                                                                                                                                         gTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGttt         <  2:153266/149‑1 (MQ=255)                                                                                                                                               ttGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCtgt   >  6:337343/1‑149 (MQ=255)                                                                                                                                               ttGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGATCTCGGTCATGGTGTTCGTGAAAGCCAGGCAAAGAGATGCCTATGATAAACACATATCCGGGCGCAGACAGTGTTTGCCtgt   >  2:168298/1‑149 (MQ=255)                                                                                                                                                tGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCtgtg  >  8:67132/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GAAAGCGACATCCAGCGCAGTGATGGCGATTGCAGGGGCGAAGCTGTCCTTGACCGGGTCGACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTG  >  NC_002947/5201296‑5201587

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