Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,255,849	C→A	6.3%	intergenic (‑112/+24)	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,849	0	C	A	6.3%	196.1 / 3.7	64	intergenic (‑112/+24)	PP_2852/PP_2853	sulfatase domain‑containing protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base A (2/2);  total (30/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 = 3.65e-01

AAGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGAA  >  NC_002947/3255702‑3255995                                                                                                                                                     |                                                                                                                                                   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                                                                                                 >  5:25981/1‑148 (MQ=255)                                                    gCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTg                                                                                                 >  8:88732/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                                                                                                               <  6:104425/125‑1 (MQ=255)                                                             tCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                               >  5:104425/1‑125 (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                                                                                  >  6:29211/1‑144 (MQ=255)                                                                       tCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGt                                                                                  <  5:29211/144‑1 (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                                                                                              <  8:1971/116‑1 (MQ=255)                                                                                       tC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                              >  7:1971/1‑116 (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                   <  5:208700/149‑1 (MQ=255)                                                                                                                                  ccAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAg                   <  7:330170/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)                                                                                                                                                    gCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGaa  <  1:201812/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AAGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATC‑AGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGAA  >  NC_002947/3255702‑3255995

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