Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,992	A→T	9.7%	F248I (TTC→ATC)	PP_4563 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,992	0	A	T	9.7%	173.3 / 8.2	62	F248I (TTC→ATC)	PP_4563	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base T (3/3);  total (30/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.46e-01

CGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTC  >  NC_002947/5182844‑5183136                                                                                                                                                     |                                                                                                                                                  cGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa                                                                                                                                                   <  3:39956/149‑1 (MQ=255) cGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa                                                                                                                                                   <  1:372/149‑1 (MQ=255) cGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGGAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa                                                                                                                                                   <  8:137398/149‑1 (MQ=255)        gTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTg                                                                                                                                            >  1:132602/1‑149 (MQ=255)         tGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGc                                                                                                                                           <  7:56992/149‑1 (MQ=255)          ggCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAAGAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCt                                                                                                                                          >  1:282435/1‑149 (MQ=255)            cGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                        <  3:252362/149‑1 (MQ=255)            cGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                        <  5:345910/149‑1 (MQ=255)                 agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCaa                                                                                                                                    <  7:57472/148‑1 (MQ=255)                 agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAAc                                                                                                                                   <  2:272370/149‑1 (MQ=255)                  gcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACg                                                                                                                                  >  7:21781/1‑149 (MQ=255)                        cTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATTGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGt                                                                                                                            <  5:2517/149‑1 (MQ=255)                           ccGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAgg                                                                                                                              >  7:169151/1‑144 (MQ=255)                           ccGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAgg                                                                                                                              <  8:169151/144‑1 (MQ=255)                           ccGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCAc                                                                                                                         >  1:12718/1‑149 (MQ=255)                             gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                                                                                       <  6:332457/149‑1 (MQ=255)                               aaCTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCtt                                                                                                                     <  2:282435/149‑1 (MQ=255)                                    cGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                  <  1:4652/147‑1 (MQ=255)                                         aGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAcc                                                                                                           >  6:17491/1‑149 (MQ=255)                                         aGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAcc                                                                                                           <  8:21781/149‑1 (MQ=255)                                              tAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCaa                                                                                                      <  7:135729/149‑1 (MQ=255)                                                   ccGTCAGGCTTTCTCGATGCGATCGTCATGGATTACGATCAGGCCCTGCTTGAACAACCAGCCAATGGCCTTCTTGAAGTTGCACTTGCTGACGTTGAACAGCTTGCTGATCAACGAAGGGTCACTCTTGTCGCACAACCCCAacacac                                                                                                 >  3:157906/1‑147 (MQ=255)                                                     gTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCg                                                                                               >  8:35066/1‑149 (MQ=255)                                                        aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCt                                                                                            >  7:313002/1‑149 (MQ=255)                                                           cTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCat                                                                                                      <  8:1188/136‑2 (MQ=255)                                                           cTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCat                                                                                                      >  7:1188/1‑135 (MQ=255)                                                                  gATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCtt                                                                                     <  5:17491/146‑1 (MQ=255)                                                                  gATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCtt                                                                                     <  8:313002/146‑1 (MQ=255)                                                                   aTGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑Ag                                                                                <  5:45611/149‑1 (MQ=255)                                                                   aTGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑Ag                                                                                <  2:132602/149‑1 (MQ=255)                                                                     gCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCACCACCCCGCCTTCGGCTTCC‑AGgc                                                                              >  8:41317/1‑149 (MQ=255)                                                                       gATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGgcgc                                                                            <  3:313289/149‑1 (MQ=255)                                                                           gTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGcc                                                                         <  5:81872/148‑1 (MQ=255)                                                                            tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGc                                                                                                                                           >  2:70071/1‑82 (MQ=255)                                                                            tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGc                                                                                                                                           <  1:70071/82‑1 (MQ=255)                                                                            tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCAt                                                                       >  7:22412/1‑149 (MQ=255)                                                                                 gATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATc                                                                   <  8:22412/148‑1 (MQ=255)                                                                                 gATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCt                                                                  >  4:35744/1‑149 (MQ=255)                                                                                       gATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                   <  7:227845/79‑1 (MQ=38)                                                                                       gATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                   >  8:227845/1‑79 (MQ=38)                                                                                               ccTGCTTGAACAACCCGCCAATGGCCTTCTTGATGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGAAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGgtg                                                    >  4:53250/1‑149 (MQ=255)                                                                                               ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                       >  4:290153/1‑146 (MQ=255)                                                                                               ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                       <  3:290153/146‑1 (MQ=255)                                                                                               ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGgtg                                                    >  3:163669/1‑149 (MQ=255)                                                                                                      gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGTTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGGCTTCGGCTTCCAAGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                 <  1:41846/146‑1 (MQ=255)                                                                                                      gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                 >  2:41846/1‑145 (MQ=255)                                                                                                      gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                 <  3:164412/145‑1 (MQ=255)                                                                                                      gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                 >  4:164412/1‑145 (MQ=255)                                                                                                         caacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAggg                                          <  7:152127/149‑1 (MQ=255)                                                                                                             ccGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGc                                               >  4:37135/1‑140 (MQ=255)                                                                                                             ccGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGc                                               <  3:37135/140‑1 (MQ=255)                                                                                                                       tcttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGc                            >  3:115663/2‑149 (MQ=255)                                                                                                                         ttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTg                          >  1:323281/1‑149 (MQ=255)                                                                                                                         ttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACAACCCGCCTTCGGATTAC‑AGGCTCGCCATGATATGCTCTTGCAGGGTGTCCGCATCGGCTGCGCCTACAGGCt                           >  1:227391/1‑148 (MQ=255)                                                                                                                             tGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAc                                <  2:54764/139‑1 (MQ=255)                                                                                                                             tGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAc                                >  1:54764/1‑139 (MQ=255)                                                                                                                                 gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                       >  1:264821/1‑112 (MQ=255)                                                                                                                                 gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                       <  2:264821/112‑1 (MQ=255)                                                                                                                                  ttGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAg                 >  5:278332/1‑149 (MQ=255)                                                                                                                                           cctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCtt        >  4:239915/5‑149 (MQ=255)                                                                                                                                            tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg       >  7:122502/1‑149 (MQ=255)                                                                                                                                            ctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAgg                      <  4:210065/131‑1 (MQ=255)                                                                                                                                            ctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAgg                      >  3:210065/4‑134 (MQ=255)                                                                                                                                              aCGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGcc     <  3:35744/149‑1 (MQ=255)                                                                                                                                                 ttGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTc  <  2:12718/149‑1 (MQ=255)                                                                                                                                                    aaCAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAgg                                           <  7:6805/105‑1 (MQ=255)                                                                                                                                                    aaCAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAgg                                           >  8:6805/1‑105 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTC  >  NC_002947/5182844‑5183136

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