Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,703,267	Δ1 bp	7.6%	coding (91/270 nt)	PP_5596 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,703,267	0	G	.	7.6%	246.9 / 13.5	79	coding (91/270 nt)	PP_5596	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (40/33);  new base . (3/3);  total (43/36)
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.04e-01

TGGCCACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGCCCAA  >  NC_002947/3703121‑3703412                                                                                                                                                   |                                                                                                                                                    tGGCCACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAg                                                                                                                                                   >  4:323560/1‑149 (MQ=255)  ggCCACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATACGCCTGGAACGCCTGGACAATGCCGAgg                                                                                                                                                  >  3:21090/1‑149 (MQ=255)   gCCACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAg                                                                                                                                                   <  3:323560/147‑1 (MQ=255)      aCCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑g                                                                                                                                             <  8:144489/149‑1 (MQ=255)           ccGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCAtt                                                                                                                                        <  3:274466/149‑1 (MQ=255)           ccGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCAtt                                                                                                                                        <  6:32898/149‑1 (MQ=255)             gATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTgg                                                                                                                                      >  5:7235/1‑149 (MQ=255)               tcaatcaaAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCC‑AGGCCTCGGCATTgtcc                                                                                                                                    <  3:105373/149‑4 (MQ=255)                 aatcaaAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGc                                                                                                                                    >  1:153358/1‑147 (MQ=255)                     aaaaTGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTg                                                                                                                                       >  1:360234/1‑140 (MQ=255)                     aaaaTGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTg                                                                                                                                       <  2:360234/140‑1 (MQ=255)                     aaaaTGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTtg                                                                                                                              <  6:118640/149‑1 (MQ=255)                     aaaaTGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGAAGTTTCATCATACGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GCCATTgggaggatcg                                                                                                                              >  5:17776/1‑142 (MQ=255)                         tGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTg                                                                                                                                       <  5:35793/136‑1 (MQ=255)                         tGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTg                                                                                                                                       >  6:35793/1‑136 (MQ=255)                         tGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGAt                                                                                                                          <  8:127890/149‑1 (MQ=255)                         tGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGAt                                                                                                                          <  8:133094/149‑1 (MQ=255)                             aTCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑g                                                                                                                      <  2:153358/148‑1 (MQ=255)                                     agaagaGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑g                                                                                                                                             >  3:38235/1‑118 (MQ=255)                                     agaagaGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑g                                                                                                                                             <  4:38235/118‑1 (MQ=255)                                         gaggagGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGcatca                                                                                                         >  7:39159/1‑149 (MQ=255)                                         gaggagGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGcatca                                                                                                         >  3:90845/1‑149 (MQ=255)                                           ggaggCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCAtt                                                                                                       >  3:175549/1‑149 (MQ=255)                                           ggaggCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCAtt                                                                                                       >  6:353346/1‑149 (MQ=255)                                            gaggCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTAGAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTtgtg                                                                                                                            >  1:261663/1‑128 (MQ=255)                                            gaggCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTtgtg                                                                                                                            <  2:261663/128‑1 (MQ=255)                                            gaggCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGAATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTc                                                                                                      >  8:71393/1‑149 (MQ=255)                                               gCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGACGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTccc                                                                                                    >  1:355336/1‑148 (MQ=255)                                                  cccTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTccc                                                                                                    >  3:232716/1‑145 (MQ=255)                                                  cccTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTccc                                                                                                    <  4:232716/145‑1 (MQ=255)                                                    cTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTc                                                                                              <  5:353346/149‑1 (MQ=255)                                                       cGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCg                                                                                             <  4:90845/147‑1 (MQ=255)                                                              tAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCgag                                                                                     >  2:100896/1‑148 (MQ=255)                                                                gCCCACAGCTGTACATAATCGATTACAAACTGCATGACAAACCGCGCAGTTTCATCATCCGCCTGGAACGCCAGGACAAAACCGAGGCCT‑GGCAATGGCAAAGTAGTGAAGC‑GGGAATAGCAAAAAATCCCAAGATCAGGCGCGAGaa                                                                                   >  2:332397/1‑148 (MQ=255)                                                                   cacaGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAAt                                                                                >  2:260773/1‑148 (MQ=255)                                                                    acaGCTGTACATATTCGATTACAAACTGCATGACCAACCGCGAAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCAT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCa                                                                              <  7:93427/149‑1 (MQ=255)                                                                       gCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCTCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGa                                                                           <  4:21090/149‑1 (MQ=255)                                                                            aCATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGt                                                                       >  5:323627/1‑148 (MQ=255)                                                                             cATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGt                                                                       >  5:7840/1‑147 (MQ=255)                                                                             cATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCa                                                                     <  7:71393/149‑1 (MQ=255)                                                                                aaTCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTc                                                                  >  5:128338/1‑149 (MQ=255)                                                                                    gATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCga                                                                                      <  1:265045/125‑1 (MQ=255)                                                                                    gATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCga                                                                                      >  2:265045/1‑125 (MQ=255)                                                                                     aTTACAAACTGCATGACCAACCGCGCAGTATCATCATCCGCCTGGAACGCCAGGACAAAGACGAGGCCT‑GGCATTGGGCAAGTTGTGATCC‑CGGCATTCGCATCAATCCCAACTTCGGTCGCGAGAAAAt                                                                                >  5:267780/1‑130 (MQ=255)                                                                                           aaCTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGcc                                                       <  6:323627/149‑1 (MQ=255)                                                                                             cTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCgg                                                     <  2:355336/149‑2 (MQ=255)                                                                                                      cAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCg                                                             <  8:64690/132‑1 (MQ=255)                                                                                                      cAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCg                                                             >  7:64690/1‑132 (MQ=255)                                                                                                      cAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACgg                                            >  1:89924/1‑149 (MQ=255)                                                                                                            cgcAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCg                                                                 <  3:192833/122‑1 (MQ=255)                                                                                                            cgcAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCg                                                                 >  4:192833/1‑122 (MQ=255)                                                                                                                   tcatcaTCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGt                                <  4:175549/148‑1 (MQ=255)                                                                                                                   tcatcaTCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGt                                >  7:170003/1‑148 (MQ=255)                                                                                                                   tcatcaTCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGt                                <  6:7235/148‑1 (MQ=255)                                                                                                                        atcCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTc                                                                                              >  5:329329/1‑81 (MQ=255)                                                                                                                        atcCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTc                                                                                              <  6:329329/81‑1 (MQ=255)                                                                                                                            gCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGCGGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCg                       >  1:91393/1‑149 (MQ=255)                                                                                                                            gCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCgg                      >  3:243654/1‑149 (MQ=255)                                                                                                                                ggAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGAttt                  >  8:201382/1‑149 (MQ=255)                                                                                                                                  aaCGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCg                >  7:291863/1‑149 (MQ=255)                                                                                                                                    cGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCTGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCggg              <  1:100896/149‑1 (MQ=255)                                                                                                                                     gCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCgggg             <  6:7840/149‑1 (MQ=255)                                                                                                                                       ttgcccAATGCC‑AGGCCTCGGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGc                                           <  5:51716/112‑1 (MQ=255)                                                                                                                                       ttgcccAATGCC‑AGGCCTCGGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGc                                           >  6:51716/6‑117 (MQ=255)                                                                                                                                        tGGACAATGCCGAGGCCT‑GGCATTGGGCAACTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAg           >  7:18965/1‑148 (MQ=255)                                                                                                                                         ggACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGc          >  5:141791/1‑148 (MQ=255)                                                                                                                                         ggACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGAGAGc          >  3:318492/1‑148 (MQ=255)                                                                                                                                          gACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCgg               >  4:317248/1‑142 (MQ=255)                                                                                                                                          gACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCgg               <  3:317248/142‑1 (MQ=255)                                                                                                                                          gACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGTGAGc          <  4:243654/147‑1 (MQ=255)                                                                                                                                          gACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGc          >  3:40584/1‑147 (MQ=255)                                                                                                                                          gACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCaa        >  6:39276/1‑149 (MQ=255)                                                                                                                                          cccAATGCC‑AGGCCTCGGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGc          >  3:162044/3‑147 (MQ=255)                                                                                                                                           ccAATGCC‑AGGCCTCGGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGa                                                     <  3:91759/102‑1 (MQ=255)                                                                                                                                           ccAATGCC‑AGGCCTCGGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGa                                                     >  4:91759/2‑103 (MQ=255)                                                                                                                                              aTGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGccc    <  7:201382/149‑1 (MQ=255)                                                                                                                                              aTGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGccc    <  4:40584/149‑1 (MQ=255)                                                                                                                                               tGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGccc    >  1:234470/1‑148 (MQ=255)                                                                                                                                                 ccGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGCCCaa  <  8:39159/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TGGCCACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGCCCAA  >  NC_002947/3703121‑3703412

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