Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,703,273:1	+C	8.2%	coding (97/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,273	1	.	C	8.2%	273.7 / 20.3	73	coding (97/270 nt)	PP_5596	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (37/30);  new base C (3/3);  total (40/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 = 7.09e-01

ACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGCCCAA  >  NC_002947/3703126‑3703412                                                                                                                                                     |                                                                                                                                             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                                                                                                                                       >  6:35793/1‑136 (MQ=255)                    tGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTg                                                                                                                                       <  5:35793/136‑1 (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                                                                                      >  2:265045/1‑125 (MQ=255)                                                                               gATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCga                                                                                      <  1:265045/125‑1 (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                                >  7:170003/1‑148 (MQ=255)                                                                                                              tcatcaTCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGt                                <  6:7235/148‑1 (MQ=255)                                                                                                              tcatcaTCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGt                                <  4:175549/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                                           >  6:51716/6‑117 (MQ=255)                                                                                                                                  ttgcccAATGCC‑AGGCCTCGGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGc                                           <  5:51716/112‑1 (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               <  3:317248/142‑1 (MQ=255)                                                                                                                                     gACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCgg               >  4:317248/1‑142 (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                                                     >  4:91759/2‑103 (MQ=255)                                                                                                                                      ccAATGCC‑AGGCCTCGGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGa                                                     <  3:91759/102‑1 (MQ=255)                                                                                                                                         aTGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGccc    <  4:40584/149‑1 (MQ=255)                                                                                                                                         aTGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGccc    <  7:201382/149‑1 (MQ=255)                                                                                                                                          tGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGccc    >  1:234470/1‑148 (MQ=255)                                                                                                                                            ccGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGCCCaa  <  8:39159/148‑1 (MQ=255)                                                                                                                                                    t‑GGCTTTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCggg              <  5:223952/128‑1 (MQ=255)                                                                                                                                                    t‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCggg              >  6:223952/1‑128 (MQ=255)                                                                                                                                                     |                                                                                                                                             ACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCT‑GGCATTGGGCAAGTTGTGATGC‑GGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGCCCAA  >  NC_002947/3703126‑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: ‑