Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,321,067	A→G	7.9%	intergenic (‑303/‑115)	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,067	0	A	G	7.9%	176.9 / 5.0	63	intergenic (‑303/‑115)	PP_3791/PP_5629	site‑specific recombinase/protease/amidase
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base G (2/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

CTGGAGTCAATTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATACTGACACCG  >  NC_002947/4320920‑4321211                                                                                                                                                     |                                                                                                                                                    cTGGAGTCAATTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGaa                                                                                                                                                    >  8:319788/1‑149 (MQ=255) cTGGAGTCAATTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGaa                                                                                                                                                    >  4:133808/1‑149 (MQ=255)    gAGTCAATTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAAt                                                                                                                                                   <  8:327939/147‑1 (MQ=255)    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:93008/1‑149 (MQ=255)                                                ccTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGc                                                                                                     <  5:372707/149‑1 (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                                                                 >  1:147828/1‑148 (MQ=255)                                                                                  ttGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGccgccg                                                                 >  3:17290/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                                             <  3:62136/141‑1 (MQ=255)                                                                                                             ccTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATgg                                             >  4:62136/1‑141 (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)                                                                                                                                                     |                                                                                                                                                    CTGGAGTCAATTGCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATT‑TCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCT‑AGCTGA‑‑ACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATACTGACACCG  >  NC_002947/4320920‑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: ‑