Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,615,347	T→C	9.2%	I573T (ATT→ACT)	PP_4085 →	rhs family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,615,347	0	T	C	9.2%	172.7 / 6.3	65	I573T (ATT→ACT)	PP_4085	rhs family protein
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base C (2/4);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCAGG  >  NC_002947/4615201‑4615493                                                                                                                                                   |                                                                                                                                                   tCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTc                                                                                                                                                  <  4:282591/149‑1 (MQ=255) tCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAAACCGATCACGAATGGCGAGTGATTc                                                                                                                                                  >  1:137625/1‑149 (MQ=255)        gCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTg                                                                                                                                           >  3:196347/1‑149 (MQ=255)         cTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTgg                                                                                                                                          <  4:196347/149‑1 (MQ=255)         cTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTgg                                                                                                                                          <  3:239886/149‑1 (MQ=255)            gCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAATCACTCGCCattcgt                                                                                                                                         <  2:43916/147‑6 (MQ=255)               gTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGc                                                                                                                                    <  6:361286/149‑1 (MQ=255)                tGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCg                                                                                                                                   <  4:61355/149‑1 (MQ=255)                     tCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATgg                                                                                                                               <  4:324108/148‑1 (MQ=255)                     tCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATg                                                                                                                                <  3:138536/147‑1 (MQ=255)                        gcgAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAATCACTCGCCattc                                                                                                                                           <  4:309183/133‑4 (MQ=255)                               cGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTACACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTAcc                                                                                                                    <  3:41330/149‑1 (MQ=255)                                     ggCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCACCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTc                                                                                                               <  6:273870/148‑1 (MQ=255)                                      gcgcgcTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTc                                                                                                               >  7:216436/1‑147 (MQ=255)                                            ttCGCCATTCCTGAAACGCGCCCTATTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACaa                                                                                                       >  5:84983/1‑149 (MQ=255)                                            ttCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAGTGGCGAATCACTCGCCattc                                                                                                                                           >  2:340795/1‑110 (MQ=255)                                            ttCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAGTGGCGAATCACTCGCCattc                                                                                                                                           <  1:340795/113‑4 (MQ=255)                                               gCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTAc                                                                                                                     >  6:130903/1‑132 (MQ=255)                                               gCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTAc                                                                                                                     <  5:130903/132‑1 (MQ=255)                                                  aTTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTgg                                                                                                  >  8:140710/1‑148 (MQ=255)                                                       tGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCg                                                                                            >  8:113268/1‑149 (MQ=255)                                                       tGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCGGGAAGCCg                                                                                            >  4:291341/1‑149 (MQ=255)                                                                cccTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTAc                                                                                                         <  3:189066/127‑1 (MQ=255)                                                                cccTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTAc                                                                                                         >  4:189066/1‑127 (MQ=255)                                                                  cTATTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTc                                                                                 <  6:84983/149‑1 (MQ=255)                                                                      tCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCa                                                                             >  3:280770/1‑149 (MQ=255)                                                                        ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACg                                                                           >  5:69380/1‑149 (MQ=255)                                                                        ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACg                                                                           >  7:146581/1‑149 (MQ=255)                                                                        ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACg                                                                           >  3:79864/1‑149 (MQ=255)                                                                         cAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACg                                                                           >  5:44148/1‑148 (MQ=255)                                                                          aaGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGa                                                                         <  7:140710/149‑1 (MQ=255)                                                                              gcCTGGACGGCACGCCGTACCGCTTGCCTCCACTGCTGGAGCCCCAACGGGATCACGAATGGCGAGTGATTGGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACTGACAAc                                                                     <  3:163284/149‑1 (MQ=255)                                                                              gcCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAAc                                                                     >  3:285795/1‑149 (MQ=255)                                                                               ccTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTGCCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAAc                                                                     >  4:75013/1‑148 (MQ=255)                                                                                cTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGaa                                                                                                >  2:151256/1‑120 (MQ=255)                                                                                cTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGaa                                                                                                <  1:151256/120‑1 (MQ=255)                                                                                 tGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGaa                                                                                                                           >  6:109227/1‑92 (MQ=255)                                                                                 tGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGaa                                                                                                                           <  5:109227/92‑1 (MQ=255)                                                                                            cccTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGa                                                       <  4:285795/149‑1 (MQ=255)                                                                                               tACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCgg                                                                                           <  8:116914/110‑1 (MQ=255)                                                                                               tACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCgg                                                                                           >  7:116914/1‑110 (MQ=255)                                                                                                  cGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCacta                                                 >  4:78217/1‑149 (MQ=255)                                                                                                        ccTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCactactact                                            <  4:79864/148‑1 (MQ=255)                                                                                                         cTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTgg                                          <  6:44148/149‑1 (MQ=255)                                                                                                          tCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTgg                                          >  4:23341/1‑148 (MQ=255)                                                                                                            cACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGtt                                                                                                                >  8:297726/1‑76 (MQ=255)                                                                                                            cACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAAccc                                                                                   <  5:121501/105‑1 (MQ=255)                                                                                                            cACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAAccc                                                                                   >  6:121501/1‑105 (MQ=255)                                                                                                            cACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTggggg                                       <  6:69380/149‑1 (MQ=255)                                                                                                            cACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTggg                                         <  8:216436/147‑1 (MQ=255)                                                                                                                   ggAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTg                                 >  8:212654/1‑148 (MQ=255)                                                                                                                          cAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCgg                                                                                                                  <  5:37884/60‑1 (MQ=255)                                                                                                                          cAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCgg                                                                                                                  >  6:37884/1‑60 (MQ=255)                                                                                                                          cAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGt                         <  4:280770/149‑1 (MQ=255)                                                                                                                            aCCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTTAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCt                       <  7:113268/149‑1 (MQ=255)                                                                                                                            aCCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCt                       <  7:212654/149‑1 (MQ=255)                                                                                                                              ccGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGAGCCACTGTACGAGGTCtat                     <  3:23341/149‑1 (MQ=255)                                                                                                                                    aCGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCactactac                                             <  6:57732/119‑1 (MQ=255)                                                                                                                                    aCGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCactactac                                             >  5:57732/1‑119 (MQ=255)                                                                                                                                     ccagtGGCGAATCACTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGcac               >  8:185023/5‑148 (MQ=255)                                                                                                                                      gAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCATCGTCACGGACAACCTCGGCCGCCGTGATCAGTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACAt             >  3:112059/1‑149 (MQ=255)                                                                                                                                          ggCGAATCACTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGa          <  7:185023/148‑1 (MQ=255)                                                                                                                                           gCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGAc         <  3:75013/148‑1 (MQ=255)                                                                                                                                              aGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCt      >  6:262775/1‑148 (MQ=255)                                                                                                                                                  aTTCGCCAATTGGGCAGGGATGGTGAAGTGTACCGGTTCGAGTACAACCTGGAAGCCGGTGTAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCAgg  <  4:112059/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCAGG  >  NC_002947/4615201‑4615493

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