Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,237,543	G→A	8.2%	M1L (TTG→TTA) †‡	PP_4614 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,237,543	0	G	A	8.2%	174.5 / 5.0	61	M1L (TTG→TTA) †‡	PP_4614	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base A (3/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CCAGTTCGCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGG  >  NC_002947/5237396‑5237688                                                                                                                                                          |                                                                                                                                                  cgaGTTCGCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGa                                                                                                                                                  <  3:22965/147‑1 (MQ=255)       cGCTGGACCAGGCAGTCGTCGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGcc                                                                                                                                               <  6:172032/146‑1 (MQ=255)        gCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGc                                                                                                                                           >  1:102785/1‑149 (MQ=255)        gCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGc                                                                                                                                           >  2:314180/1‑149 (MQ=255)         cTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCg                                                                                                                                          <  4:212149/149‑1 (MQ=255)            caCCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCTCAAGGCCGGGCGttt                                                                                                                                       <  4:16084/148‑1 (MQ=255)              ccAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTc                                                                                                                                     <  7:310947/149‑1 (MQ=255)               cAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCCATTGATCGGCACAAACGAAAACGCCCGGCTCAAGGCCGGGCGtt                                                                                                                                        >  3:129201/1‑149 (MQ=25)               cAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCTCAAGGCCGGGCGttttga                                                                                                                                    >  7:182129/1‑147 (MQ=255)               cAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCTCAAGGCCGGGCGtttt                                                                                                                                      <  8:182129/147‑1 (MQ=255)                  gCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGttt                                                                                                                                 <  6:331802/149‑1 (MQ=255)                         aCGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGa                                                                                                                          >  4:223766/1‑149 (MQ=255)                             tGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCaa                                                                                                                      <  8:6440/149‑1 (MQ=255)                                    accaAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑GAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGgcg                                                                                                                 >  6:143020/1‑147 (MQ=255)                                    accaAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑GAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGgcg                                                                                                                 <  5:143020/147‑1 (MQ=255)                                       aaGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGaa                                                                                                            <  3:223766/149‑1 (MQ=255)                                        aGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAAc                                                                                                           <  8:10980/149‑1 (MQ=255)                                          cTCAATCACGCGCTGGAGCACTTCTACGTCAATCTGCGCAGCATTAGCGGCGCCTTGGAACTTTGTGTCCCGGCCATGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAAc                                                                                                           <  1:275187/149‑1 (MQ=255)                                             aaTCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTc                                                                                                                                     >  6:271055/1‑118 (MQ=255)                                             aaTCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGgcg                                                                                                                 <  6:22242/138‑1 (MQ=255)                                             aaTCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGgcg                                                                                                                 >  5:22242/1‑138 (MQ=255)                                             aaTCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACa                                                                                                      >  5:89200/1‑149 (MQ=255)                                                 aCGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGc                                                                                                   >  7:64271/1‑148 (MQ=255)                                                         gAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCg                                                                                          >  2:104878/1‑149 (MQ=255)                                                          aGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCgg                                                                                         <  7:258212/149‑1 (MQ=255)                                                              cTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCaa                                                                                     >  6:139119/1‑149 (MQ=255)                                                               ttCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAAc                                                                                    >  4:207506/1‑149 (MQ=255)                                                                tCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAAc                                                                                    <  6:55318/148‑1 (MQ=255)                                                                            tGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCAtcact                                                                        >  6:5348/1‑148 (MQ=255)                                                                                   caAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAg                                                                <  6:89200/149‑1 (MQ=255)                                                                                   caAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAg                                                                <  1:314180/149‑1 (MQ=255)                                                                                      gcgcGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAg                                                             <  8:64271/149‑1 (MQ=255)                                                                                      gcgcGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAg                                                             >  3:111958/1‑149 (MQ=255)                                                                                               cTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCAt                                                    >  2:246629/1‑149 (MQ=255)                                                                                                 tGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATg                                                   >  3:218130/1‑148 (MQ=255)                                                                                                   gAACTGAGTGCCCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGaa                                                 <  7:149395/148‑1 (MQ=255)                                                                                                    aaCTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGaac                                                >  3:245300/1‑148 (MQ=255)                                                                                                         aGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCa                                          >  7:297694/1‑149 (MQ=255)                                                                                                           tGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCAtt                                        <  3:207506/149‑1 (MQ=255)                                                                                                               ccGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTcgaa                                    <  4:245300/149‑1 (MQ=255)                                                                                                                        cc‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCa                             <  1:104878/149‑1 (MQ=255)                                                                                                                            tGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGt                           >  5:80668/1‑149 (MQ=255)                                                                                                                            tGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGt                           >  7:227369/1‑149 (MQ=255)                                                                                                                               tCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGtgct                        <  8:297694/149‑1 (MQ=255)                                                                                                                                    a‑aTGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTcgaac                                   <  8:32255/133‑1 (MQ=255)                                                                                                                                    a‑aTGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTcgaac                                   >  7:32255/1‑133 (MQ=255)                                                                                                                                      aaacgaaaaCGCCCGGCTCAAGGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCAtt                   >  1:184577/6‑149 (MQ=255)                                                                                                                                      aTGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTc                  <  1:246629/149‑1 (MQ=255)                                                                                                                                           aaaaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGAAAGACCATGAACAACACATTCGAACGAATCAGTGCTGCATTCGcac              >  8:216715/1‑147 (MQ=255)                                                                                                                                            aaaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGcccc             <  4:218130/149‑1 (MQ=255)                                                                                                                                             aaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCAtt                                        <  5:123494/121‑1 (MQ=255)                                                                                                                                             aaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCAtt                                        >  6:123494/1‑121 (MQ=255)                                                                                                                                             aaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGcc               >  5:345667/1‑146 (MQ=255)                                                                                                                                             aaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGcc               <  6:345667/146‑1 (MQ=255)                                                                                                                                             aaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATCGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCAtt                                        >  2:92784/1‑121 (MQ=255)                                                                                                                                             aaCGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATCGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCAtt                                        <  1:92784/121‑1 (MQ=255)                                                                                                                                                 cccGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGaa         >  7:126874/1‑148 (MQ=255)                                                                                                                                                   cGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATAc      >  5:296765/1‑149 (MQ=255)                                                                                                                                                      ccTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCgg   <  5:5348/149‑1 (MQ=255)                                                                                                                                                       cTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCggg  <  7:216715/149‑1 (MQ=255)                                                                                                                                                       cTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCggg  <  8:227369/149‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                  CCAGTTCGCTGGACCAGGCAGTCGACGGTGCCGCCACCAAGCTCAATCACGCGCTGGAGCACTTCTACGGCAAGCTGCGCAGCAAGCGCGGCGCCTTGGAACTGAGTGACCCGG‑‑ACGCC‑‑TGATCGGCA‑ATGC‑AAAACGCCCGGCCTTGAGCCGGGCGTTTTCGTTTGTGCCGATCAATGGCGCTGAACCAACATGCGAACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGG  >  NC_002947/5237396‑5237688

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