Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,379,146	G→C	9.2%	A116P (GCT→CCT)  P65A (CCT→GCT)	PP_5639 → PP_3855 ←	hypothetical protein hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,379,146	0	G	C	9.2%	198.6 / 8.8	65	A116P (GCT→CCT)  P65A (CCT→GCT)	PP_5639 PP_3855	hypothetical protein hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base C (3/3);  total (33/32)
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.72e-01

TACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTT  >  NC_002947/4379008‑4379293                                                                                                                                           |                                                                                                                                                    tACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCg                                                                                                                                             >  1:70603/1‑147 (MQ=255)   cGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCAGAGCCTGCGCCGGAGc                                                                                                                                         <  3:208493/149‑1 (MQ=255)         tCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCg                                                                                                                                    >  7:148394/1‑148 (MQ=255)           tGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCaa                                                                                                                                 >  4:239616/1‑149 (MQ=255)             aTCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGc                                                                                                                                         <  3:93157/139‑1 (MQ=255)             aTCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAAt                                                                                                                                >  8:34655/1‑148 (MQ=255)             aTCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTACCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCAGAGCCTGCGCCGGAGc                                                                                                                                         <  2:323666/139‑1 (MQ=255)             aTCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTACCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCAGAGCCTGCGCCGGAGc                                                                                                                                         >  1:323666/1‑139 (MQ=255)              tCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGa                                                                                                                                           <  6:75275/136‑1 (MQ=255)              tCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGa                                                                                                                                           >  5:75275/1‑136 (MQ=255)                gACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGAt                                                                                                                             <  7:26137/148‑1 (MQ=255)                    tGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGaca                                                                                                                        <  3:260551/149‑1 (MQ=255)                          cGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTa                                                                                                                   <  8:148394/148‑1 (MQ=255)                          cGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTa                                                                                                                   <  5:253014/148‑1 (MQ=255)                          cGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTa                                                                                                                   <  3:239616/148‑1 (MQ=255)                           gTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAgc                                                                                                                 >  1:39149/1‑149 (MQ=255)                                         acacAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCtt                                                                                                     <  2:70603/147‑1 (MQ=255)                                         acacAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCtt                                                                                                     <  6:297379/147‑1 (MQ=255)                                         acacAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTc                                                                                                    <  2:39149/148‑1 (MQ=255)                                             aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCGGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTa                                                                                                                   <  4:89924/129‑1 (MQ=255)                                             aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTa                                                                                                                   >  3:89924/1‑129 (MQ=255)                                             aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCa                                                                                                   >  7:271321/1‑145 (MQ=255)                                             aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCa                                                                                                   <  8:271321/145‑1 (MQ=255)                                             aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTAc                                                                                               >  2:122208/1‑149 (MQ=255)                                                          tACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCtt                                                                                                     >  8:24287/1‑130 (MQ=255)                                                          tACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCtt                                                                                                     <  7:24287/130‑1 (MQ=255)                                                             acCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATAcc                                                                               >  2:14537/1‑149 (MQ=255)                                                               cAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTg                                                                             >  5:338282/1‑149 (MQ=255)                                                                       aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCg                                                                     >  3:246496/1‑149 (MQ=255)                                                                        gCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                    <  1:14537/149‑1 (MQ=255)                                                                        gCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                    <  6:338282/149‑1 (MQ=255)                                                                        gCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGAGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                    >  4:2784/1‑149 (MQ=255)                                                                               aaGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgc                                                                            <  6:249763/134‑1 (MQ=255)                                                                               aaGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgc                                                                            >  5:249763/1‑134 (MQ=255)                                                                                 gCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACt                                                           >  4:335282/1‑149 (MQ=255)                                                                                    aGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCt                                                                                                                     <  2:240196/88‑1 (MQ=255)                                                                                    aGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCt                                                                                                                     >  1:240196/1‑88 (MQ=255)                                                                                          tatGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATcgc                                                  <  1:122208/149‑1 (MQ=255)                                                                                               tGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCg                                                                                        <  4:149232/106‑1 (MQ=255)                                                                                               tGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCg                                                                                        >  3:149232/1‑106 (MQ=255)                                                                                                   aCCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGtt                                                                                     >  4:116857/1‑105 (MQ=255)                                                                                                   aCCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGtt                                                                                     <  3:116857/105‑1 (MQ=255)                                                                                                    ccAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACgg                                        >  2:172421/1‑149 (MQ=255)                                                                                                      aaGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATACTTCAGTACCTGGTCGGTTAATACCTACGAAGTCGGGAACTGACTCAGGATCGCGCGGGCACgcg                                       >  8:312747/1‑146 (MQ=255)                                                                                                         cTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATAc                                                                                <  1:324583/104‑1 (MQ=255)                                                                                                         cTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATAc                                                                                >  2:324583/1‑104 (MQ=255)                                                                                                         cTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATc                                   <  3:2784/149‑1 (MQ=255)                                                                                                            cAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCa                                  <  7:34655/147‑1 (MQ=255)                                                                                                             aGCCCGCAGATGCATCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTAGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATAAGcc                               >  2:140955/1‑149 (MQ=255)                                                                                                              gCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGccc                              >  2:266485/1‑149 (MQ=255)                                                                                                               cccGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGcc                               >  3:128602/1‑147 (MQ=255)                                                                                                               cccGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGcc                               <  4:128602/147‑1 (MQ=255)                                                                                                                     gATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCg                                                                     <  6:312923/103‑1 (MQ=255)                                                                                                                     gATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCg                                                                     >  5:312923/1‑103 (MQ=255)                                                                                                                     gATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCtg                       <  3:196200/149‑1 (MQ=255)                                                                                                                      aTGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGt                                                                                      <  3:236539/85‑1 (MQ=255)                                                                                                                      aTGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGt                                                                                      >  4:236539/1‑85 (MQ=255)                                                                                                                         ccTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTgga                   <  1:172421/149‑1 (MQ=255)                                                                                                                          cgaaGGCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggag                  >  6:51179/4‑149 (MQ=255)                                                                                                                           gaaGGCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCtgt                      >  3:198333/3‑144 (MQ=255)                                                                                                                           gaaGGCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCtgt                      <  4:198333/142‑1 (MQ=255)                                                                                                                                 tCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACg                                         <  7:180030/119‑1 (MQ=255)                                                                                                                                 tCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACg                                         >  8:180030/1‑119 (MQ=255)                                                                                                                                   cGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGc         >  8:223269/1‑149 (MQ=255)                                                                                                                                       gcAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCtgtt     >  1:221917/1‑149 (MQ=255)                                                                                                                                          ggCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTgttgtt  <  7:223269/149‑1 (MQ=255)                                                                                                                                           gCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCt        >  7:3649/1‑142 (MQ=255)                                                                                                                                           gCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCt        <  8:3649/142‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                    TACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTT  >  NC_002947/4379008‑4379293

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