Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,332	C→A	8.0%	R17R (CGG→CGT)	trmE ←	GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,332	0	C	A	8.0%	193.6 / 2.6	61	R17R (CGG→CGT)	trmE	GTPase
Reads supporting (aligned to +/- strand):  ref base C (24/32);  new base A (2/3);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAG  >  NC_002947/6184‑6481                                                                                                                                                      |                                                                                                                                                        aCAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGcccc                                                                                                                                                        >  2:322473/1‑150 (MQ=255)  cAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGCCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCCGACAACCTGACAATTCCCACACCGCCCCg                                                                                                                                                       >  2:191935/1‑150 (MQ=255)    gagCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGc                                                                                                                                                     >  2:1061592/1‑150 (MQ=255)        aaTCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGTAATGGACCGGAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGCGATCCACA‑GCCGGCCTGGTCCAGCTAGACGTTCTGATTACCTGACAATTCCCACACCGCCCCGGCCC‑Tg                                                                                                                                                >  2:1592469/1‑150 (MQ=255)          taccTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGACGGCCTTG‑GCAGCCAGAGGGCCGGATAACCTGACCATTCCCACACCGCCCCGGCCC‑Tggg                                                                                                                                              >  2:204140/3‑150 (MQ=255)          tCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑Tggg                                                                                                                                              <  2:918245/150‑1 (MQ=255)           cccTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGc                                                                                                                                             <  2:271800/150‑1 (MQ=255)            ccTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGcc                                                                                                                                            >  2:1216973/1‑150 (MQ=255)               tCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTg                                                                                                                                         <  2:1571008/150‑1 (MQ=255)                 gTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGc                                                                                                                                       <  2:672023/150‑1 (MQ=255)                      gCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATgg                                                                                                                                  >  1:1110368/1‑150 (MQ=255)                       cACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCAGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCC‑AGGGCCGGGGCGGTgtg                                                                                                                                 <  2:85188/150‑3 (MQ=255)                        aCCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACATCGCCACGGCCCAGGGGCCGGGGCggtgtg                                                                                                                                 <  2:660909/150‑6 (MQ=255)                                  cATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGttt                                                                                                                      <  2:1110368/150‑1 (MQ=255)                                       tcgCGGAATGGATCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGc                                                                                                                 >  2:828428/1‑150 (MQ=255)                                           ggAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGTCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGc                                                                                                                                             >  2:1408493/1‑118 (MQ=255)                                           ggAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGTCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGc                                                                                                                                             <  1:1408493/118‑1 (MQ=255)                                                     gTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGtt                                                                                                         >  1:397872/1‑144 (MQ=255)                                                     gTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGtt                                                                                                         <  2:397872/144‑1 (MQ=255)                                                      tAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                  <  1:240670/150‑1 (MQ=255)                                                      tAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                  <  1:1216973/150‑1 (MQ=255)                                                            gCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                        <  2:974361/138‑1 (MQ=255)                                                            gCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                        >  1:974361/1‑138 (MQ=255)                                                            gCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                            <  2:1554029/150‑1 (MQ=255)                                                            gCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                            >  1:649861/1‑150 (MQ=255)                                                               tGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACACCCTGACAAATACCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGTCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCctc                                                                                         <  1:524011/150‑1 (MQ=255)                                                                  tgcgGGGTGAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                      <  1:1096010/149‑1 (MQ=255)                                                                  cgcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                      <  2:1539124/150‑1 (MQ=255)                                                                  cgcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                      <  2:1205701/150‑1 (MQ=255)                                                                   gcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACg                                                                                     >  2:1060617/1‑150 (MQ=255)                                                                              gTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGAc                                                                          <  1:662671/150‑1 (MQ=255)                                                                                             aaCAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGc                                                           >  2:1153201/1‑150 (MQ=255)                                                                                             aaCAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGc                                                           >  2:1273055/1‑150 (MQ=255)                                                                                             aaCAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGc                                                           >  1:1173728/1‑150 (MQ=255)                                                                                              aCAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                          >  1:996866/1‑150 (MQ=255)                                                                                              aCAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                          >  1:335093/1‑150 (MQ=255)                                                                                                  gccggccTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                        <  2:591483/100‑1 (MQ=255)                                                                                                  gccggccTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                        >  1:591483/1‑100 (MQ=255)                                                                                                         ttG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑gggg                                              >  2:592273/1‑150 (MQ=255)                                                                                                             gCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCC‑CTAGTGGGGCGt                                           <  1:381709/150‑1 (MQ=255)                                                                                                                gCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGtttt                                        <  2:649861/150‑1 (MQ=255)                                                                                                                   agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GTGGCGTTTTTTa                                     <  1:1060617/150‑1 (MQ=255)                                                                                                                   agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTa                                     <  1:322473/150‑1 (MQ=255)                                                                                                                    gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTAt                                    <  1:1061592/150‑1 (MQ=255)                                                                                                                          cGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAg                              <  1:98430/150‑1 (MQ=255)                                                                                                                           ggACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑g                                                 >  2:49992/1‑130 (MQ=255)                                                                                                                           ggACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑g                                                 <  1:49992/130‑1 (MQ=255)                                                                                                                                aCCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGc                        <  2:996866/150‑1 (MQ=255)                                                                                                                                   tGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTa                     <  2:1173728/150‑1 (MQ=255)                                                                                                                                    gACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                              <  2:723282/61‑1 (MQ=255)                                                                                                                                    gACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                              >  1:723282/1‑61 (MQ=255)                                                                                                                                    gACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTAt                    >  2:243027/1‑150 (MQ=255)                                                                                                                                      cAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTa                  <  1:828428/150‑1 (MQ=255)                                                                                                                                      cAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTa                  <  1:1592469/150‑1 (MQ=255)                                                                                                                                        aTTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGc                <  1:592273/150‑1 (MQ=255)                                                                                                                                             cacaCCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc           <  1:243027/150‑1 (MQ=255)                                                                                                                                              acacCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCa          >  1:123075/1‑150 (MQ=255)                                                                                                                                                atcGCCACGGCCC‑AGGGCCGGGGCGGTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag        >  2:575248/3‑150 (MQ=255)                                                                                                                                                atcGCCACGGCCC‑AGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTa                                     <  2:84631/119‑1 (MQ=255)                                                                                                                                                atcGCCACGGCCC‑AGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGACGGTTTTTa                                     >  1:84631/3‑121 (MQ=255)                                                                                                                                                 ccGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                            >  1:1187961/1‑66 (MQ=255)                                                                                                                                                 ccGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                            <  2:1187961/66‑1 (MQ=255)                                                                                                                                                      ccGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagcag  >  1:428037/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                        ACAGAGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAG  >  NC_002947/6184‑6481

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