Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,688,218	Δ1 bp	5.8%	coding (292/1107 nt)	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,218	0	T	.	5.8%	210.5 / 6.9	69	coding (292/1107 nt)	PP_3253	carboxylate‑amine ligase
Reads supporting (aligned to +/- strand):  ref base T (36/29);  new base . (2/2);  total (38/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 = 6.69e-01

TCCGCAGCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGC  >  NC_002947/3688069‑3688364                                                                                                                                                       |                                                                                                                                                    tCCGCAGCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGt                                                                                                                                                     <  2:683638/150‑1 (MQ=255)    gCAGCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGgc                                                                                                                                                  >  1:607369/1‑150 (MQ=255)        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                                                                                                                                             <  1:1557423/131‑1 (MQ=255)                            cGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGctgctg                                                                                                                                             >  2:1557423/1‑131 (MQ=255)                            cGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAg                                                                                                                         <  2:1242015/150‑1 (MQ=255)                              gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGGTTGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAggg                                                                                                                       <  1:604547/150‑1 (MQ=255)                              gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGctgctg                                                                                                                                             >  1:1106303/1‑129 (MQ=255)                              gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGctgctg                                                                                                                                             <  2:1106303/129‑1 (MQ=255)                                   cACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACta                                                                                                                  >  1:1607459/1‑150 (MQ=255)                                      cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGTCGAAGAGGGGATGGGGCCATACTCTGCGGCAA‑GCCACCCCTGCGCAGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATaa                                                                                                               >  1:1754135/1‑150 (MQ=255)                                      cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACtat                                                                                                                 >  1:818224/1‑148 (MQ=255)                                      cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACtat                                                                                                                 <  2:818224/148‑1 (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                                                                                                                               <  1:103743/127‑1 (MQ=255)                                             tATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCG‑GGCTTCTGCAGCAGCCACGCgg                                                                                                                               >  2:103743/1‑127 (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                                                                                                   >  1:1769645/1‑144 (MQ=255)                                                        tGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcga                                                                                                   <  2:1769645/144‑1 (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                                                                                                           <  1:605525/63‑1 (MQ=255)                                                                                                                                 gCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                           >  2:605525/1‑63 (MQ=255)                                                                                                                                       ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACCAGGGCATTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCa               >  1:1883918/1‑150 (MQ=255)                                                                                                                                       ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagca                                                                                                            >  2:363477/1‑57 (MQ=255)                                                                                                                                       ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagca                                                                                                            <  1:363477/57‑1 (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                        >  2:121763/1‑136 (MQ=255)                                                                                                                                            ccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGa                        <  1:121763/136‑1 (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)                                                                                                                                                       |                                                                                                                                                    TCCGCAGCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGC  >  NC_002947/3688069‑3688364

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