Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,688,223	G→T	6.5%	L99L (CTG→CTT)	PP_3253 →	carboxylate‑amine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,688,223	0	G	T	6.5%	192.4 / 2.4	61	L99L (CTG→CTT)	PP_3253	carboxylate‑amine ligase
Reads supporting (aligned to +/- strand):  ref base G (33/24);  new base T (2/2);  total (35/26)
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.63e-01

CCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCC  >  NC_002947/3688076‑3688372                                                                                                                                                     |                                                                                                                                                       ccAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGgctgct                                                                                                                                                      <  1:1052095/150‑1 (MQ=255) ccAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGgctgct                                                                                                                                                      <  1:1690977/150‑1 (MQ=255)  cAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGctgctg                                                                                                                                                     <  2:1821175/150‑1 (MQ=255)       cGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGaa                                                                                                                                                <  1:1107967/150‑1 (MQ=255)         aGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGTCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGc                                                                                                                                              >  2:981157/1‑150 (MQ=255)         aGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGc                                                                                                                                              >  2:1284209/1‑150 (MQ=255)             ggCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCTTGTGCCGCG‑GGCTTCTGCAGCAGccacgc                                                                                                                                         <  2:1374182/150‑5 (MQ=255)              gCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑cgcg                                                                                                                                        <  2:607369/150‑1 (MQ=255)                     cGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGctgctg                                                                                                                                                     >  2:1557423/1‑131 (MQ=255)                     cGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGctgctg                                                                                                                                                     <  1:1557423/131‑1 (MQ=255)                     cGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAg                                                                                                                                 <  2:1242015/150‑1 (MQ=255)                       gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGGTTGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAggg                                                                                                                               <  1:604547/150‑1 (MQ=255)                       gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGctgctg                                                                                                                                                     <  2:1106303/129‑1 (MQ=255)                       gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGctgctg                                                                                                                                                     >  1:1106303/1‑129 (MQ=255)                            cACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACta                                                                                                                          >  1:1607459/1‑150 (MQ=255)                               cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGTCGAAGAGGGGATGGGGCCATACTCTGCGGCAA‑GCCACCCCTGCGCAGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATaa                                                                                                                       >  1:1754135/1‑150 (MQ=255)                               cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACtat                                                                                                                         <  2:818224/148‑1 (MQ=255)                               cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACtat                                                                                                                         >  1:818224/1‑148 (MQ=255)                               cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATaa                                                                                                                       >  1:1013538/1‑150 (MQ=255)                                 aCCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAag                                                                                                                      >  1:1795360/1‑149 (MQ=255)                                 aCCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAag                                                                                                                      <  2:1795360/149‑1 (MQ=255)                                   cTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                   <  2:1134571/150‑1 (MQ=255)                                    tGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCGTTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcagc                                                                                                                  >  1:48504/1‑150 (MQ=255)                                      tATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCG‑GGCTTCTGCAGCAGCCACGCgg                                                                                                                                       >  2:103743/1‑127 (MQ=255)                                      tATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCG‑GGCTTCTGCAGCAGCCACGCgg                                                                                                                                       <  1:103743/127‑1 (MQ=255)                                       aTGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                   >  2:918223/1‑146 (MQ=255)                                       aTGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                   <  1:918223/146‑1 (MQ=255)                                          aaGCACGTGAGTTTCTCCAGCGCAGCCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcg                                                                                                            >  1:599578/1‑150 (MQ=255)                                           aGCACGTGAGTTTCTCCAGCGCAACAGTAAGAGGCTGAGAGTGTCATTGGTCGAAGAGGGGATGGGGCCATACTCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcga                                                                                                           <  2:1754135/150‑1 (MQ=255)                                                gTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTAGACGACt                                                                                                      >  1:392814/1‑150 (MQ=255)                                                 tGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcga                                                                                                           <  2:1769645/144‑1 (MQ=255)                                                 tGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcga                                                                                                           >  1:1769645/1‑144 (MQ=255)                                                      ttCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGcc                                                                                                >  2:1756039/1‑150 (MQ=255)                                                      ttCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGcc                                                                                                >  2:224497/1‑150 (MQ=255)                                                        ctcCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCAt                                                                                              <  1:1658318/150‑1 (MQ=255)                                                         tcCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCAtg                                                                                             >  1:225009/1‑150 (MQ=255)                                                           cAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACGCCGTCTCTATAAGCAGCTGTTCGACGACTACCGCCAtgtg                                                                                           >  1:1039649/1‑150 (MQ=255)                                                             gcgcAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGc                                                                                         >  2:1371756/1‑150 (MQ=255)                                                               gcAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGGTTACGGCTG‑GCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCaa                                                                                       <  1:797812/150‑1 (MQ=255)                                                                         gCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                   <  2:1629104/112‑1 (MQ=255)                                                                         gCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                   >  1:1629104/1‑112 (MQ=255)                                                                            gCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCt                                                                          <  1:1011127/150‑1 (MQ=255)                                                                              tGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCa                                                                        >  1:735168/1‑150 (MQ=255)                                                                                   gTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGc                                                                   >  2:378501/1‑150 (MQ=255)                                                                                      gCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTg                                                                >  2:300294/1‑150 (MQ=255)                                                                                         ttGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCAt                                                             <  2:735168/150‑1 (MQ=255)                                                                                               gAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAg                                                                                      >  1:1544947/1‑119 (MQ=255)                                                                                               gAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAg                                                                                      <  2:1544947/119‑1 (MQ=255)                                                                                                          ggggCCATACGCTGCGGCAA‑GCCACCGCGGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgcc                                            <  2:439797/150‑1 (MQ=255)                                                                                                          ggggCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgcc                                            >  2:176104/1‑150 (MQ=255)                                                                                                            ggCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGc                                                                               <  2:785029/113‑1 (MQ=255)                                                                                                            ggCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGc                                                                               >  1:785029/1‑113 (MQ=255)                                                                                                                 tACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGTACTATAAGCAGCCGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCt                                     >  2:373455/1‑150 (MQ=255)                                                                                                                          gCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                   >  2:605525/1‑63 (MQ=255)                                                                                                                          gCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                   <  1:605525/63‑1 (MQ=255)                                                                                                                                ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACCAGGGCATTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCa                       >  1:1883918/1‑150 (MQ=255)                                                                                                                                ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagca                                                                                                                    <  1:363477/57‑1 (MQ=255)                                                                                                                                ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagca                                                                                                                    >  2:363477/1‑57 (MQ=255)                                                                                                                                ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCa                       >  1:502776/1‑150 (MQ=255)                                                                                                                                 cACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCaa                      >  2:1863199/1‑150 (MQ=255)                                                                                                                                 cACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCaa                      <  2:48504/150‑1 (MQ=255)                                                                                                                                  aCCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAAc                     >  1:259837/1‑150 (MQ=255)                                                                                                                                    ccttgtgccgcG‑GGCTTCTGCAGCAGCCACGCGGCGCAGGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTg                   >  2:1673461/7‑150 (MQ=37)                                                                                                                                     ccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGa                                <  1:121763/136‑1 (MQ=255)                                                                                                                                     ccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGa                                >  2:121763/1‑136 (MQ=255)                                                                                                                                     ccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGa                  <  1:247529/150‑1 (MQ=255)                                                                                                                                          gccgcGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAAc             >  2:686074/1‑150 (MQ=255)                                                                                                                                             gcgTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGc          <  1:300294/150‑1 (MQ=255)                                                                                                                                                   ctgctgCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTg    >  2:575005/1‑150 (MQ=255)                                                                                                                                                     gctgcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGcc  <  2:1619407/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       CCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCC  >  NC_002947/3688076‑3688372

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