Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,760,265	C→A	6.2%	D94E (GAC→GAA)	pvdM →	dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,760,265	0	C	A	6.2%	206.1 / 2.4	64	D94E (GAC→GAA)	pvdM	dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base A (2/2);  total (32/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 = 7.77e-01

TTCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTGCGCTACAACGC  >  NC_002947/4760116‑4760414                                                                                                                                                      |                                                                                                                                                      ttCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGAc                                                                                                                                                       >  2:1729623/1‑150 (MQ=255)    gACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTg                                                                                                                                                    >  2:1579252/1‑150 (MQ=255)    gACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTg                                                                                                                                                    <  1:1647949/150‑1 (MQ=255)      cAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGc                                                                                                                                                  >  2:1674415/1‑150 (MQ=255)         ggACGGCTACCCCATAGAGGTCGTGTGCTTGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACCGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCaa                                                                                                                                               >  1:1484546/1‑150 (MQ=255)             ggCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGCCCAAAGcc                                                                                                                                           >  1:1797857/1‑150 (MQ=255)              gCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGATGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcg                                                                                                                                          <  2:982777/150‑1 (MQ=255)                    ccATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTCACAACGATGGCCCGGGCCAGTTCGACCTCGCCACAGCcgcgcgcg                                                                                                                                    >  1:55695/1‑150 (MQ=255)                     cATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCgg                                                                                                                                   >  1:1018469/1‑150 (MQ=255)                     cATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCgg                                                                                                                                   >  1:1017911/1‑150 (MQ=255)                     cATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGGTCGAACTGGcccgg                                                                                                                                              <  2:1598045/139‑4 (MQ=255)                     cATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGGTCGAACTGGcccgg                                                                                                                                              >  1:1598045/1‑136 (MQ=255)                     cATAGAGGGCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCgg                                                                                                                                   <  1:1460953/150‑1 (MQ=255)                      aTAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCggg                                                                                                                                  <  1:1729623/150‑1 (MQ=255)                       tAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGc                                                                                                                                 <  1:1218787/150‑1 (MQ=255)                         gagGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGa                                                                                                                               <  1:1674415/150‑1 (MQ=255)                                    ccAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGgcgc                                                                                                                    >  1:623895/1‑150 (MQ=255)                                     cAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGt                                                                                                                           >  1:1332369/1‑142 (MQ=255)                                     cAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGt                                                                                                                           <  2:1332369/142‑1 (MQ=255)                                        gCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGcc                                                                                                                <  2:442604/150‑1 (MQ=255)                                         ccAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGATTTTTACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGccc                                                                                                               >  1:118504/1‑150 (MQ=255)                                          cAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTTGGAACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCt                                                                                                              <  2:58790/150‑1 (MQ=255)                                          cAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCt                                                                                                              <  2:616766/150‑1 (MQ=255)                                           aaTGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGa                                                                                                                               >  2:905212/1‑132 (MQ=255)                                           aaTGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGa                                                                                                                               <  1:905212/132‑1 (MQ=255)                                            aTGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGa                                                                                                            <  1:1579252/150‑1 (MQ=255)                                             tGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGCCTGTCCGGCGCTGCCCTGAc                                                                                                           >  1:819512/1‑150 (MQ=255)                                             tGATCTGCAAGAGCGGATCATTCCGTACTACAGCCATCCCATCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGAc                                                                                                           <  1:375/150‑1 (MQ=255)                                              gATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGAcc                                                                                                          <  1:1704840/150‑1 (MQ=255)                                                        agCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGCCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCtt                                                                                                <  2:1797857/150‑1 (MQ=255)                                                           ggATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGTGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGc                                                                                             >  1:422716/1‑150 (MQ=255)                                                               cATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGa                                                                                         >  1:736474/1‑150 (MQ=255)                                                                   tcgttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTcc                                                                                     <  1:1325184/150‑1 (MQ=255)                                                                    cgttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCt                                                                                    <  2:319677/150‑1 (MQ=255)                                                                     gttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGa                                                                                                            >  1:1165976/1‑125 (MQ=255)                                                                     gttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGa                                                                                                            <  2:1165976/125‑1 (MQ=255)                                                                     gttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTg                                                                                   <  1:174234/150‑1 (MQ=255)                                                                            aGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGc                                                                            <  2:1025180/150‑1 (MQ=255)                                                                              ccATCTCACCCTGCCGCTGGACTTCGGCACCGAAGCCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGACTTCGCTTGGCCCGAGTCCTGGACCGGCGa                                                                          >  2:1018914/1‑150 (MQ=255)                                                                                  ctcACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAAc                                                                      >  2:435144/1‑150 (MQ=255)                                                                                  ctcACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAAc                                                                      >  1:980566/1‑150 (MQ=255)                                                                                  ctcACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAAc                                                                      >  1:980565/1‑150 (MQ=255)                                                                                     aCCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCCGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGc                                                                   >  1:612275/1‑150 (MQ=255)                                                                                        cTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCg                                                                <  2:736474/150‑1 (MQ=255)                                                                                          gCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCCTTGGGCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGca                                                              <  2:643768/150‑1 (MQ=255)                                                                                           ccGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGcac                                                             >  1:1630177/1‑150 (MQ=255)                                                                                            cGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGTCTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGcacc                                                            >  2:1038120/1‑150 (MQ=255)                                                                                                      cGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGa                                                  <  1:1018914/150‑1 (MQ=255)                                                                                                           ccGAAGGCCCCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACgaaa                                             <  1:534039/150‑4 (MQ=255)                                                                                                            cGAAGGCCGCGAAGCTGACAAGGATGTCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGttt                                            >  1:521548/1‑150 (MQ=255)                                                                                                            cGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACg                                                <  2:1570359/146‑1 (MQ=255)                                                                                                            cGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACg                                                >  1:1570359/1‑146 (MQ=255)                                                                                                                 gCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTtcgtcg                                       >  1:632760/1‑150 (MQ=255)                                                                                                                  ccGCGCAGTTGACGTGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGa                                      <  2:612275/150‑1 (MQ=255)                                                                                                                         gCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCggg                               <  2:1331698/150‑1 (MQ=255)                                                                                                                             aCAAGGATGGCCCGGGCCAGTTCGACCTGGCCTAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGAc                           >  2:1068209/1‑150 (MQ=255)                                                                                                                             aCAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCGCGCGCCACCCCTCCGACGGTTTCGGCGGTGCCGGGCGAc                           >  2:576318/1‑150 (MQ=255)                                                                                                                                 ggATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGc                       <  2:623895/150‑1 (MQ=255)                                                                                                                                       cccGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGa                 >  2:1328313/1‑150 (MQ=255)                                                                                                                                           ggCCAGGTCGAACTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGccc                                                                                                               >  2:561764/1‑52 (MQ=38)                                                                                                                                           ggCCAGGTCGAACTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGccc                                                                                                               <  1:561764/52‑1 (MQ=38)                                                                                                                                            gCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTgc            <  2:1630177/150‑1 (MQ=255)                                                                                                                                            cccAATTCGACCTGGCCAACGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGTTGCCGGGCGACATGCCCAGGAAGTgc            <  1:667092/149‑1 (MQ=255)                                                                                                                                             ccAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTgcg           >  2:399567/1‑150 (MQ=255)                                                                                                                                                 ttCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACAGTCTTCGCTTGGCCAGAGTCCTGGACCGGCGAGAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGGGACATGGCCAGGAAGTGGGCTAc       >  1:646150/1‑150 (MQ=255)                                                                                                                                                      ccTGGCCAAAGCCGCGTGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACACCGGCCCGCCCCACCCCACCGACCGTTTCGTCGATGCCGGGCGCCATGCCCTGGAAGTGCGGAACACCGc  <  2:422716/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TTCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTGCGCTACAACGC  >  NC_002947/4760116‑4760414

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