Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	152,158	A→T	10.6%	intergenic (‑52/+42)	PP_0143 ← / ← PP_0144	membrane protein/insulinase family metalloprotease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	152,158	0	A	T	10.6%	177.6 / 10.3	66	intergenic (‑52/+42)	PP_0143/PP_0144	membrane protein/insulinase family metalloprotease
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base T (4/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTA  >  NC_002947/152012‑152304                                                                                                                                                   |                                                                                                                                                   ccGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGcc                                                                                                                                                 >  1:1124569/1‑150 (MQ=255)     tACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCcgcg                                                                                                                                             >  1:342562/1‑150 (MQ=255)         ggatagGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaaga                                                                                                                                         <  1:1621870/146‑3 (MQ=255)         ggAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAg                                                                                                                                         <  2:1285480/150‑1 (MQ=255)         ggAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAg                                                                                                                                         <  1:178061/150‑1 (MQ=255)          gagaGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGGAAgg                                                                                                                                        >  2:854248/1‑150 (MQ=255)          gagaGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAgg                                                                                                                                        <  1:770636/150‑1 (MQ=255)               gAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCgg                                                                                                                                   <  1:237854/150‑1 (MQ=255)                    cgccgAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTg                                                                                                                              <  1:109621/150‑1 (MQ=255)                     gccgAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTCGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAACGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTgg                                                                                                                             <  2:1004548/150‑1 (MQ=255)                        gAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCaa                                                                                                                          >  1:1813874/1‑150 (MQ=255)                          gcgcAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCaaaa                                                                                                                        <  1:1597338/150‑1 (MQ=255)                           cgcAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAg                                                                                                                       <  1:1777984/150‑1 (MQ=255)                               atgatgCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATa                                                                                                                   <  1:1072540/150‑1 (MQ=255)                                          gTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTa                                                                                                        <  2:1115797/150‑1 (MQ=255)                                          gTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTa                                                                                                        >  2:241221/1‑150 (MQ=255)                                             gcgcTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCg                                                                                                     <  1:30491/150‑1 (MQ=255)                                                         ttCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGc                                                                                         >  1:1159739/1‑150 (MQ=255)                                                         ttCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGAAGAGCGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGc                                                                                         >  1:593046/1‑150 (MQ=255)                                                         ttCGATGATGCGCTGGACGATATGCGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGc                                                                                         >  1:1484759/1‑150 (MQ=255)                                                          tCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCGCTTCGCGGGCAACCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGca                                                                                        <  1:909032/150‑1 (MQ=255)                                                              tgatgCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcaccac                                                                                    >  2:178008/1‑150 (MQ=255)                                                                 tgCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTgg                                                                                 <  2:517593/150‑1 (MQ=255)                                                                          cGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAg                                                                        <  1:710512/150‑1 (MQ=255)                                                                          cGATACGGGGGGTTTGGGGCACTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAg                                                                        <  1:553930/150‑1 (MQ=255)                                                                              atGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTcc                                                                    <  2:693288/150‑1 (MQ=255)                                                                                gggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCAt                                                                  <  1:1235471/150‑1 (MQ=255)                                                                                 ggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTc                                                                     >  2:1437503/1‑146 (MQ=255)                                                                                 ggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTc                                                                     <  1:1437503/146‑1 (MQ=255)                                                                                 ggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCAGTACTGGCAAAAGCATACCGCCGATTTACCGCTGCCGTGCAGCACCACTGGGCGGTTCAGGTCCATg                                                                 >  2:1031602/1‑150 (MQ=255)                                                                                  gggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCgtg                                                                                             >  2:1292331/1‑121 (MQ=255)                                                                                  gggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCgtg                                                                                             <  1:1292331/121‑1 (MQ=255)                                                                                  gggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATgc                                                                <  1:220067/150‑1 (MQ=255)                                                                                  gggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATgc                                                                >  1:1477679/1‑150 (MQ=255)                                                                                      tttGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGtt                                                                           <  2:1112625/135‑1 (MQ=255)                                                                                      tttGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGtt                                                                           >  1:1112625/1‑135 (MQ=255)                                                                                          gggCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCa                                                        >  1:1634311/1‑150 (MQ=255)                                                                                            gCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAgc                                                      >  1:582160/1‑150 (MQ=255)                                                                                            gCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGCGCAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAgc                                                      >  1:398493/1‑150 (MQ=255)                                                                                                  ggggCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACa                                                >  1:1896809/1‑150 (MQ=255)                                                                                                                   gTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATAccgcc                                                                                                              >  1:284897/1‑71 (MQ=255)                                                                                                                   gTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATAccgcc                                                                                                              <  2:284897/71‑1 (MQ=255)                                                                                                                   gTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGccc                               <  2:1896809/150‑1 (MQ=255)                                                                                                                         tGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGTACCACTGGGCCGTTCAGGTCCACGCCGCGGCCAGCGCAACAGCCCCGCCACCACGCCCAGcccc                         >  1:28614/1‑150 (MQ=255)                                                                                                                             tgtgCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGAGTAGGTCGTGGGAATTTCATACCGCCGATTTACCGCTGCTTTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCACCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGa                     >  1:798785/1‑150 (MQ=255)                                                                                                                               tgCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCGGGCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGc                   >  1:1452669/1‑150 (MQ=255)                                                                                                                                 ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGa                                                                                                            <  2:938860/59‑1 (MQ=255)                                                                                                                                 ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGa                                                                                                            >  1:938860/1‑59 (MQ=255)                                                                                                                                 ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCaa                 >  2:673627/1‑150 (MQ=255)                                                                                                                                  cAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCacg                >  2:131649/1‑148 (MQ=255)                                                                                                                                  cAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAg                <  2:1835711/150‑1 (MQ=255)                                                                                                                                  cAGCCTCTTCCCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAg                <  1:854248/150‑1 (MQ=255)                                                                                                                                  aacccTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAg                <  2:897262/147‑1 (MQ=255)                                                                                                                                   gaccctttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAg                             <  2:1324135/130‑1 (MQ=255)                                                                                                                                   gaccctttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAg                             >  1:1324135/7‑136 (MQ=255)                                                                                                                                    gCCTCTTCGCGGGCAGGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGc              >  1:1139842/1‑150 (MQ=255)                                                                                                                                       ctttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCAt           >  1:1016279/3‑150 (MQ=255)                                                                                                                                        tttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGcccc                                           <  1:177679/116‑1 (MQ=255)                                                                                                                                        tttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGcccc                                           >  2:177679/2‑117 (MQ=255)                                                                                                                                        tttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATc          >  2:678851/2‑150 (MQ=255)                                                                                                                                        ctTCGCGGGCAAGCCCGCGAAAGGGTGGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATc          <  2:593046/150‑1 (MQ=255)                                                                                                                                        ctTCGCGGGCAAGCCCGCGAAAGGGTCGGGACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATc          >  1:815950/1‑150 (MQ=255)                                                                                                                                         ttCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGcccc                                           <  2:260981/116‑1 (MQ=255)                                                                                                                                         ttCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGcccc                                           >  1:260981/1‑116 (MQ=255)                                                                                                                                         ttCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCa         >  1:1030118/1‑150 (MQ=255)                                                                                                                                                gCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTa  <  2:582160/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTA  >  NC_002947/152012‑152304

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