Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,355,094:1	+C	13.4%	coding (250/681 nt)	modB →	molybdate ABC transporter permease ModB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,355,094	1	.	C	13.4%	218.4 / 33.9	67	coding (250/681 nt)	modB	molybdate ABC transporter permease ModB
Reads supporting (aligned to +/- strand):  ref base . (25/33);  new base C (5/4);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

TGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGA  >  NC_002947/4354947‑4355234                                                                                                                                                     |                                                                                                                                             tggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑Gc                                                                                                                                            >  2:1180940/1‑150 (MQ=255)  ggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑Gcc                                                                                                                                           <  2:386129/150‑1 (MQ=255)    cTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTg                                                                                                                                         <  1:1122735/150‑1 (MQ=255)     tCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTgg                                                                                                                                        >  1:465218/1‑150 (MQ=255)      cgcgcgCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAC‑GCCTggg                                                                                                                                       <  2:89199/150‑1 (MQ=255)       gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCGGCTGGTACTGCCACAGACGGTGCTAGGTATTTACGTGATGATCGCGCTTGGCCCGCAGGGCTGGCTTGGTCAGGCTACTCCAGCGCTAG‑GCCTGGGc                                                                                                                                      >  1:453182/1‑150 (MQ=255)           gcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAgtg                                                                                                                                  <  1:93823/150‑1 (MQ=255)              cgcgCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTgtgg                                                                                                                               <  2:336313/150‑1 (MQ=255)                    ctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCa                                                                                                                         <  2:430943/150‑1 (MQ=255)                            cgcgTGCCGGGGGGTGCGGTGTTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTcaca                                                                                                                 <  2:720824/150‑1 (MQ=255)                                   cGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTgg                                                                                                          <  2:1306784/150‑1 (MQ=255)                                         gtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTg                                                                                                    >  2:395670/1‑150 (MQ=255)                                          tgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTgg                                                                                                   <  1:1102287/150‑1 (MQ=255)                                            cGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCt                                                                                                 <  2:618169/150‑1 (MQ=255)                                            cGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCt                                                                                                 <  1:419303/150‑1 (MQ=255)                                             ggtggtGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTc                                                                                                <  1:1405605/150‑1 (MQ=255)                                              gtggtgGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCTGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCg                                                                                               >  2:470228/1‑150 (MQ=255)                                                ggtggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGAc                                                                                             <  2:1074372/150‑1 (MQ=255)                                                      aCTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCGCCCAGGC‑CTAGCGCCTgggtggcc                                                                                                                                  >  2:996131/1‑102 (MQ=255)                                                      aCTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCGCCCAGGC‑CTAGCGCCTgggtggcc                                                                                                                                  <  1:996131/107‑6 (MQ=255)                                                      aCTGCCGCTGGTACTGCCACCGACGGTGATCGGGTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTa                                                                                       >  1:1282363/1‑150 (MQ=255)                                                         gCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTc                                                                                    <  1:947167/150‑1 (MQ=255)                                                          ccGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGtt                                                                                                                           >  1:1312695/1‑110 (MQ=255)                                                          ccGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGtt                                                                                                                           <  2:1312695/110‑1 (MQ=255)                                                          ccGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGACTAGCTCGACACGGTGACCGTTTCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTcc                                                                                   >  2:692284/1‑150 (MQ=255)                                                            gCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCCCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGC‑CTAGCGCCTgggtg                                                                                                                                     >  1:1078752/1‑96 (MQ=255)                                                            gCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCCCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGC‑CTAGCGCCTgggtg                                                                                                                                     <  2:1078752/98‑3 (MQ=255)                                                             cTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATg                                                                                >  2:1137037/1‑150 (MQ=255)                                                                    tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTg                                                                                                                                         <  1:745455/86‑1 (MQ=255)                                                                    tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTg                                                                                                                                         >  2:745455/1‑86 (MQ=255)                                                                    tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCa                                                                                                                         <  1:191040/102‑1 (MQ=255)                                                                    tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCa                                                                                                                         >  2:191040/1‑102 (MQ=255)                                                                    tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCg                                                                         >  1:1077179/1‑150 (MQ=255)                                                                     gCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACGGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgt                                                                        >  1:1058957/1‑150 (MQ=255)                                                                      ccACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg                                                                       >  1:919352/1‑150 (MQ=255)                                                                            aCGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCaa                                                                 <  2:952010/150‑1 (MQ=255)                                                                            aCGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCaa                                                                 >  2:1278123/1‑150 (MQ=255)                                                                             cGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAAc                                                                <  1:152237/150‑1 (MQ=255)                                                                              ggTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGATGTACTCCATGCCGTTCGTGGTGaaaaa                                                               >  1:1003106/1‑148 (MQ=255)                                                                               gTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAAccc                                                              <  1:496524/150‑1 (MQ=255)                                                                               gTGATCGGTTTTTACCTGTTGATCGCCCCTGGCCCGCACGGCTGGCTTGCTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTCGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAAccc                                                              <  1:519706/150‑1 (MQ=255)                                                                                   tCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGc                                                                   >  1:200712/1‑141 (MQ=255)                                                                                   tCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGc                                                                   <  2:200712/141‑1 (MQ=255)                                                                                    cGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGATCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCa                                                         >  2:1171075/1‑150 (MQ=255)                                                                                        ttttACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAAt                                                     <  2:465218/150‑1 (MQ=255)                                                                                                    aTCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAtt                                                                                                     >  2:352216/1‑90 (MQ=255)                                                                                                    aTCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAtt                                                                                                     <  1:352216/90‑1 (MQ=255)                                                                                                        cccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAt                                                                                                      <  2:148983/85‑1 (MQ=255)                                                                                                        cccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAt                                                                                                      >  1:148983/1‑85 (MQ=255)                                                                                                                     cGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCt                        <  2:1282363/150‑1 (MQ=255)                                                                                                                          ggctTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAg                   <  2:191702/150‑1 (MQ=255)                                                                                                                               ggTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgcc              <  1:1087689/150‑1 (MQ=255)                                                                                                                                gTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACGGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTAGAACTACCGGGTCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccg             <  2:1058957/150‑1 (MQ=255)                                                                                                                                 tCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccg             <  2:1420483/149‑1 (MQ=255)                                                                                                                                 tCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccg             >  1:1420483/1‑149 (MQ=255)                                                                                                                                     gCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTAct        >  2:1228973/1‑150 (MQ=255)                                                                                                                                     gCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTAct        >  1:1185677/1‑150 (MQ=255)                                                                                                                                     actgcccAGGC‑CTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTg                       >  1:1136443/5‑135 (MQ=255)                                                                                                                                     actgcccAGGC‑CTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTg                       <  2:1136443/131‑1 (MQ=255)                                                                                                                                     actgcccAGGC‑CTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTAct        >  2:475346/5‑150 (MQ=255)                                                                                                                                     actgcccAGGC‑CTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTAct        >  2:1391474/5‑150 (MQ=255)                                                                                                                                       tgcccAGGC‑CTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTActct      <  1:475346/148‑1 (MQ=255)                                                                                                                                        aCCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCATGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTg     <  2:919352/150‑1 (MQ=255)                                                                                                                                           cAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACCCTGCGa  >  2:224190/1‑150 (MQ=255)                                                                                                                                           cAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTAGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGa  <  2:453182/150‑1 (MQ=255)                                                                                                                                            aGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAtt                                                                                                     <  2:560857/50‑1 (MQ=255)                                                                                                                                            aGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAtt                                                                                                     >  1:560857/1‑50 (MQ=255)                                                                                                                                                     |                                                                                                                                             TGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGA  >  NC_002947/4354947‑4355234

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