Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,235,837	G→C	6.8%	P316P (CCG→CCC)	uvrB →	uvrABC excinuclease protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,235,837	0	G	C	6.8%	188.6 / 3.9	59	P316P (CCG→CCC)	uvrB	uvrABC excinuclease protein B
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base C (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

GAGCGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTCC  >  NC_002947/2235691‑2235985                                                                                                                                                     |                                                                                                                                                      gAGCGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCccggc                                                                                                                                                     >  6:126245/1‑149 (MQ=255)    cGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCccggccgg                                                                                                                                                  <  2:222710/149‑1 (MQ=255)     gTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCggg                                                                                                                                                 <  8:94040/149‑1 (MQ=255)     gTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCggg                                                                                                                                                 <  2:7782/149‑1 (MQ=255)      tCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCgggg                                                                                                                                                >  2:99427/1‑149 (MQ=255)          gAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCCGGCCGGGcgcc                                                                                                                                            <  7:320799/149‑1 (MQ=255)             tACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGcgccgcc                                                                                                                                         >  4:39523/1‑149 (MQ=255)                  gCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGgc                                                                                                                                               <  2:59506/138‑1 (MQ=255)                  gCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGgc                                                                                                                                               >  1:59506/1‑138 (MQ=255)                  gCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCa                                                                                                                                    >  2:237865/1‑149 (MQ=255)                   cacaAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGccc                                                                                                                                     <  5:126245/147‑1 (MQ=255)                      aaGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCt                                                                                                                                >  8:29711/1‑149 (MQ=255)                        ggCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTg                                                                                                                               >  2:259464/1‑148 (MQ=255)                            aacaaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTAc                                                                                                                            >  1:90531/1‑148 (MQ=255)                            aacaaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGa                                                                                                                          >  2:229440/1‑149 (MQ=255)                             acaaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGcgccgc                                                                                                                                          <  4:262329/132‑1 (MQ=255)                             acaaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGcgccgc                                                                                                                                          >  3:262329/1‑132 (MQ=255)                              caaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGAtt                                                                                                                        <  8:16505/149‑1 (MQ=255)                                aGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTAc                                                                                                                      <  5:11555/149‑1 (MQ=255)                                   tGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTg                                                                                                                   >  7:6142/1‑149 (MQ=255)                                        gAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTgccgc                                                                                                               <  5:321704/148‑1 (MQ=255)                                           gcgcAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGa                                                                                                           <  1:202951/149‑1 (MQ=255)                                             gcAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATg                                                                                                         <  6:150787/149‑1 (MQ=255)                                               aGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCATGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATgc                                                                                                        <  1:99427/149‑1 (MQ=255)                                                       caGCAACGCTCACGCTTCGATCTGG‑AGATGATCCTCGATCTGGGGTAGTGCAACGGCATCGAAAACTACTAGCGTTACC‑TGTCCGGGCTCCCGGCCGGGGCGCCGCCGCCCACCCTGTTCGATTACCTGCCGCCCGATGCGct‑gctg                                                                                                <  4:252978/146‑1 (MQ=255)                                                        aGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTACAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGt                                                                                              >  8:177261/1‑148 (MQ=255)                                                             cGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATcga                                                                                        <  2:90531/149‑1 (MQ=255)                                                                  ccGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATgcgc                                                                                                      <  6:177792/131‑1 (MQ=255)                                                                  ccGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATgcgc                                                                                                      >  5:177792/1‑131 (MQ=255)                                                                         gATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGt                                                                            <  8:6142/149‑1 (MQ=255)                                                                               g‑agaTGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCa                                                                          <  1:237865/145‑1 (MQ=255)                                                                                        ccTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCGTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAg                                                               >  1:76562/1‑149 (MQ=255)                                                                                        ccTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAg                                                               <  2:76562/148‑1 (MQ=255)                                                                                         cTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGt                                                             <  1:229440/149‑1 (MQ=255)                                                                                          tCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTc                                                            >  5:174438/1‑149 (MQ=255)                                                                                          tCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCCGGCCGGGcgccc                                                                                                                                           <  3:199053/71‑2 (MQ=39)                                                                                          tCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCCGGCCGGGcgccc                                                                                                                                           >  4:199053/1‑70 (MQ=39)                                                                                            gAGCTGGGCTACTACAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTc                                                            <  7:177261/147‑1 (MQ=255)                                                                                                   gCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGct‑gctg                                                                                                <  1:304154/104‑1 (MQ=255)                                                                                                   gCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGct‑gctg                                                                                                >  2:304154/1‑104 (MQ=255)                                                                                                    cTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCAGCGCCATGt                                                  <  7:63199/149‑1 (MQ=255)                                                                                                        tGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTAc                                                >  8:240886/1‑147 (MQ=255)                                                                                                          cAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAgg                                            <  7:240886/149‑1 (MQ=255)                                                                                                               gCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGAc                                       >  4:294566/1‑149 (MQ=255)                                                                                                               gCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGAc                                       >  5:133741/1‑149 (MQ=255)                                                                                                                  tCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGt                                    >  4:30147/1‑149 (MQ=255)                                                                                                                  tCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGt                                    >  3:149663/1‑149 (MQ=255)                                                                                                                     aaaaCTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCa                                 <  4:149663/149‑1 (MQ=255)                                                                                                                     aaaaCTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCa                                 >  7:43393/1‑149 (MQ=255)                                                                                                                            ctCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAg                           <  4:254420/148‑1 (MQ=255)                                                                                                                            ctCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAg                           <  4:258837/148‑1 (MQ=255)                                                                                                                               gcgTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGgaga                       <  3:30147/149‑1 (MQ=255)                                                                                                                                    aCC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTg                  >  1:25125/1‑149 (MQ=255)                                                                                                                                         gcggcggCGCCCCGGCCGGGCGCCGCCGCCCCCCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGCCGAATCACACGTCAGCGGTCCGCCCGTCGGGGCCATGTACAAGGGCGCCCGTTCACGCAAGGAGACCCTGGtgg              >  4:331037/6‑147 (MQ=255)                                                                                                                                          tCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTc               >  7:300674/1‑147 (MQ=255)                                                                                                                                          tCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGa             <  3:39523/149‑1 (MQ=255)                                                                                                                                             gggCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTa          >  6:252998/1‑149 (MQ=255)                                                                                                                                              ggCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTAc         <  6:174438/149‑1 (MQ=255)                                                                                                                                              ggCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTAc         >  4:267894/1‑149 (MQ=255)                                                                                                                                                 gcCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACgg       >  1:95044/1‑149 (MQ=255)                                                                                                                                                 gcCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAg                           >  3:109366/1‑128 (MQ=255)                                                                                                                                                    cggccggGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTc   <  7:29711/149‑1 (MQ=255)                                                                                                                                                     ggccggGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTcc  >  3:336359/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GAGCGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGG‑AGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACC‑TGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCT‑GCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTCC  >  NC_002947/2235691‑2235985

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