Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,054,237:1	+T	6.3%	coding (767/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,237	1	.	T	6.3%	246.7 / 11.1	63	coding (767/1320 nt)	PP_3573	monooxygenase
Reads supporting (aligned to +/- strand):  ref base . (29/30);  new base T (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGAC  >  NC_002947/4054089‑4054385                                                                                                                                                       |                                                                                                                                                     tACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑c                                                                                                                                                     >  1:1175774/1‑150 (MQ=255)   ccGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CAt                                                                                                                                                   <  1:181479/150‑1 (MQ=255)       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:412896/1‑150 (MQ=255)                                                              ccAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAcc                                                                                        <  2:672464/150‑1 (MQ=255)                                                                aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTg                                                                                      <  2:663520/150‑1 (MQ=255)                                                                aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTg                                                                                      >  1:62562/1‑150 (MQ=255)                                                                aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACAGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCa                                                                                                <  2:1049612/140‑1 (MQ=255)                                                                aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACAGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCa                                                                                                >  1:1049612/1‑140 (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:62562/150‑1 (MQ=255)                                                                     cGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCga                                                                                 >  2:1117248/1‑150 (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                                                                     <  1:112661/150‑1 (MQ=255)                                                                                 tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcg                                                                     >  2:424106/1‑150 (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:1129027/1‑150 (MQ=255)                                                                                                                                        cTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGc              >  2:1161257/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)                                                                                                                                                       |                                                                                                                                                     TACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCA‑GGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGAC  >  NC_002947/4054089‑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: ‑