Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,255,846	T→G	6.2%	intergenic (‑109/+27)	PP_2852 ← / ← PP_2853	sulfatase domain‑containing protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,255,846	0	T	G	6.2%	203.3 / 4.0	65	intergenic (‑109/+27)	PP_2852/PP_2853	sulfatase domain‑containing protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base G (2/2);  total (31/34)
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.79e-01

GGCAAGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATA  >  NC_002947/3255699‑3255990                                                                                                                                                     |                                                                                                                                                 cccAAGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCGAAGCCAGGCCCGTTc                                                                                                                                                 <  8:122393/147‑1 (MQ=255)    aaGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGca                                                                                                                                              >  6:228413/1‑149 (MQ=255)                  ggtggtCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCTGTGCGAACGGGCCCGGctttg                                                                                                                                <  2:205471/149‑5 (MQ=255)                  ggtggtCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGAt                                                                                                                                <  2:326635/149‑1 (MQ=255)                   gtggtCGGACACGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATA‑AGGACAGCGCAGGGAAGGCGGCCAGTTACAAAGCTTGTAACCAAAGCCGGGCCCCTTCCCACAGGCCCGGCCGATg                                                                                                                               >  8:301095/1‑149 (MQ=255)                           aCAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTAt                                                                                                                       >  8:176228/1‑149 (MQ=255)                               gTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAg                                                                                                                    >  2:170309/1‑148 (MQ=255)                               gTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATA‑AGGCCAGGGCAGGGAAGGAGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAg                                                                                                                    >  7:320983/1‑148 (MQ=255)                                ttcttcACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGaa                                                                                                                  <  4:27186/149‑1 (MQ=255)                                     cacaAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGt                                                                                                                <  7:355783/146‑1 (MQ=255)                                      acaAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaat                                                                                                             >  5:227083/1‑148 (MQ=255)                                       caAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaataa                                                                                                           <  2:344022/149‑1 (MQ=255)                                       caAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaat                                                                                                             >  1:344022/1‑147 (MQ=255)                                        aaaaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGATGGCGGCCAGTTACAAGGATTGTAAGCAAATCCGGGCCCGTTCGCTATGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                          <  7:301095/149‑1 (MQ=255)                                        aaaaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                          >  5:127806/1‑149 (MQ=255)                                        aaaaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGACAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAAATATCAGAAGTAATAAg                                                                                                          >  5:183837/1‑149 (MQ=255)                                          aaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTAt                                                                                                                       >  3:163172/1‑134 (MQ=255)                                          aaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTAt                                                                                                                       <  4:163172/134‑1 (MQ=255)                                                 ccTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAAttt                                                                                                   >  3:319539/1‑147 (MQ=255)                                                       gCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTg                                                                                            >  8:88732/1‑148 (MQ=255)                                                       gCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTg                                                                                            >  5:25981/1‑148 (MQ=255)                                                         cATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                         <  7:109309/149‑1 (MQ=255)                                                         cATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                         <  5:201475/149‑1 (MQ=255)                                                          aTTAAATCGTCATGCATCTGTGCTTGATGATATCGAGGCCAGGGCAGGTAAGGCGCCCAGTTACAAGGATTGTAACAAAAGCCGGGCCCTTTCGCACAGGCCCGCCCGATGGAAACTATCAGAAGAAATAAGGGAATTACAGGCTgaa                                                                                          >  1:243046/1‑148 (MQ=255)                                                                tCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                          >  5:104425/1‑125 (MQ=255)                                                                tCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                          <  6:104425/125‑1 (MQ=255)                                                                  gTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCa                                                                                <  4:319539/149‑1 (MQ=255)                                                                     aTGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAgg                                                                              <  5:228413/148‑1 (MQ=255)                                                                      tGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTg                                                                            <  7:88732/149‑1 (MQ=255)                                                                        cATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCg                                                                          <  7:176228/149‑1 (MQ=255)                                                                         aTCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTTCAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGt                                                                         <  6:127806/149‑1 (MQ=255)                                                                         aTCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGt                                                                         >  5:361689/1‑149 (MQ=255)                                                                          tCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGt                                                                             <  5:29211/144‑1 (MQ=255)                                                                          tCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGt                                                                             >  6:29211/1‑144 (MQ=255)                                                                                   gatgatATC‑TGGCCAGGGCATGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTc                                                                 <  8:320983/147‑1 (MQ=255)                                                                                    atgatATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGa                                                              >  5:195576/1‑149 (MQ=255)                                                                                         atC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAg                                                                               <  5:60453/127‑1 (MQ=255)                                                                                         atC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAg                                                                               >  6:60453/1‑127 (MQ=255)                                                                                          tC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                         >  7:1971/1‑116 (MQ=255)                                                                                          tC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                         <  8:1971/116‑1 (MQ=255)                                                                                          tC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATg                                                        <  6:25981/149‑1 (MQ=255)                                                                                             aGGCCAGGGCAGTGAAGGCTGCCAGTTACTAGGCTTGTAACCAAAGCAGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAgg                                                                                               <  2:158170/108‑1 (MQ=255)                                                                                             aGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAgg                                                                                               >  1:158170/1‑108 (MQ=255)                                                                                                 cagggcagggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAgg                                                  >  4:374349/1‑149 (MQ=255)                                                                                                    ggcagggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTg                                               >  8:215600/1‑149 (MQ=255)                                                                                                     gcagggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTgg                                              <  6:227083/149‑1 (MQ=255)                                                                                                         ggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCAc                                          >  6:131853/1‑149 (MQ=255)                                                                                                          gAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCAc                                          <  5:131853/148‑1 (MQ=255)                                                                                                           aaGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGa                                        <  1:170309/149‑1 (MQ=255)                                                                                                             ggcggcCAGTTACAAGGGTTGTAACCAAAGCCGGGCGCGTTAGCTCATGCCCGGCCGATGGCAACTATCAGAAGTAATATGGTAATTTCAGGCTGAAAGAAAAGTCAGGGGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGAtt                                      <  6:183837/149‑2 (MQ=255)                                                                                                             ggcggcCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATg                                      <  1:63667/149‑1 (MQ=255)                                                                                                                 tttggTTACAAGCCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGt                                    >  1:378742/5‑147 (MQ=255)                                                                                                                   cAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAgg                                <  6:361689/149‑1 (MQ=255)                                                                                                                   cAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGATGGGCACGATGGTCAgg                                >  6:208700/1‑149 (MQ=255)                                                                                                                           aGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGt                        >  2:201812/1‑149 (MQ=255)                                                                                                                             ccTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTgg                      <  2:378742/148‑1 (MQ=255)                                                                                                                              cTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTggc                     >  1:161204/1‑149 (MQ=255)                                                                                                                                 gTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTggcggc                  >  8:330170/1‑149 (MQ=255)                                                                                                                                     ccAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAg              <  7:330170/149‑1 (MQ=255)                                                                                                                                     ccAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAg              <  5:208700/149‑1 (MQ=255)                                                                                                                                        aaGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTc           <  6:195576/149‑1 (MQ=255)                                                                                                                                            cGGGCCTGTGCGAACGGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGt                               <  2:14971/125‑1 (MQ=255)                                                                                                                                            cGGGCCTGTGCGAACGGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGt                               >  1:14971/1‑125 (MQ=255)                                                                                                                                              ggCCTGTGCGAACGGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGcc      >  6:319661/1‑148 (MQ=255)                                                                                                                                                 ccGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAata  <  2:243046/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGCAAGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATA  >  NC_002947/3255699‑3255990

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑