Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,758,664	C→A	5.5%	Q312K (CAA→AAA)	hemBB →	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,758,664	0	C	A	5.5%	234.5 / 2.7	72	Q312K (CAA→AAA)	hemBB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base A (2/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

CCCTATCTGGACATCATTCAAAGGGTCAAGGCAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCACCCCTGCCCCCCATACCAACGAAGCC  >  NC_002947/3758526‑3758811                                                                                                                                           |                                                                                                                                                    cccTATCTGGACATCATTCAAAGGGTCAAGGCAAACTTTGCCGTGACGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGATGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAg                                                                                                                                          <  2:419796/150‑1 (MQ=255) cccTATCTGGACATCATTCAAAGGGTCAAGACAACATTTGCCGTGCCGACATTCGCCTACCATGTCAGCGGTGAATACTCGATGCTCAAGGCATCGATCAATCATGGCTGTGTGGATCCGCGCATGTTTGAACGCCATCAATGCCTcctg                                                                                                                                          >  2:208848/1‑148 (MQ=255)       ttGGACATCATGCAAATGGTCAAGGCAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAg                                                                                                                                    <  2:1327109/149‑1 (MQ=255)        tGGACATCATTCAAAGGGTCAAGGCAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGc                                                                                                                                   >  2:1244126/1‑150 (MQ=255)                atTCAAAGGGTCAAGGCAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGt                                                                                                                           >  2:434612/1‑150 (MQ=255)                 ttCAAAGGGTCAAGGCAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCAGGTTTGAACGCCATCAATGACTCCAGGATGAGCTGCCCGGTa                                                                                                                          >  2:311287/1‑150 (MQ=255)                 ttCAAAGGGTCAAGGCAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTa                                                                                                                          >  2:1266619/1‑150 (MQ=255)                  tCAAAGGGTCAAGGCAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTaa                                                                                                                         >  1:113550/1‑150 (MQ=255)                         gtaaaGGCAACCTTTGACGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACCTTGAAAGGCCATCAATGCCTCTAGGAGGGGCTACACGGTAAATGcgcc                                                                                                                  >  2:357152/4‑150 (MQ=255)                         gTCAAGGCAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGcgcc                                                                                                                  >  2:1177183/1‑150 (MQ=255)                         gTCAAGGCAACCTTTGCCGGGCCGACAGTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGAGCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGcgcc                                                                                                                  <  1:586095/150‑1 (MQ=255)                              ggCAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTg                                                                                                             >  2:938409/1‑150 (MQ=255)                                cAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATTAGCTGCACGGTAAATGCGCCGCCTGTa                                                                                                           <  1:357152/150‑1 (MQ=255)                                cAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTa                                                                                                           <  1:1177183/150‑1 (MQ=255)                                cAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTa                                                                                                           <  1:166210/150‑1 (MQ=255)                                  aCCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACt                                                                                                         >  1:1546527/1‑150 (MQ=255)                                    cTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATg                                                                                                                      >  1:431876/1‑135 (MQ=255)                                    cTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATg                                                                                                                      <  2:431876/135‑1 (MQ=255)                                        gccgtgccgACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGt                                                                                                   <  1:208848/150‑1 (MQ=255)                                         ccgtgccgACATTCGCCTACCAGGTCAGCGGTGAATACGCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTg                                                                                                  >  2:270362/1‑150 (MQ=255)                                          cgtgccgACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGa                                                                                                 <  2:1269639/150‑1 (MQ=255)                                              ccgACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCa                                                                                             >  1:1058361/1‑150 (MQ=255)                                              ccgACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCa                                                                                             >  1:434155/1‑150 (MQ=255)                                                gACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAgg                                                                                           >  1:1133837/1‑150 (MQ=255)                                                   aTTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGa                                                                                        <  1:167389/150‑1 (MQ=255)                                                   aTTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGa                                                                                        >  1:1438780/1‑150 (MQ=255)                                                   aTTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGa                                                                                        <  1:311287/150‑1 (MQ=255)                                                   aTTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGa                                                                                        <  1:1266619/150‑1 (MQ=255)                                                   aTTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGa                                                                                        <  2:1058361/150‑1 (MQ=255)                                                          tACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGa                                                                                 >  1:90884/1‑150 (MQ=255)                                                           aCCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAg                                                                                <  1:946819/150‑1 (MQ=255)                                                              aGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATAGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCag                                                                             <  2:1186814/150‑1 (MQ=255)                                                               ggTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCaga                                                                            <  1:189378/150‑1 (MQ=255)                                                                  cAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGcc                                                                         <  1:1244126/150‑1 (MQ=255)                                                                      ggTGAATACTCGATGCTCAAGGCAGCGATGAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCcgcac                                                                     <  2:113550/150‑1 (MQ=255)                                                                      ggTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCcgcac                                                                     >  1:1150609/1‑150 (MQ=255)                                                                               tCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGca                                                            >  1:1371393/1‑150 (MQ=255)                                                                                    gCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAg                                                                                                                                    >  2:462662/1‑73 (MQ=255)                                                                                    gCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAg                                                                                                                                    <  1:462662/73‑1 (MQ=255)                                                                                      tCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGa                                                                                     >  2:495438/1‑118 (MQ=255)                                                                                      tCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGa                                                                                     <  1:495438/118‑1 (MQ=255)                                                                                      tCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCt                                                     <  2:90884/150‑1 (MQ=255)                                                                                          ggCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTAt                                                 <  1:502021/150‑1 (MQ=255)                                                                                               cGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGa                                            <  2:1133837/150‑1 (MQ=255)                                                                                                   caatcaTGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCt                                        <  2:1150609/150‑1 (MQ=255)                                                                                                     atcaTGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTc                                      >  1:1301273/1‑150 (MQ=255)                                                                                                     atcaTGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTc                                      <  1:938409/150‑1 (MQ=255)                                                                                                           gctggctgGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATg                                <  2:434155/150‑1 (MQ=255)                                                                                                             tggctggATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAg                              <  1:1038637/150‑1 (MQ=255)                                                                                                                 tggATCCCCGCACGTTTGACCGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCAcc                          >  1:329302/1‑150 (MQ=255)                                                                                                                 tggATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCAcc                          >  2:117683/1‑150 (MQ=255)                                                                                                                    aTCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTACCGATGTGTACCCCt                       >  2:1274825/1‑150 (MQ=255)                                                                                                                     tCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATgcgc                                                                                                                   <  2:1032724/57‑1 (MQ=255)                                                                                                                     tCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATgcgc                                                                                                                   >  1:1032724/1‑57 (MQ=255)                                                                                                                       cccGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCa                                                                              >  1:1337604/1‑92 (MQ=255)                                                                                                                       cccGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCa                                                                              <  2:1337604/92‑1 (MQ=255)                                                                                                                        ccGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGa                                                                                                 >  2:1542058/1‑72 (MQ=255)                                                                                                                        ccGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGa                                                                                                 <  1:1542058/72‑1 (MQ=255)                                                                                                                        ccGCACGTTTGAACGCCATAAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGa                                                                                     <  1:1152358/84‑1 (MQ=255)                                                                                                                        ccGCACGTTTGAACGCCATAAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGa                                                                                     >  2:1152358/1‑84 (MQ=255)                                                                                                                        ccGCACGTTTGAACGCCATAAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCAc                           >  1:1082960/1‑142 (MQ=255)                                                                                                                        ccGCACGTTTGAACGCCATAAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCAc                           <  2:1082960/142‑1 (MQ=255)                                                                                                                         cGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGca                                                            <  2:206284/108‑1 (MQ=255)                                                                                                                         cGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGca                                                            >  1:206284/1‑108 (MQ=255)                                                                                                                         cGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCACCCCTGcccc                  <  1:434612/150‑1 (MQ=255)                                                                                                                          gCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCACCCCTGccccc                 >  2:975220/1‑150 (MQ=255)                                                                                                                              gTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCACCCCTGCCCCCCATa             <  2:1546527/150‑1 (MQ=255)                                                                                                                                      gCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTcc                                     <  1:1168042/118‑1 (MQ=255)                                                                                                                                      gCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTcc                                     >  2:1168042/1‑118 (MQ=255)                                                                                                                                      gCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCACCCCTGCCCCCCATACCAACGaa     >  1:1266699/1‑150 (MQ=255)                                                                                                                                       ccATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCACCCCTGCCCCCCATACCAACGAAg    <  2:329302/150‑1 (MQ=255)                                                                                                                                         aTCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCACCCCTGCCCCCCATACCAACGAAGcc  <  2:1438780/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                    CCCTATCTGGACATCATTCAAAGGGTCAAGGCAACCTTTGCCGTGCCGACATTCGCCTACCAGGTCAGCGGTGAATACTCGATGCTCAAGGCAGCGATCAATCATGGCTGGCTGGATCCCCGCACGTTTGAACGCCATCAATGCCTCCAGGATGAGCTGCACGGTAAATGCGCCGCCTGTACTCGCCGTGATGCAGGCGAGGACTGAGCAGAGCCGCACGCATGCGCACCACCCTTTATCATGAAGCTTCCCGATGAGCACCCCTGCCCCCCATACCAACGAAGCC  >  NC_002947/3758526‑3758811

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