Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,707,113	G→C	5.4%	L377L (CTG→CTC)	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,707,113	0	G	C	5.4%	220.9 / 3.9	74	L377L (CTG→CTC)	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/39);  new base C (2/2);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

GACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGG  >  NC_002947/5706967‑5707257                                                                                                                                                     |                                                                                                                                                  gACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCTGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTggg                                                                                                                                                 <  8:86837/149‑1 (MQ=255)    cTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGa                                                                                                                                              <  1:86938/149‑1 (MQ=255)      gCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATg                                                                                                                                            >  7:75364/1‑149 (MQ=255)        cTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTc                                                                                                                                          <  8:206479/149‑1 (MQ=255)         ttCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTc                                                                                                                                          >  6:306442/1‑148 (MQ=255)           cGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCCCAGCGATGtgca                                                                                                                                        <  5:215064/148‑4 (MQ=255)            gCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGaa                                                                                                                                       <  7:147627/148‑1 (MQ=255)              aGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGa                                                                                                                                     >  8:262140/1‑148 (MQ=255)                  acgaacTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGa                                                                                                                                        >  4:239311/1‑141 (MQ=255)                  acgaacTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGa                                                                                                                                        <  3:239311/141‑1 (MQ=255)                    gaacTGCGCCCGCTGTCAGCCCACCGCGATGGCGTGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCg                                                                                                                               >  1:47464/1‑149 (MQ=255)                    gaacTGCGCCCGCTGTCAGCCCACCGCGATGGCGTGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCg                                                                                                                               >  1:53295/1‑149 (MQ=255)                    gaacTGCGCCCGCAGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGt                                                                                                                              <  5:385010/149‑1 (MQ=255)                     aacTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGt                                                                                                                              >  3:81451/1‑148 (MQ=255)                       cTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑At                                                                                                                          <  5:306442/149‑1 (MQ=255)                       cTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑At                                                                                                                          >  2:232756/1‑149 (MQ=255)                         gcgcCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTa                                                                                                                        >  1:22917/1‑149 (MQ=255)                             ccGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTAcgcc                                                                                                                    >  3:177652/1‑149 (MQ=255)                                   tCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTAcgccgctcgc                                                                                                              <  4:81451/149‑1 (MQ=255)                                     accccACCGCGATGGA‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTCCACATCGCAGGGAGATCTCGAAGATATCCGTG‑ATTAcgccgctcgccg                                                                                                            >  2:304034/3‑149 (MQ=255)                                     aGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTAcgccgctcgccg                                                                                                            >  5:146068/1‑149 (MQ=255)                                     aGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTAcgccgctcgcc                                                                                                             <  3:195329/148‑1 (MQ=255)                                      gCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAg                                                                                                                                      >  2:338051/1‑123 (MQ=255)                                        ccACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCtt                                                                                                         <  3:87911/149‑1 (MQ=255)                                         cACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACCGGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCtt                                                                                                         >  1:64314/1‑149 (MQ=255)                                          aCCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGATCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGc                                                                                                       <  8:75364/149‑1 (MQ=255)                                            cgcgATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGt                                                                                                     <  2:22917/149‑1 (MQ=255)                                                   gc‑tgctCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGAcc                                                                                               <  2:221653/148‑1 (MQ=255)                                                    c‑tgctCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCt                                                                                             <  3:52731/149‑1 (MQ=255)                                                      tgctCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTg                                                                                            >  6:164504/1‑149 (MQ=255)                                                       gctCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTg                                                                                            >  3:131098/1‑148 (MQ=255)                                                        ctCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAGCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCa                                                                                          >  2:393100/1‑149 (MQ=255)                                                        ctCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTg                                                                                            <  7:112233/147‑1 (MQ=255)                                                         tCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTg                                                                                            >  4:345868/1‑146 (MQ=255)                                                           aaTGCCAGCTCGGAGGACC‑GCTACGCCTATCACATGGCACTTCTCGATCTCGCGTCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGTGCGATGTCGAAGATATCCGTG‑ATTACGACGCTCGCCGCTTGCGt                                                                                                     <  8:330879/135‑1 (MQ=255)                                                           aaTGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTga                                                                                                   <  5:68320/137‑1 (MQ=255)                                                           aaTGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTga                                                                                                   >  6:68320/1‑137 (MQ=255)                                                           aaTGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGAc                                                                                       <  7:262140/149‑1 (MQ=255)                                                           aaTGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGAc                                                                                       >  2:385200/1‑149 (MQ=255)                                                           aaTGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTAGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGt                                                                                                     >  7:330879/1‑135 (MQ=255)                                                            aTGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGAcc                                                                                      <  1:232756/149‑1 (MQ=255)                                                             tGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCa                                                                                     <  4:177652/149‑1 (MQ=255)                                                               ccAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAg                                                                                    <  2:64314/148‑1 (MQ=255)                                                                cAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAGCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCt                                                                                  <  1:393100/149‑1 (MQ=255)                                                                          gACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGGGATCTCGAAGATATCCGTG‑ATAACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCAATCCAGAc                                                                          >  3:104688/1‑147 (MQ=255)                                                                            cc‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAgg                                                                      <  1:385200/149‑1 (MQ=255)                                                                                cTACGCCAATCACACGGCACTGCTCGATGTGGCGGCGAACCTAGGTCGGTCGAATGGCTGCACATCGCTGGGCGATGTCGAATATTTGCGTGTTTTACGCCGCTCGCCGCTTGCGTGAGACCCTGCCGACCAGCTATCCAGACAAGGCa                                                                    <  1:305306/149‑1 (MQ=255)                                                                                   cGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGt                                                                  <  2:47464/147‑1 (MQ=255)                                                                                   cGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGt                                                                  <  2:53295/147‑1 (MQ=255)                                                                                   cGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTac                                                                <  3:373158/149‑1 (MQ=255)                                                                                   cGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTa                                                                 >  7:331816/1‑148 (MQ=255)                                                                                       tATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTaca                                                               <  3:224649/146‑1 (MQ=255)                                                                                       tATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTaca                                                               >  4:224649/1‑146 (MQ=255)                                                                                         ccacaCGCAACTGCTCGATCTGGCGGCTAACCAAGGTCGGTCGACGGGCTGCACATCGCGTGGCGATGTCGAAGATATCCGTT‑CTTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGAcc                                                          <  8:163850/148‑1 (MQ=255)                                                                                          cacaCGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCAGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGCAAAGGCAGTACACGAccc                                                         >  6:332334/1‑149 (MQ=255)                                                                                           acacGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGAccc                                                         >  7:41297/1‑148 (MQ=255)                                                                                              cGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCa                                                                    <  5:141004/134‑1 (MQ=255)                                                                                              cGGCAATGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCa                                                                    >  6:141004/1‑134 (MQ=255)                                                                                                   cTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAtgt                                                <  7:225351/149‑1 (MQ=255)                                                                                                      cTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAtgtg                                               <  4:131098/147‑1 (MQ=255)                                                                                                      cTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGc                                              <  3:345868/148‑1 (MQ=255)                                                                                                           tCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATcc                                        <  2:122638/149‑1 (MQ=255)                                                                                                               gcggcgAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGCCCCGAACGATGTGGACATCCGCAt                                    >  7:242188/1‑149 (MQ=255)                                                                                                                 ggcgAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCAttt                                  >  4:65940/1‑149 (MQ=255)                                                                                                                         aaGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCATACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGcc                          <  8:331816/149‑1 (MQ=255)                                                                                                                                   gAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACATGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGgcg                <  1:304034/149‑1 (MQ=255)                                                                                                                                     aGGGCGGCACATCGCTTGGCCATGTCGAATATTTCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGcc              <  8:242188/149‑1 (MQ=255)                                                                                                                                      gggCTGCACATCGCTGGGCGATGTCGATGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCg             <  3:167006/149‑1 (MQ=255)                                                                                                                                      gggCTGCACATCGCTGGGCGATGTCGATGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCg             <  3:172982/149‑1 (MQ=255)                                                                                                                                        cttcgacaTCGCCCAGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGCGCACATCCGCATTTCCCAAGCCATCATCGGCGCCgg            >  3:136869/6‑148 (MQ=255)                                                                                                                                         ttcgacaTCGCCCAGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAtgtg                                               >  7:212842/5‑112 (MQ=255)                                                                                                                                         ttcgacaTCGCCCAGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAtgtg                                               <  8:212842/108‑1 (MQ=255)                                                                                                                                              caTCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACaa     <  6:146068/149‑1 (MQ=255)                                                                                                                                                  gCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCGGCTTGCGTGATACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAggg  <  5:332334/148‑1 (MQ=255)                                                                                                                                                  gCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAggg  <  5:164504/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGC‑TGCTCAATGCCAGCGCTGAGGACC‑GCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTG‑ATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGG  >  NC_002947/5706967‑5707257

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