Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,054,242	Δ1 bp	7.9%	coding (772/1320 nt)	PP_3573 →	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,054,242	0	A	.	7.9%	186.4 / 11.0	63	coding (772/1320 nt)	PP_3573	monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base . (3/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

AACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGAC  >  NC_002947/4054095‑4054385                                                                                                                                                      |                                                                                                                                                aaCCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGTCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGAcc                                                                                                                                               >  2:1040652/1‑150 (MQ=255)   ccGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCt                                                                                                                                             >  1:1145599/1‑150 (MQ=255)    cgcgACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCAACA‑GGG‑CATGACCCTg                                                                                                                                            <  1:307962/150‑1 (MQ=255)     gcgACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTgg                                                                                                                                           <  2:1272812/150‑1 (MQ=255)       gACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATAGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGGTCATG‑CCCTGGtg                                                                                                                                         <  2:1093897/150‑2 (MQ=255)          ttCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGc                                                                                                                                      <  1:1144304/150‑1 (MQ=255)            ctgctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCt                                                                                                                                    <  2:1211471/150‑1 (MQ=255)               ctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTACGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGAc                                                                                                                                 >  2:879048/1‑150 (MQ=255)                    ggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGTCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGt                                                                                                                            <  1:1040652/150‑1 (MQ=255)                     gCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTc                                                                                                                           >  2:311027/1‑150 (MQ=255)                         ggCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTc                                                                                                                       <  2:307488/150‑1 (MQ=255)                                           gACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCgcgcgc                                                                                                     >  1:1160366/1‑150 (MQ=255)                                            aCGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCt                                                                                                    >  2:165093/1‑150 (MQ=255)                                                 gAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGGTCATG‑CCCTGGTGGGc                                                                                                                                      <  2:1014586/111‑1 (MQ=255)                                                 gAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGGTCATG‑CCCTGGTGGGc                                                                                                                                      >  1:1014586/1‑111 (MQ=255)                                                    aTTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGA‑CCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACa                                                                                           >  1:353912/1‑150 (MQ=255)                                                        ccAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAcc                                                                                        <  2:672464/150‑1 (MQ=255)                                                        ccAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAcc                                                                                        >  2:412896/1‑150 (MQ=255)                                                          aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTg                                                                                      <  2:663520/150‑1 (MQ=255)                                                          aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTg                                                                                      >  1:62562/1‑150 (MQ=255)                                                          aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACAGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCa                                                                                                >  1:1049612/1‑140 (MQ=255)                                                          aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACAGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCa                                                                                                <  2:1049612/140‑1 (MQ=255)                                                           aGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTgg                                                                                     <  1:625696/150‑1 (MQ=255)                                                             cccGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTc                                                                                   <  2:1175774/150‑1 (MQ=255)                                                              ccGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGA‑CCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCga                                                                                 <  2:353912/150‑1 (MQ=255)                                                              ccGGTCGCGAGCATGCGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGGCGCGCGCTTCCAGGACAACCTGGTCg                                                                                  >  1:258525/1‑150 (MQ=255)                                                               cGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCga                                                                                 >  2:1117248/1‑150 (MQ=255)                                                               cGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCga                                                                                 <  2:62562/150‑1 (MQ=255)                                                                 gTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGGCTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCATTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCgaga                                                                               <  2:931000/150‑1 (MQ=255)                                                                 gTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCACGCGCTTCCAGGACAACCTGGTCgaga                                                                               >  1:1261717/1‑150 (MQ=255)                                                                  tCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGaa                                                                              <  2:630635/150‑1 (MQ=255)                                                                      gAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATc                                                                          >  1:133324/1‑150 (MQ=255)                                                                           tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcg                                                                     >  2:424106/1‑150 (MQ=255)                                                                           tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcg                                                                     <  1:112661/150‑1 (MQ=255)                                                                                 cATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcga                                                               >  2:867095/1‑150 (MQ=255)                                                                                  aTCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGTACAACCTGGTCGAGAACATCAATCGCGGcgac                                                              >  2:565958/1‑150 (MQ=255)                                                                                       ggTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACaa                                                         >  1:597768/1‑150 (MQ=255)                                                                                            gTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAcc                                                    <  1:412896/150‑1 (MQ=255)                                                                                                      tggtggCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgct                                          <  2:76607/150‑1 (MQ=255)                                                                                                        gtggCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTgg                                        >  2:194572/1‑150 (MQ=255)                                                                                                             cacaCCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCg                                   >  1:251743/1‑150 (MQ=255)                                                                                                               cacCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTc                                                                                   >  2:934637/1‑100 (MQ=255)                                                                                                               cacCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTc                                                                                   <  1:934637/100‑1 (MQ=255)                                                                                                                acCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                                     >  2:273441/1‑129 (MQ=255)                                                                                                                acCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                                     <  1:273441/129‑1 (MQ=255)                                                                                                                acCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGcc                                >  1:784663/1‑150 (MQ=255)                                                                                                                     ggACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGc                           <  1:867095/150‑1 (MQ=255)                                                                                                                      gccaaccGAGCCCACACCCACTA‑TGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGcc                          <  2:594187/145‑1 (MQ=255)                                                                                                                       aCTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCt                         >  1:594187/1‑150 (MQ=255)                                                                                                                         ttCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGTACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTAc                       <  1:565958/150‑1 (MQ=255)                                                                                                                           ccGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACAt                     <  2:258525/150‑1 (MQ=255)                                                                                                                              tGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGa                  >  1:618793/1‑150 (MQ=255)                                                                                                                               gCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGaa                 <  2:618793/150‑1 (MQ=255)                                                                                                                                  cTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGc              >  2:1161257/1‑150 (MQ=255)                                                                                                                                  cTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGc              >  2:1129027/1‑150 (MQ=255)                                                                                                                                    cgcccgccA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACAt                                                                           >  2:361278/1‑87 (MQ=255)                                                                                                                                    cgcccgccA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACAt                                                                           <  1:361278/87‑1 (MQ=255)                                                                                                                                     gCCGACCAGGGG‑CATG‑CCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGAACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAAc           >  2:113674/1‑150 (MQ=255)                                                                                                                                     gCCGACCA‑GGGTCATG‑CCCTGGTGGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAAc           >  2:549430/1‑150 (MQ=255)                                                                                                                                        cacca‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGc        <  1:165093/150‑1 (MQ=255)                                                                                                                                          cca‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCt      <  2:784663/150‑1 (MQ=255)                                                                                                                                           ca‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTg     <  1:1117248/150‑1 (MQ=255)                                                                                                                                               gg‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGAc  >  2:693626/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                AACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGAC  >  NC_002947/4054095‑4054385

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