Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,173,052	T→A	7.6%	Q165H (CAA→CAT) ‡	PP_1925 ←	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,173,052	0	T	A	7.6%	278.6 / 11.3	93	Q165H (CAA→CAT) ‡	PP_1925	monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (42/43);  new base A (4/3);  total (47/46)
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.36e-01

CCGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCC  >  NC_002947/2172905‑2173193                                                                                                                                                    |                                                                                                                                              ccGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTg                                                                                                                                             <  1:849692/150‑1 (MQ=255)   gTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGa                                                                                                                                           >  2:472727/1‑150 (MQ=255)       gAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCt                                                                                                                                               >  2:1282903/1‑142 (MQ=255)       gAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCt                                                                                                                                               <  1:1282903/142‑1 (MQ=255)       gAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGc                                                                                                                                       >  1:1627135/1‑150 (MQ=255)       gAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTAACCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCAAGGAAAGc                                                                                                                                       >  1:806792/1‑150 (MQ=255)        aaCAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGcc                                                                                                                                      >  2:1753204/1‑150 (MQ=255)         aCAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCAAGGAAAGcgc                                                                                                                                     <  2:806792/150‑3 (MQ=255)             cacaCGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCggg                                                                                                                                 >  2:1495653/1‑150 (MQ=255)                 cGGCCATAAACGTCATCAGGCATCACCGTCGGTTTCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGtgat                                                                                                                             <  2:587366/150‑1 (MQ=255)                  ggCCATCAACGTCATCAGGCAGAAACGCCGGTTCCGCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGtgatt                                                                                                                            <  2:123857/150‑1 (MQ=255)                     cATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTGCGCATAGATTTGCGCCCCCGAGTTGCCACAACCAACGACCAGCACGCGCTTTCCTTGGAAACCCTCGGGtgattggt                                                                                                                         >  2:518201/1‑148 (MQ=255)                     cATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGtgattgat                                                                                                                         >  2:866493/1‑150 (MQ=255)                     cATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGtgattgat                                                                                                                         >  2:1442107/1‑150 (MQ=255)                      aTCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATa                                                                                                                        >  2:11253/1‑150 (MQ=255)                       tCAACGTCATCAGGCATAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGACaa                                                                                                                       >  2:610968/1‑150 (MQ=255)                        cAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAAt                                                                                                                      <  1:1495653/150‑1 (MQ=255)                        cAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGAGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAAt                                                                                                                      >  1:1246291/1‑150 (MQ=255)                         aaCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATg                                                                                                                     <  2:1736169/150‑1 (MQ=255)                         aaCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATg                                                                                                                     >  1:1450229/1‑150 (MQ=255)                         aaCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTAGGACAGCTCAGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGCAAAGGACGCGCTTTCCTTGGAAAGACTCGGGTGATTGACAATg                                                                                                                     >  1:1451606/1‑150 (MQ=255)                          aCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATgg                                                                                                                    >  2:1170385/1‑150 (MQ=255)                          aCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATgg                                                                                                                    >  2:853849/1‑150 (MQ=255)                          aCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCAAGGAAAGcgc                                                                                                                                     <  2:1742373/133‑3 (MQ=255)                          aCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCAAGGAAAGcgc                                                                                                                                     <  2:1742329/133‑3 (MQ=255)                          aCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCAAGGAAAGcgc                                                                                                                                     >  1:1742329/1‑131 (MQ=255)                          aCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCAAGGAAAGcgc                                                                                                                                     >  1:1742373/1‑131 (MQ=255)                               atcaGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGa                                                                                                               <  2:535270/150‑1 (MQ=255)                                   ggCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGa                                                                                                           <  1:853849/150‑1 (MQ=255)                                   ggCAGAAACGCCGGTCCCCCCGTGGTTACCCACCTTGCGTCGCCGCCCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCACCGAACGGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGa                                                                                                           <  1:1078962/150‑1 (MQ=255)                                    gCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGaa                                                                                                          <  1:1175206/150‑1 (MQ=255)                                      aGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATc                                                                                                        >  1:1194580/1‑150 (MQ=255)                                        aaaCGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATctg                                                                                                      <  2:1526496/150‑1 (MQ=255)                                            gCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGGTAATGGGCGGAGTGAATCTGCTGt                                                                                                  >  1:1310866/1‑150 (MQ=255)                                               ggTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACAAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTccc                                                                                               <  2:640052/150‑1 (MQ=255)                                                gTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCt                                                                                              >  2:753285/1‑150 (MQ=255)                                                   cctccGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGa                                                                                           <  2:1627135/150‑1 (MQ=255)                                                    ctccGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGaa                                                                                          <  2:1450229/150‑1 (MQ=255)                                                       cGTGGTTAACCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAATGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAg                                                                                       <  2:1041534/150‑1 (MQ=255)                                                        gTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGAGTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGc                                                                                      >  1:1587224/1‑150 (MQ=255)                                                         tGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGct                                                                                     <  2:1119505/150‑1 (MQ=255)                                                               cccACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATctg                                                                                                      <  2:1700154/127‑1 (MQ=255)                                                               cccACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATctg                                                                                                      >  1:1700154/1‑127 (MQ=255)                                                                     ttGCGTCGGAGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAAt                                                                         <  2:1246291/150‑1 (MQ=255)                                                                       gCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAg                                                                       <  1:472727/150‑1 (MQ=255)                                                                       gCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAg                                                                       <  1:836987/150‑1 (MQ=255)                                                                         gTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGACAGCCTCGGGTGATTGATACTGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGc                                                                     <  2:1857199/150‑1 (MQ=255)                                                                           cGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCg                                                                   <  1:753285/150‑1 (MQ=255)                                                                               gACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAAc                                                               >  2:917765/1‑150 (MQ=255)                                                                               gACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCCGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAAc                                                               <  2:1093120/150‑1 (MQ=255)                                                                                  cTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATCATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATa                                                            >  2:879035/1‑150 (MQ=255)                                                                                  cTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATa                                                            >  1:1280789/1‑150 (MQ=255)                                                                                  cTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATa                                                            >  1:302406/1‑150 (MQ=255)                                                                                  cTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATa                                                            >  2:1388945/1‑150 (MQ=255)                                                                                   tCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAg                                                           >  1:1604369/1‑150 (MQ=255)                                                                                   tCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAg                                                           <  1:866493/150‑1 (MQ=255)                                                                                   tCGGACAGCTCCGCATAGATTTGCGCCCACGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAg                                                           >  2:1664519/1‑150 (MQ=255)                                                                                     ggACAGCTCCGCATAGATTTGCGCCCACGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAggg                                                         <  1:1664519/150‑1 (MQ=255)                                                                                      gACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAgggg                                                        <  1:31644/150‑1 (MQ=255)                                                                                       aCAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCAACACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGt                                                       <  2:268667/150‑1 (MQ=255)                                                                                          gCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGc                                                    <  1:970805/150‑1 (MQ=255)                                                                                           cTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAggg                                                         >  1:1235155/1‑144 (MQ=255)                                                                                           cTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAggg                                                         <  2:1235155/144‑1 (MQ=255)                                                                                               gCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGGTAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAg                                               <  2:1310866/150‑1 (MQ=255)                                                                                               gCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAg                                               <  2:139965/150‑1 (MQ=255)                                                                                               gCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCGTGGACAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTCGAACAGCTCTATCCCCGAATAGGCCGGGCCATAGGGGttcccccgg                                               >  2:231122/1‑143 (MQ=255)                                                                                                cATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAgg                                              >  1:645686/1‑150 (MQ=255)                                                                                                 aTAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCg                                                                   <  2:1215343/128‑1 (MQ=255)                                                                                                 aTAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCg                                                                   >  1:1215343/1‑128 (MQ=255)                                                                                                  tAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTAt                                                                                 <  2:825490/113‑1 (MQ=255)                                                                                                  tAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTAt                                                                                 >  1:825490/1‑113 (MQ=255)                                                                                                           gcCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTgcg                                   >  2:1861952/1‑150 (MQ=255)                                                                                                             ccccGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCt                                 >  2:286734/1‑150 (MQ=255)                                                                                                             ccccGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTACCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCt                                 <  1:1388945/150‑1 (MQ=255)                                                                                                               ccGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCTGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCtga                               >  1:350577/1‑150 (MQ=255)                                                                                                                cGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTAGCTTGATGGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCtgat                              >  1:1788524/1‑150 (MQ=255)                                                                                                                      gCCACCACCAACGCCCAGCACGCGCTCTGCTTGGAAAGCCTCGGGTGATTGACAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATgacgtt                        <  2:1451606/150‑3 (MQ=255)                                                                                                                         accaccaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAggg                                                         <  2:1919769/114‑1 (MQ=255)                                                                                                                         accaccaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAggg                                                         >  1:1919769/1‑114 (MQ=255)                                                                                                                         accaccaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACg                     >  2:618603/1‑150 (MQ=255)                                                                                                                           caccaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAg                   >  2:1343106/1‑150 (MQ=255)                                                                                                                           caccaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAg                   >  1:298577/1‑150 (MQ=255)                                                                                                                            accaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGc                  <  2:302406/150‑1 (MQ=255)                                                                                                                              caACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAgg                                                          <  1:1052835/108‑1 (MQ=255)                                                                                                                              caACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTACTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAgg                                                          >  2:1052835/1‑108 (MQ=255)                                                                                                                                 cGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTccc             <  2:1280789/150‑1 (MQ=255)                                                                                                                                 cGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTccc             <  2:1486083/150‑1 (MQ=255)                                                                                                                                  gACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCa            >  2:378234/1‑150 (MQ=255)                                                                                                                                    ccAGCAGGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGt          >  1:643630/1‑150 (MQ=255)                                                                                                                                          gaggcTTTCCAAGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGCCTGCTCCAGGTCCCCGGTGCGCTGATGACTACACGAGCCTCCCAGTCGCGGt    >  2:860123/4‑150 (MQ=255)                                                                                                                                           gcgcTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATctg                                                                                                      <  1:1548717/51‑1 (MQ=255)                                                                                                                                           gcgcTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATctg                                                                                                      >  2:1548717/1‑51 (MQ=255)                                                                                                                                           gcgcTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTc   <  1:1442107/150‑1 (MQ=255)                                                                                                                                           gcgcTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTc   <  1:906600/150‑1 (MQ=255)                                                                                                                                            cgcTTTCCTTGGAATGCCTCGGGTGATTGATAATGGGCGGTGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTcc  <  1:1843375/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              CCGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCC  >  NC_002947/2172905‑2173193

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