Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,321,070	C→T	8.4%	intergenic (‑306/‑112)	PP_3791 ← / → PP_5629	site‑specific recombinase/protease/amidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,321,070	0	C	T	8.4%	169.3 / 4.5	60	intergenic (‑306/‑112)	PP_3791/PP_5629	site‑specific recombinase/protease/amidase
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base T (2/3);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTCAATTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATACTGACACCG  >  NC_002947/4320923‑4321211                                                                                                                                                     |                                                                                                                                                 gAGTCAATTGAGTGTTATATAAGATGCGTCAGCTTCTCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCTACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATc                                                                                                                                                  <  3:249038/148‑1 (MQ=255)  aGTCAATTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATTGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGTAGCACAGTATTCTTAATT‑TCGATATTGAATcc                                                                                                                                                 >  4:59664/1‑148 (MQ=255)   gTCAATTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCaa                                                                                                                                               >  6:372707/1‑149 (MQ=255)    tCAATTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCCTAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCaa                                                                                                                                               >  2:131489/1‑148 (MQ=255)       aTTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATc                                                                                                                                           >  1:345640/1‑149 (MQ=255)           cGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGGATTCAATATcgaaa                                                                                                                                       <  8:75897/149‑4 (MQ=255)           cGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGa                                                                                                                                       >  4:318016/1‑149 (MQ=255)              gtTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCg                                                                                                                                    >  2:313258/1‑149 (MQ=255)                 aTCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGGATTCAATATcgaaatt                                                                                                                                     <  5:367676/145‑6 (MQ=255)                 aTCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGGATTCAATATcgaaatt                                                                                                                                     >  6:367676/1‑140 (MQ=255)                 aTCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCggg                                                                                                                                  <  1:313258/148‑1 (MQ=255)                   cTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGaa                                                                                                                                >  4:169743/1‑148 (MQ=255)                    tAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGa                                                                                                                                 >  4:228571/1‑146 (MQ=255)                    tAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGa                                                                                                                                 <  3:228571/146‑1 (MQ=255)                    tAAGATGCGTCAGCTTCGCTTTCGCCATCATAATCGTAGTCAATCAAGCTTGCCCACACTTCACAAATTTCAAAATACCCTAACGCCCTGGTTCCAGCACCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAAca                                                                                                                              >  2:106352/1‑149 (MQ=255)                      aGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAAc                                                                                                                               <  7:319788/146‑1 (MQ=255)                       gATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAAcatca                                                                                                                           <  3:318016/149‑1 (MQ=255)                         tGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACatcatc                                                                                                                         <  3:133808/149‑1 (MQ=255)                                   tCGCTTTCGCCCTCCTAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTa                                                                                                               <  1:131489/149‑1 (MQ=255)                                           gCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑GCGATATTGGATTCAATATcgaaa                                                                                                                                       <  8:283059/117‑4 (MQ=255)                                           gCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑GCGATATTGGATTCAATATcgaaa                                                                                                                                       >  7:283059/1‑114 (MQ=255)                                             ccTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGc                                                                                                     <  5:372707/149‑1 (MQ=255)                                             ccTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGc                                                                                                     >  5:93008/1‑149 (MQ=255)                                                 aTAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑a                                                                                                  <  4:206764/147‑1 (MQ=255)                                                  tAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCt                                                                                               <  1:342107/149‑1 (MQ=255)                                                       gTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑Ac                                                                                         <  3:266311/148‑1 (MQ=255)                                                          gTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCtt                                                                                      >  6:173587/1‑148 (MQ=255)                                                             aatcaaGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGAt                                                                                   <  6:54580/148‑1 (MQ=255)                                                             aatcaaGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATg                                                                                  >  3:78066/1‑149 (MQ=255)                                                              atcaaGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGt                                                                                 >  2:354138/1‑149 (MQ=255)                                                                caaGCTTGCCGACACTTGACAAATTTC‑AATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAAcatca                                                                                                                           <  6:158160/107‑1 (MQ=255)                                                                caaGCTTGCCGACACTTGACAAATTTC‑AATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAAcatca                                                                                                                           >  5:158160/1‑107 (MQ=255)                                                                   gCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGtt                                                                            <  2:345640/149‑1 (MQ=255)                                                                            cacTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGTAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAgccgc                                                                   <  3:59664/149‑1 (MQ=255)                                                                               ttGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGccgccg                                                                 >  3:17290/1‑148 (MQ=255)                                                                               ttGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGccgccg                                                                 >  1:147828/1‑148 (MQ=255)                                                                                    aaaTTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTActc                                                            >  4:363822/1‑146 (MQ=255)                                                                                     aaTTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTa                                                          >  8:238012/1‑149 (MQ=255)                                                                                       tttCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATc                                                        >  1:82592/1‑149 (MQ=255)                                                                                              ttACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTGTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCAtt                                                  >  1:256587/1‑149 (MQ=255)                                                                                                aCCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCa                                                <  3:169743/148‑1 (MQ=255)                                                                                                 cccTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATg                                              <  6:135518/149‑1 (MQ=255)                                                                                                   cTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATgga                                            <  4:17290/149‑1 (MQ=255)                                                                                                    tAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATgg                                             >  8:289117/1‑147 (MQ=255)                                                                                                    tAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATgg                                             <  7:289117/147‑1 (MQ=255)                                                                                                        gCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAAt                                       <  5:173587/149‑1 (MQ=255)                                                                                                          ccTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATgg                                             >  4:62136/1‑141 (MQ=255)                                                                                                          ccTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATgg                                             <  3:62136/141‑1 (MQ=255)                                                                                                           cTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTg                                    <  4:78066/149‑1 (MQ=255)                                                                                                                 ccAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGAt                              <  2:82592/149‑1 (MQ=255)                                                                                                                   agcatcACATTTTTATTAATT‑TCGATATTGAATACAATATAGCGACCGGGAACATCATCCTCCCTTCTATTAGCACCGCT‑AGCTTATCTCCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGAt                              <  1:106352/149‑1 (MQ=255)                                                                                                                   agcagcACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATgg                            >  7:172294/1‑149 (MQ=255)                                                                                                                       gcaCAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCAcg                        >  7:162960/1‑149 (MQ=255)                                                                                                                               ttCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCAc                         <  7:302045/140‑1 (MQ=255)                                                                                                                               ttCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCAc                         >  8:302045/1‑140 (MQ=255)                                                                                                                               ttCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGtt                <  6:93008/149‑1 (MQ=255)                                                                                                                                  ttAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGc                                  >  7:319886/1‑128 (MQ=255)                                                                                                                                  ttAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGc                                  <  8:319886/128‑1 (MQ=255)                                                                                                                                    aaTT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTg              <  7:238012/146‑1 (MQ=255)                                                                                                                                      tt‑tCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATACt         >  7:94913/1‑149 (MQ=255)                                                                                                                                      tt‑tCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATACt         <  1:354138/149‑1 (MQ=255)                                                                                                                                               ttGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAAGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATACTGACACCg  >  1:280943/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GAGTCAATTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATACTGACACCG  >  NC_002947/4320923‑4321211

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